





PHYSICAL REVIEW 
LETTERS 












VoLuME 4 MARCH 15, 1960 NuMBER 6 














SOME PROPERTIES OF THE VAN ALLEN RADIATION 


A. J. Dessler and Robert Karplus* 
Lockheed Aircraft Corporation, Missiles and Space Division, Palo Alto, California 
(Received February 23, 1960) 


The Van Allen radiation belt (energetic particle 
radiation trapped in the geomagnetic field) has 
been shown!~* to consist of a hard proton com- 
ponent centered at about 10x10° km from the 
earth’s magnetic axis (the proton belt) and an 
electron component centered at about 22 x 10° km 
from the earth’s magnetic axis (the electron belt). 
The electron belt extends through the region 
occupied by the proton belt. It is generally ac- 
cepted that the proton belt is produced by the 
decay of cosmic-ray albedo neutrons, ** although 
some doubt has been expressed on this point.° 
The electron belt is usually held to be of solar 
origin.° The purposes of this Letter are two- 
fold: first, to point out that the results of the 
observations of the outer zone of the Van Allen 
radiation belt made with the Explorer IV and 
Explorer VI satellite systems are inconsistent 
with the solar injection hypothesis; and second, 
to show that the electrons released in the decay 
of cosmic-ray neutron albedo may represent a 
satisfactory source for the outer zone. Source 
strength and trapping lifetime are not discussed." 

An illustrative sketch of particle flux vs radial 
distance in the equatorial plane is shown in Fig.1. 
This sketch is based principally on the measure- 
ments made by Simpson and his co-workers.” 
Hydromagnetic scattering probably limits the 
radial extent of the proton belt."* The measure- 
ments show that the electron belt consists of two 
regions separated by a relative minimum near 
19x10° km from the geomagnetic axis. We know 
of no satisfactory explanation for this gap in the 
electron belt other than the one offered by the 






















Capetown anomaly (a region several thousand kilo- 
meters in extent, centered southwest of Capetown, 
South Africa, and having an abnormally weak mag- 
netic field strength). 

It has been pointed out that the Capetown anomaly 
lowers the mirror altitude of the trapped radia- 
tion reflected near the anomaly by about 1000 km.™ 
Since the atmospheric density is almost constant 
above about 1300 km where hydrogen predom- 
inates and increases very rapidly below this alti- 
tude where heavier atmospheric constituents pre- 
dominate, *® we should expect radiation mirroring 
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FIG. 1. Illustrative sketch of the Van Allen radia- 
tion belt particle flux vs radial distance in the equa- 
torial plane. This sketch is based principally on data 
presented by Simpson et al.'* The electron belt under- 
goes changes during magnetic storms which produce 
large variations in the counting rate of instruments 
flown through the belt. 
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below about 1300 km to be rapidly removed from 
the radiation belt. As the trapped radiation drifts 
around the earth’s magnetic axis, a zone of re- 
latively low radiation intensity will be developed 
at the magnetic latitude of the Capetown anomaly 
because that part of the Van Allen radiation 
which normally mirrors at altitudes up to 2300 
km will be removed by atmospheric scattering 

at the lower mirror altitude over the Capetown 
anomaly. The part of the radiation belt affected 
by the anomaly has been determined from the 
surfaces given by the two adiabatic invariants 
(magnetic moment and integral invariant) of a 
charged particle in the earth’s field as calculated 
by Vestine and Sibley** using the 48-term har- 
monic analysis of the geomagnetic field. With 
these data plus data on the altitudes of surfaces 
of constant field strength, '* the effect of the 
Capetown anomaly is illustrated in Fig. 2. In 
this figure, mirror altitude is plotted against 
longitude for 3 different shells (the surface gen- 
erated by a particle as it drifts arount the mag- 
netic axis) as defined by the adiabatic invariants. 
It is apparent from this figure that the greatest 
decrease from a “normal” mirror altitude occurs 
over the Capetown anomaly for the field line or 
shell which crosses the equatorial plane at 

17x 10° km from the earth’s center. Adjacent 
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FIG. 2. The altitude of a given mirroring field, 
By» along 3 different traces determined from integral 
invariant calculations for magnetic field lines which 
extend a distance 5 from the earth’s center in the 
equatorial plane. 

Trace (1): By, = 0.18 gauss, b= 13 10° km. 

Trace (2): B,, = 0.20 gauss, b=17 x 10° km. 

Trace (3): By, = 0.30 gauss, b= 28 x 10° km. 

Note that the average mirror point altitude for trace 
(2) in the northern hemisphere is lowered by about 
1000 km over the Capetown anomaly. Traces (1) and 
(3) are not so severely affected. 
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shells on either side of this one show a less 
pronounced effect. Further, it is expected that 
the sharpness of the anomaly is underestimated 
because the spherical harmonic analysis contains 
too few terms to describe its features adequately," 
This view is supported by recent magnetic field 
measurements made by Heppner et al., with the 
Vanguard III satellite.‘7 These measurements 
show that at altitudes of 2000 - 3000 km, the 
Capetown anomaly is much sharper and slightly 
deeper than predicted. Thus, we may conclude 
that the effect of the Capetown anomaly is much 
more localized than shown in Fig. 2. It therefore 
appears that the observed gap in the electron 
belt is in satisfactory agreement with the location 
of the Capetown anomaly. 

As we have said, the Capetown anomaly can 
lead to the rapid removal only of electrons with 
mirror points below 2300 km altitude which will 
be reflected below 1300 km over the anomaly. 
Electrons with high mirror points will not be 
affected and will be responsible for a counting 
rate in the region where the gap is observed. 
The electron source, therefore, if reasonably 
isotropic, must inject electrons predominantly 
at low altitudes so as to assure low mirror 
points. This is quite a severe requirement. It 
seems obvious that particles of solar origin 
would not satisfy it. Electrons from neutron 
decay, on the other hand, do have the desired 
property of mirroring principally at low altitudes. 
Hess® finds that the neutron decay density near 
the earth varies approximately as r~*(1+3 sin’), 
where is the latitude. This injection distribu- 
tion leads to a flux decrease of about 10% near 
the equator if we use an altitude of 1300 km as 
the lower cutoff for injection but retain only 
particles mirroring above 2300 km. On the basis 
of this model we predict that the depth of the gap 
in the electron belt will show an altitude depend- 
ence: the gap will be more pronounced at low 
altitudes than near the equator on the magnetic 
field shell that passes through the Capetown 
anomaly. 

The spectrum of the electrons resulting from 
neutron decay is peaked near 300 kev. The elec- 
trons which form the radiation belt are “cooled” 
by electron-electron collisions" to the observed 

E~* spectrum.” These interactions will also 
affect the depth of the gap as measured by de- 
tectors of different thresholds. 

We turn finally to a discussion of the magnetic- 
storm -induced fluctuations in the radiation belt.” 
These fluctuations (during which the counting 
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rate in the electron belt may change by 10*) are 
the principal evidence cited in support of the 
solar origin hypothesis of the trapped electrons. 
There are four difficulties with the solar injection 
hypothesis: (1) Solar particles could most prob- 
ably be injected during the active phase of a mag- 
netic storm, when the dipole character of the 
earth’s field is disturbed. Yet, at this time, a 
decrease in counting rate is observed. That is, 
when additional electrons could be injected, they 
are not. (2) After the active part of a storm is 
over (recovery phase—low K index) the observed 
counting rate increases by several orders of 
magnitude. However, this is the time when the 
geomagnetic field is restored so that the entry 

of any but cosmic-ray particles into the field is 
most unreasonable. Thus, we argue that the 
counting rate increase occurs too late for solar 
injection to have been effective. (3) The particles 
have a lifetime of a year or more in the very 
rare atmosphere at high altitude. This lifetime 
is in accord with the fact that the intensity of the 
belt is stable during magnetically quiet periods. 
Nevertheless, with the solar injection hypothesis 
it is proposed that additional electrons are in- 
jected during a magnetic storm and that most 

of them disappear with a lifetime of a few days. 
This short lifetime seems inadmissible. (4) As 
discussed earlier, the existence of the gap in 

the electron component, and especially its per- 
sistence during and after magnetic storms, seem 
impossible to reconcile with the solar injection 
hypothesis. 

The 10-kev solar wind protons, whose injection 
into the earth’s field gives rise to the stresses 
that produce the magnetic field decrease during 
the main phase of a magnetic storm, have a 
trapping lifetime that is comparable to the time 
required for the enhanced counting rate to return 
to prestorm values. These protons are removed 
by charge exchange with the neutral hydrogen 
telluric corona.’® We propose, therefore, that 
it is the magnetic field changes rather than 
changes in the number of trapped particles that 
give rise to the varying counting rates. There 
are two mechanisms by which this may happen. 
First, since the field changes are slow compared 
to the electron cyclotron frequency, the dis- 
placement of field lines will cause an accompany- 
ing displacement of particles. Many particles 
may therefore appear where only few had been 
before. Second, changes in the field intensity 
will result in an energy change of the particles; 
because the radiation counters have a threshold, 





the measured flux will be strongly affected. Both 
mechanisms can cause large changes in the 
counting rate, the first because the particle 
density varies rapidly with position, the second 
because the energy spectrum is very steep (e.g., . 
an integral spectrum? as steep as E~®). It seems 
clear that an appreciable fraction of the mag- 
netic storm energy can be temporarily given to 
the radiation belt by the mutual inductance be- 
tween the ring currents: (1) the solar wind protons 
and (2) the Van Allen belt electrons. 

There are some incompletely treated problems 
that must be resolved before we can claim a 
complete understanding of the radiation belt. 

The source strength, the trapping lifetime, and 
the details of the magnetic storm effects” are 

the most important. Still, we feel that the pre- 
sent analysis leads to a fairly satisfactory pic- 
ture in which the entire Van Allen radiation belt 

is due to neutron decay: The decay of fast neutrons 
produces the proton belt while the electron belt 

is attributed to the large number of slow neutrons 
which decay relatively near the earth. 

We wish to acknowledge the generous assistance 
of Dr. E. H. Vestine and Dr. W. N. Hess in making 
the results of their respective calculations avail- 
able to us prior to publication. 
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The 14.4-kev y ray emitted without recoil by 
0.1-ysec Fe” in metallic iron'* excited great 
interest as the most precisely defined electro- 
magnetic frequency yet discovered. It may be 
adequately well defined to allow measurement of 
the influence of a gravitational potential on fre- 
quency® and of other small effects hitherto be- 
yond the sensitivity available in the laboratory. 
As a preliminary step in the operation of an ex- 
perimental system designed to measure the 
gravitational effect, we have been making tests 
to find out whether other influences than the one 
intended might lead to systematic errors by in- 
troducing important frequency shifts not taken 
into account. 

So far the largest such effect found is that of 
temperature. That temperature should influence 
the frequency exactly as we observe is very 
simply explained. Thermally excited vibrations 
cause little broadening through first order 
Doppler effect under the conditions obtaining in 
the solid because the value of any component of 
the nuclear velocity averages very nearly to zero 
over the nuclear lifetime. The precision of the 
y ray of Fe*” requires the second order Doppler 
effect also to be considered. A shift to lower 
frequency with increased temperature results 
from this because the also well-defined average 
of the square of the velocity of the particle in- 
creases in direct proportion to the average kine- 
tic energy. As a consequence one would expect 
a temperature coefficient of frequency in a 
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VARIATION WITH TEMPERATURE OF THE ENERGY OF RECOIL- FREE 
GAMMA RAYS FROM SOLIDS* 


R. V. Pound and G. A. Rebka, Jr. 
Harvard University, Cambridge, Massachusetts 
(Received February 17, 1960) 


homogeneous solid, 


(8v/8T) =- vC ,/2Mc?, 


where C7, is the specific heat of the lattice and 
M is the gram atomic weight of iron. In the high- 
temperature classical limit where C; =3R, 


(8v/8T),_ = -2.44x10"7'5y per °K. 


At lower temperatures one would expect a coef- 
ficient reduced by the value of the appropriate 
normalized Debye specific heat function. For 
iron, at 300°K one should find about 0.9 times, 
and at 80°K about 0.3 times, the above classical 
value. 

The temperature dependence has been meas- 
ured by counting the y rays from our 0.4-curie 
Co™” source transmitted through enriched Fe” 
absorbing films (0.6 mg Fe*’/cm?). The Co™ of 
the source is distributed in about 3x10™*° cm 
thickness below the surface of a 2-in. diameter 
iron disk, made in the manner described ear- 
lier.’ Small frequency shifts that result when 
the source and absorber are held at different 
temperatures were measured by using a trans- 
ducer to move the source sinusoidally at ten cps 
toward and away from the absorber at a peak 
speed of about 0.01 cm/sec. A gate pulse and 
mercury relays were used to make one counter 
record during 25 milliseconds of the modulation 
period symmetrically disposed about the time of 
maximum velocity toward the absorber. Another 
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counter recorded the corresponding counts with 
the source going away from the absorber. The 
difference of the counts in the two registers 
should be proportional to the relative frequency 
shift of the absorber and source for shifts small 
compared to the line width. Quantitative know- 
ledge of the parameters of the system that are 
involved in determining the constant of propor- 
tionality is rendered unnecessary by adding 
through a clock-driven hydraulic system a con- 
tinuous relative motion of 6.3 x10™* cm/sec di- 
rected oppositely during each of the two halves 
of the time for a given datum point. In this way 
the sensitivity to frequency shift originating in 
the Doppler effect is measured simultaneously 
with the shift sought. The algebraic sum of the 
counting rate differences for the two halves of 
the run are proportional to the shift and the dif- 
ference to the sensitivity. 

The shift at liquid nitrogen relative to room 
temperature is comparable to the line width and 
for that point the two counting rates were re- 
corded at a series of values of the sinusoidal 
modulation amplitude. From these a value of the 
shift and of the apparent line width could be ob- 
tained although difficulties of calibration under 
the conditions of operation have contributed 
strongly to the uncertainties. There is evidence 
that the line appears to broaden with such a tem- 
perature difference by perhaps a factor of 2.3 
which might be evidence that the hyperfine struc - 
ture splittings are temperature sensitive to 
some extent, as must be expected. 

The data are plotted in Fig. 1. A solid line re- 
presenting the effect expected with a Debye tem- 
perature of 420°K is also drawn. The agreement 
can be regarded as an experimental demonstra- 
tion of the second order Doppler effect using 
thermal velocities rather than a centrifuge. It 
might be remarked that crystalline anisotropy 
might make this source of high velocities useful 
for experiments to the end of detecting such spa- 
tial anisotropies as might accompany ether drift 
or an inertial frame. 

The temperature sensitivity at room tempera- 
ture [experimentally (-2.09 + 0.24) x107'5 per 
degree C, theoretically -2.21x10-'* per degree 
C]is highly relevant to the interpretation of data 
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FIG. 1. Fractional shift of energy of 14.4-kev 
gamma-ray absorption of Fe’ vs absolute tempera- 
ture of the metal. The solid line is derived from as- 
suming a Debye temperature of 420°K. 


from our experiment on the effect of gravitational 
potential. A temperature difference of 1°C be- 
tween the top and the bottom of our 22- meter 
tower would result in a shift about equal to that 
predicted by the principle of equivalence. For 
smaller height differences correspondingly 
smaller temperature differences would be re- 
quired. It is now clear that correction for or 
control of the temperature difference and perhaps 
other parameters must be included in the instru- 
mentation of experiments intending to utilize the 
extreme frequency discrimination available with 
gamma rays of this type. 
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ELECTRON SPIN RESONANCE IN BISMUTH AND ANTIMONY 


G. E. Smith, J. K. Galt, and F. R. Merritt 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received February 17, 1960) 


Spin resonance of conduction electrons has 
been observed previously in the alkali metals’ 
with a g factor approximately that of the free 
electron, g=-2. More recently”»*, observations 
and calculations of the g factor for conduction 
electrons in the semiconductor InSb have shown 
it to be very large (g= -50) because of the small 
energy gap and large spin-orbit coupling involved. 
The same situation— small energy gap and large 
spin-orbit coupling— exists in the semimetals 
bismuth and antimony and, for the case of bis- 
muth, Cohen and Blount* have calculated the g 
factor in terms of the effective mass tensor. It 
is expected that their theory should be equally 
valid for antimony. This Letter reports the ex- 
perimental observations of spin resonance in 
these materials and generally confirms the pre- 
dictions of Cohen and Blount. 

The measurements were made on crystallo- 
graphically oriented plane surfaces of zone- 
refined single-crystal samples which had been 
acid polished. The experimental arrangement 
was such that the sample formed one end of a 
microwave cavity. A signal proportional to the 
variations in the Q of the cavity, and therefore 
to variations in the power absorption coefficient 
of the sample surface, was observed as a function 
of magnetic field. The measurements were made 
at 72 kMc/sec and 1.3°K. Details of the apparatus 
have been reported elsewhere. °® 

The results for antimony were obtained with 
the magnetic field perpendicular to the plane of 
the sample and with circularly polarized radi- 
ation incident upon it. The bisectrix axis of the 
crystal was oriented parallel to the magnetic 
field. A plot of power absorption vs magnetic 
field, shown in Fig. 1, exhibits a sharp reso- 
nance at 1.14 kilogauss and a sharp rise preceded 
by a dip and with some structure upon it at 2.28 
kilogauss. An upper limit on the line width of 
the sharp resonance at 1.14 kilogauss is 30 gauss. 
It is felt that this line could not be attributed to 
cyclotron resonance because its sharpness in- 
dicated too large a relaxation time. Furthermore, 
neither this line nor the dip and other structure 
on the sharp rise are understandable in terms 
of the theory of such a curve taking account of 
cyclotron resonance, while the sharp rise itself 


276 


is understandable in these terms as a “dielectric 
anomaly.”® 

The bismuth data shown in Fig. 2 were obtained 
using linearly polarized radiation with the static 
magnetic field in the plane of the sample and 
parallel to the trigonal axis. In the top curve, 
the microwave H field and the static H field are 
perpendicular, while for the bottom curve they 
are parallel to one another. The appearance of 
resonances in the top curve and not the bottom 
indicates a magnetic dipole transition and con- 
vinces us that the sharp resonances observed 
arise from spin flips. These effects have also 
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FIG. 1. Microwave absorption in antimony with 
circularly polarized radiation. The static field is 
normal to the sample surface and along a bisectrix 
axis. 
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FIG. 2. 
wave absorption in 
bismuth using plane 
polarized radiation. 
Curve A has the 
microwave H field 
parallel to the static 
H field and curve B 
has them perpendic- 
war to one another. 
The static field is 
along a trigonal 
axis. 
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been seen with the magnetic field along twofold 
and bisectrix axes in bismuth. 

The essential result of Cohen and Blount’s 
theory is that the spin effective mass msg [defined 
ty m,/m=2/g, where m is the free electron mass] 
should be equal to the orbital effective mass m, 
vhen the orbital effective mass is sufficiently 
small. An energy level diagram showing the 
sin splitting nearly equal to the orbital splitting 
isshown in Fig. 3. Higher Landau levels (n=3, 
{, etc.) are similar to the first two. Type 1 
tansitions are the usual cyclotron resonances, 
ad type 2 are pure spin resonances. Types 3 
ad 4 combine an orbital transition with a spin 
lip and have not previously been considered. 

The theory indicates that the energy difference 
uarked E can vary with orientation from zero 
Onearly the Landau spacing, depending upon 

te mass parameters involved. Our data have 
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FIG. 3. Energy level diagram showing two Landau 
levels and the effect of spin splitting for a group of 
conduction electrons in bismuth and antimony. 
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from cyclotron resonance in bismuth? predict 
resonances of type 1 at 2.1 and 1.7 kilogauss, 
respectively. The peak observed at 2.0 kilogauss 
is interpreted as a type 2, or spin, resonance, 
and the close agreement between these fields 
suggests that E is very small in this orientation. 
The fact that the experimental peak appears 
doubled is thought to be due to slight sample 
misorientation. There are three other peaks, 
one of which is tentatively believed to arise 

from hole spin resonance and the other two from 
type-4 transitions, one for electrons and the 
other for holes. The exact positions of the reso- 
nances are in some doubt, both because of their 
close spacing and because the resonance field is 
shifted from the peak, according to the theory 

of conduction electron spin resonance.’»® In 
spite of this uncertainty, however, the data lend 
considerable credence to the above interpretation. 
Transitions of type 3 as well as of types 2 and 

4 have been observed in bismuth at other orien- 
tations. The observation of type-3 transitions 
implies that E+0 is bismuth, at least in some 
orientations, and the value of this separation 
determined from our data as taken with the mag- 
netic field along a dyad axis agrees with a split- 
ting seen earlier in de Haas-van Alphen type 
oscillations in the electronic specific heat by 
Kunzler and Boyle.” The g factor observed with 
the field along this axis is about -200 for the 
low-mass electrons. The resonance peak at 2.1 
kilogauss with the magnetic field along a three- 
fold axis suggests a g factor of about -30 for 
electrons. 

The antimony data show cyclotron resonance 
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at 2.28 kilogauss, a field which is obtained by 
extrapolating the sharply sloping portion of the 
curve to zero absorption, and which corresponds 
to an electron cyclotron mass of 0.088m,. This 
method for identifying the cyclotron resonance 

is justified in reference 5. The sharp resonance 
at 1.14 kg is identified as a type-4 transition and 
the fact that it is at exactly half the cyclotron 
resonance field indicates that E is zero within 
the accuracy of the experiment. A dip preceding 
the sharp rise at cyclotron resonance and a kink 
in the slope of the sharp rise indicates that the 
pure spin resonance line is present at exactly 
the cyclotron frequency, as is expected if E =0. 


————__ 


The g factor for this resonance is 23. 
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‘mM. H. Cohen and E. I. Blount, Phil. Mag. (to be 
published) . 

5J. K. Galt, W. A. Yager, F. R. Merritt, B. B. 
Cetlin, and A. D. Brailsford, Phys. Rev. 114, 1396 
(1959). 

SF. J. Dyson, Phys. Rev. 98, 349 (1955). 

'w. S. Boyle, F. S. L. Hsu, and J. E. Kunzler, 
following Letter [Phys. Rev. Letters 4, (1960)). 





SPIN SPLITTING OF THE LANDAU LEVELS IN BISMUTH 
OBSERVED BY MAGNETOTHERMAL EXPERIMENTS 


W.S. Boyle, F. S. L. Hsu, and J. E. Kunzler 
Bell Telephone Laboratories, Murry Hill, New Jersey 
(Received February 17, 1960) 


In two previous oral presentations’ we have 
described a new type of manifestation of the 
oscillatory dependence on magnetic field of the 
density of states of an electron gas at the Fermi 
surface. Specifically, there is a periodic (in 
1/H) dependence of the electronic heat capacity 
on magnetic field as well as an orientation de- 
pendence at constant field if the Fermi surface 
is nonspherical. The previous reports described 
some observations that have been made on a 
single crystal of bismuth and it was pointed out 
that this type of experiment offers certain ad- 
vantages in sensitivity, simplicity, and inter - 
pretation over other types of magneto-oscillatory 
experiments. In this Letter new data on bismuth 
at larger magnetic fields are presented which 
we interpret to give information on the spin split - 
ting of the Landau levels. 

Only an essential outline of the experiment 
necessary for this report will be given here. 
However, a more complete description and the 
detailed experimental results are in the course 
of preparation. The oscillatory component of 
the total heat capacity is small (usually less 
than 0.1% in bismuth) at 1.4°K. Therefore, 
direct heat capacity observations are imprac- 
tical and the approach adopted is to observe the 
changes in heat capacity as a function of the 
magnitude or orientation of the magnetic field. 
Experimentally these appear as temperature 
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changes under essentially adiabatic conditions. 
This is accomplished by mounting the oriented 
bismuth crystal in a suitable powder’ selected 
to provide a sample to helium bath thermal re- 
laxation time of several minutes. The tempera- 
ture differences are continuously recorded by 
using a carbon thermometer in good thermal 
contact® with the bismuth crystal. 

These cooling peaks occur when the bottom of 
a Landau sub-band passes through the Fermi 
surface and result from the increased density 
of states at the bottom of each band and hence 
an increased entropy of the electron gas. With 
a nonspherical Fermi surface this condition de- 
pends on both magnitude and direction of the 
magnetic field as given by the following equation: 


E =(n,+4+4 A)heH /m*c. (1) 


ie | 
4 is a parameter that determines the spin split- 
ting and is given by 


A=m"*g/4m, (2) 


where m* and gare, respectively, the cyclotron 
mass and the spectroscopic splitting factor for 
the pertinent orientation and m is the free elec- 
tron mass. my is the quantum number of the 
Landau level, Ey is the Fermi energy, and all 
other symbols have their usual significance. 
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In cases where m*<« m, using a free spin g 
factor of 2 leads to a value of A=0. Examination 
of Eq. (1) then shows that a plot of nf against 
1/H should give an intercept near a half-integral 
_ value of nf. However, the de Haas—van Alphen 
‘ev. | susceptibility data of Dhillon and Shoenberg* 
show that the intercept is considerably different 
than 1/2. Our magnetothermal data’ show that 
the intercept is within a few percent of being an 
integer when the magnetic field is in a low-mass 
6 direction. This is consistent with Cohen and 
Blount® who have shown theoretically that for the 
low-mass directions, the g factor in bismuth is 


g=2m/m"*, (3) 


which makes A~1/2. The effect of the various 
factors discussed above on the separation of the 
Landau levels is illustrated schematically in 
Fig. 1. 

When we consider data at high enough magnetic 
fields so as to include one of the  =0 levels, the 
splitting of the sub-bands can be seen as shown 
in Fig. 2. The value of A can be derived from 
the pair of peaks near 15 kilogauss. Since the 
separation of the peaks varies as H”, the lower 
i field peaks are not resolved. The resolution is 
| somewhat better on rotational diagrams. Figure 
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FIG. 1. Effect of various situations on the spin split- 
, ting of Landau levels for a given value of crystal mo- 
ro mentum along the magnetic field. (1) Spin is neglected 
r ad therefore A=0. (2) g=2 is used for electrons of 
co mass much smaller than the free electron mass. (3) 





The g factor for electrons in bismuth is approximately 
equal to but less than 2m/m*. (4) g=2m/m* and 
4=1/2, 


3 is a rotational diagram covering 180° taken at 
15.5 kilogauss with the magnetic field perpendicu- 
lar to the threefold axis. The pronounced peak 
with the field along the binary axis at 0° (and at 

+ 60°) is the 15.5-kilogauss peak of Fig. 2 and 
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MAGNETIC FIELD-—> 


FIG. 2. Cooling peaks as a function of magnetic 
field for the field along the binary axis. The values 
of magnetic field noted for each peak are from rota- 
tional diagrams which allow more precise assignments. 
The spin splitting is evident for the two peaks near 
15 kilogauss but not for peaks at lower fields. 
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FIG. 3. Rotational diagram about the threefold 
axis at 15.5 kilogauss showing angular dependence of 
cooling peaks. 
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represents contributions from four equivalent 
groups of electrons. Each of the major peaks at 
+ 8° originates from two of the four groups of 
electrons that are responsible for the 14.2-kilo- 
gauss peak of Fig. 2. The other peaks in Fig. 3 
have been unambiguously accounted for but are 
not pertinent to the present report. A series of 
rotational diagrams at fields near the 7.2-kilo- 
gauss peak of Fig. 2 provides just enough 
resolution to demonstrate the splitting of that 
pertinent level. 

From the positions of the resolved peaks of 
Fig. 2, 14.2 and 15.5 kilogauss, and assuming 
that A is less than 1/2, we obtain 


A =0.46, (4) 


for the low-mass electrons and the field along 
the binary axis. 

A value of A different from 1/2 suggests that 
an effect should be observed in cyclotron reso- 
nance and spin resonance observations. Smith, 
Galt, and Merritt® have observed transitions 
which they interpret in a manner consistent with 
these ideas. 





When the 14.2- and 15.5-kilogauss peaks are 
included on the plot of 1/H vs peak number of the 
unresolved peaks at the lower fields, one would 
expect that one peak would fall about 4% above 
and the other 4% below the curve. However, it 
is found that the actual values for the magnetic 
fields are about 4% higher for the peaks. This 
is well outside the experimental uncertainty of 
the positions of the peaks at lower field (+ 1%) 
and is ascribed to oscillations in the Fermi 
energy which become large for the last level. 





1J. E. Kunzler and W. S. Boyle, Bull. Am. Phys. 
Soc. 4, 168 (1959); J. E. Kunzler, Bull. Am. Phys. 
Soc. 4, 319 (1959). 

23. E. Kunzler (to be published). 

3w. S. Boyle and K. F. Rodgers, J. Opt. Soc. Am. 
49, 66 (1959). 

J.S. Dhillon and D. Shoenberg, Phil. Trans. Am. 
248, 1 (1955). 
~ 5M. Cohen and E. I. Blount, Phil. Mag. (to be pub- 
lished). 

8G. E. Smith, J. K. Galt, and F. R. Merritt, pre- 
ceding Letter [Phys. Rev. Letters 4, (1960)]. 





F-BAND STRUCTURE IN CESIUM BROMIDE 
Herbert Rabin and James H. Schulman 





Naval Research Laboratory, Washington, D. C. 


The F center in NaCl-type alkali halides has 
been extensively studied,’ but comparatively 


little is known of the F center in CsCl-type alkali 


halides. It has been reported recently that the 
principal optical absorption bands in two CsCl- 
type alkali halides, CsCl and CsBr, show com- 
ponent band structure at liquid helium tempera- 
ture.” It is the purpose of this Letter to report 
current results obtained on CsBr which indicate 
that this structure is due to a single center, 
most likely the F center. 

Figure 1 shows the optical absorption spectra 
of two CsBr single crystals, one a Harshaw 
crystal, additively colored with potassium,® and 
the other an MIT-grown crystal,‘ x-ray irradia- 
ted for 90 minutes (43 kvp and 20 ma) at room 
temperature. The spectra over the range 350 to 
750 mu are shown at room temperature, curves 
(a) and (c), as well as at liquid helium tempera- 
ture, curves (b) and (d). The structure in the 
principal absorption band of curve (b) is shown 
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in detail in Fig. 2 over the range 550 to 700 mu. 
It is clear that the absorption spectra in both the 
additively colored crystal and the x-rayed crys- 
tal are quite similar. At room temperature the 
principal absorptions occur at the same peak 
position, 680 my, and have nearly the same 
half-width of about 0.36 ev. At liquid helium 
temperature corresponding structure in both 
crystals appears at approximately 387, 490, 
580, 620, and 640 mu. Figure 2 shows that the 
main components in the additively colored crys- 
tal can be resolved® into bands which peak at 
616 and 642 mu. A detailed comparison of the 
x-rayed and additively colored crystals indicates 
that this major structure is present in substan- 
tially the same relative proportion in both crys- 
tals. No particular significance is attached to 
the fact that the 642-my band in the additively 
colored crystal (Fig. 2) is slightly displaced 
from the 640-my absorption in the x-rayed crys 
tal (Fig. 1); a small underlying absorption could 
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FIG. 1. The optical absorption spectra of Harshaw 
CsBr, additively colored with potassium, at room 
temperature (a) and liquid helium temperature (b), 
and MIT-grown CsBr x-rayed for 90 minutes at room 
temperature, measured optically at room temperature 
(c) and liquid helium temperature (d). 


easily account for this difference. For a similar 
reason, the relative proportion of the weak ab- 
sorptions at shorter wavelengths in Fig. 1 has 

not been analyzed in detail, although it appears 
that these bands are present in roughly the same 
relative proportion to the major components. It 
has also been observed that the 580-, 620-, and 
640-my bands occur in both Harshaw and MIT- 
grown CsBr that have been x-rayed directly at 
liquid helium temperature, and again they appear 
inthe same relative proportions as the data of 
Figs. 1 and 2, even after successive optical 
bleachings. The presence of short-wavelength 
structure corresponding to the 387- and 490-my 
bands has not been verified after liquid helium 
temperature x-raying due to the introduction of 
other overlapping bands by the low-temperature 
irradiation. 

From what is known of the behavior of the F 
tenter in NaCl-type alkali halides the identifica- 
tion of the room-temperature 680-my absorption 
tand in CsBr as the F band is almost certainly 
torrect.*»? It is the principal optical absorption; 
appears to be a simple band; it bleaches under 
visible light; it is identical in both the x-rayed 
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FIG. 2. The major absorption structure of additively 
colored Harshaw CsBr at liquid helium temperature 
(solid line) and its resolved substructure (dashed line). 


and additively colored specimens; and its peak 
position obeys a Mollwo relation® with what ap- 
pears to be the F band in cesium halides of the 
CsCl structure. Furthermore, as the tempera- 
ture is lowered a progressive shift to shorter 
wavelengths and a sharpening is observed.* The 
fact that the band does not remain simple but 
shows substructure at liquid helium temperature 
is contrary to the normal behavior of the F band 
in NaCl-like alkali halides.* The evidence, how- 
ever, for believing that this structure (or at 
least a portion of it) is related to the F center is 
the following: First, the main components are 
always present in definite proportion independent 
of the manner in which the absorption is pro- 
duced. This would be unlikely if only one of the 
components was the F band and all the others 
were unrelated to the F band. Second, the same 
structure appears in both x-rayed and additively 
colored specimens. This situation is again most 
easily understood if a single excess-electron 
center like the F center is involved.*° The pos- 
sibility of associating the band structure with 
other centers, like the M and R centers, is 
contradictory to the known relationship of these 
centers to the F center in the NaCl-type systems. 
In general they appear as less intense bands on 
the long-wavelength side of the F band, they are 
not necessarily in fixed proportion relative to 
one another, and they are not formed on irra- 
diating at liquid helium temperature. 

The possibility that these sub-bands represent 
different excited states of the F center imme- 
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diately suggests itself. The general decreasing 
strength of the first three major components in 
moving to shorter wavelengths is in accord with 
this suggestion,’ despite the fact that the oscil- 
lator strength” of the main band is considerably 
lower than values found for the oscillator strength 
of the F center in NaCl-type alkali halides. Ad- 
ditional measurements, particularly measure- 
ments of spin resonance and the temperature de- 
pendence of photoconductivity, could check on 

the possibility of multiple excited F-center states 
and give further insight into the nature of these 
states. 





‘See F. Seitz, Revs. Modern Phys. 26, 7 (1954). 

2H. Rabin, J. H. Schulman, and R. F. Marzke, 
Bull. Am. Phys. Soc. 5, 48 (1960). 

8This crystal was kindly provided by Dr. W. G. 
Maisch of the University of Maryland. 

‘Kindly provided by Professor A. Smakula and Dr. P. 
Avakian at the Massachusetts Institute of Technology 
who are also currently investigating color center for- 
mation in the cesium halides. 

‘This resolution is performed by assuming the chief 
band peaks at 642 my, shorter wavelength bands do not 
extend to wavelengths longer than 642 my, and the 
band is symmetrical on an energy scale. Similarly, 
it is assumed the second band which peaks at 616 my 
is not overlapped to the long-wavelength side of 616 my 
by shorter wavelength structure, namely the 580-myp 
band, and it is also symmetrical. It is clear from 
Fig. 2 that the structure remaining does not give a 
simple symmetrical band at 580 mu. However, it is 
to be remembered that the background absorption and 


———_ 


errors in separating the components are lumped in this 
short-wavelength remainder. 

Sp. Avakian and A. Smakula, Bull. Am. Phys. Soc, 
5, 48 (1960). 

‘Cc. J. Rauch and C. V. Heer, Phys. Rev. 105, 914 
(1957). 

SAt liquid nitrogen temperature there appears to be 
an increased distortion in the general structure of the 
band, but no clearly resolved substructure appears as 
it does at liquid helium temperature. 

*See, for example, W. H. Duerig and J. J. Markham, 
Phys. Rev. 88, 1043 (1952). 

The possible role of impurity atoms can not be 
totally ignored. H. W. Etzel and J. H. Schulman, J, 
Chem. Phys. 22, 1549 (1954), have shown that the 
same optical impurity bands produced in x-rayed NaC] 
doped with silver can be produced by electrolytic 
coloration. It is well known, of course, that the addi- 
tive coloration process is quite analogous to that of 
electrolytic coloration. Since both Harshaw and MIT- 
grown CsBr crystals have been used in the present 
study, their impurity content would not necessarily be 
expected to be identical. 

‘lf. Liity has recently discussed possible excited F- 




























band states in NaCl-type alkali halides in which the 
relative magnitudes of the bands are not of successively 
decreasing absorption strengths in passing to shorter 
wavelengths. This work was presented at the Interna- 
tional Symposium on Color Centers in Alkali Halides, 
Oregon State College, September, 1959 (unpublished). 

'2Rauch and Heer’ give the F-center oscillator 
strength as 0.26 assuming a single Gaussian band at 
room temperature. Since the 642-my band constitutes 
about 69 % of the total area of the absorption structure 
of Fig, 2, the oscillator strength of the first excited 
state would be 0.18; similarly the second excited state 
(616-my band) would be about 0.05. 
















In recent years there has been considerable 
interest in group III-group V compound semi- 
conductors as well as large band gap semicon- 
ductors. Cubic boron phosphide is such a ma- 
terial. It has been prepared and identified pre- 
viously’ and the present paper describes prelim- 
inary measurements of its semiconducting pro- 
perties. 

Cubic BP has been prepared by two basically 
different methods: (1) by the vapor phase reac- 
tion of boron compounds with phosphorus com- 
pounds or phosphorus vapors; and (2) by the 
slow cooling of a melt containing dissolved boron 
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SEMICONDUCTING PROPERTIES OF CUBIC BORON PHOSPHIDE 


B. Stone and D. Hill 
Research and Engineering Division, Monsanto Chemical Company, Dayton, Ohio 
(Received February 15, 1960) 



















and phosphorus. 

The BP prepared by the vapor phase reactions’ 
has very small crystal size and varies from 
reddish tan to black in color. However, this 
material could be deposited in the form of films 
which were useful for optical absorption meas- 
urements since they could be made very thin. A 
variation of the vapor phase reaction resulted in 
small “needle” or “blade” shaped crystals. 
These crystals are 1 mm or 2 mm in length and 
reddish brown in color and transparent. 

The method involving slow cooling of a melt 
results in crystals which are on the order of 
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1 mm in all dimensions and of a similar color. 
These crystals are also transparent and occa- 
sionally a crystal can be found with a pair of 
flat, parallel faces. Measurements of the dis- 
placement of an image by such a crystal give a 
rough value of index of refraction of between 3.0 
and 3.5. 

The x-ray powder diffraction pattern of these 
materials is identical to that reported for cubic 
BP by other workers.* Elemental analyses give 
boron to phosphorus ratios of unity within about 
1% which is the accuracy of the method. 

BP is quite inert chemically. It is not attacked 
by mineral acids or any of the usual semicon- 
ductor etches. It is attacked slightly by molten 
NaOH. BP has a hardness about equal to that of 
silicon carbide (Knoop hardness number of 3200 
for 100-g load). 

Ruehrwein and Williams’ of this laboratory 
have found by means of optical absorption meas- 
urements on thin polycrystalline films, that the 
energy gap of BP is almost 6 ev. Figure 1(a) 
shows representative measurements’ on films of 
different thicknesses. If one neglects scattering 
effects, the absorption coefficient at 0.3 micron 
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FIG. 1. (a) and (b) Percent transmission of poly- 
ttystalline films of boron phosphide. 












is about 5x10‘ cm™*. These films were deposited 
on quartz plates. The absorption of the quartz is 
negligible in this region. The value of energy 
gap of about 6 ev is consistent with what one 
would expect on the basis of lattice parameter 
considerations. Note the shoulder around 0.6 
micron accounting for the reddish color which 
this material exhibits by transmission. 

Absorption measurements on a polycrystalline 
film in longer wavelength regions are shown in 
Fig. 1(b). For this measurement the film was 
stripped off its quartz backing. The gradual in- 
crease in transmission out to about 10 microns 
is probably due to a decrease in scattering at 
longer wavelengths. The sharp absorption at 
12.1 microns is presumably a lattice absorption. 

Optical absorption measurements of single- 
crystal samples are more difficult because of the 
small size of the crystals available and the dif- 
ficulty in producing thin samples. The sensitivity 
is too poor to show the fundamental absorption 
and the only feature observed is the absorption 
responsible for the color of the crystals. The 
highest absorption coefficient measured on a 
single-crystal sample is 10° cm“?. 

Preliminary Hall measurements have been 
made on crystals grown by method 2. The sam- 
ple geometry was determined by their natural 
crystal faces. Electrical contact was made 
using silver paint, tungsten, or platinum points. 
The lead resistance could usually be decreased 
considerably by the use of a capacitor discharge. 
The measurements show that these crystals are 
p-type with a carrier concentration ranging from 
1 to 5x10** carriers/cc. The Hall coefficient is 
essentially constant from about 900°K down to 
about 160°K. Difficulties with leads prevented 
Hall measurements at lower temperatures. The 
shape of the resistivity curve at low tempera- 
tures suggests that the material is still in the 
exhaustion range at 78°K. Figure 2 gives Hall 
mobility on such p-type BP as a function of 
1000/T. At room temperature the mobility is 
decreasing with decreasing temperature, char- 
acteristic of impurity scattering. Thus the room 
temperature lattice mobility of p-type BP must 
be considerably greater than the measured value, 
perhaps as large as 300 or 500 cm?/volt second. 

Thermoelectric power and rectification meas- 
urements on the needle-shaped crystals indicate 
that these crystals are n-type. Crude resistivity 
measurements indicate that the impurity level 
may be around 10” carriers/cc if we assume 
the electron mobility to be of the same order of 
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e with both n-type and p-type BP. Rectification 
ratios of greater than 10° to 1 have been observed 
at both room temperature and 400°C. 

It is necessary to comment concerning the 

OO} color of the BP obtained to date. A semicon- 

os 7 ductor with an energy gap of nearly 6 ev should 
2 7 be transparent and colorless, while the BP ob- 
= r / tained has a reddish color. The optical absorp- 
S 7 / tion responsible for this color occurs at a - 
= / po wave 
“ F length corresponding to an energy of about 2 ev, 
§ { Y Such a level would not be ionized at the tempera- 
2 ture of the measurements, so it would not con- 
b / tribute carriers, even if this were possible. 
a r There are several possible explanations for this 
absorption but more work will be necessary be- 
‘ fore any conclusions can be reached. 
10 4 In summary it has been shown that it is possi- 
¥ \ ble to produce single crystals of cubic BP of 

- measurable size, and although the material pre- 
7 pared so far is impure, by semiconductor stand- 

{ ards, and although the measurements just pre- 

L 1 l Sea es 2 
| 9 3456810 5 sented are preliminary in nature, the properties 
of BP appear to be interesting and promising. 
l000/T (%y! 

FIG. 2. Hall mobility of single-crystal p-type boron 

phosphide. Ip, Popper and T. A. Ingles, Nature 179, 1075 
"' 3 Willi ted at A i 
. . llllams, paper presen a merican 
magnitude as the hole mobility. Chemical Society Sicten, Seaten, Massachusetts, 

The thermoelectric power of both the ”- and April, 1959 (unpublished). 
p-type material is about 300 microvolts/°C. 3R. A. Ruehrwein and F. V. Williams (unpublished 
Point contact rectification has been observed data). 

MECHANISMS OF VOLUME SELF-DIFFUSION IN a-Fe AND y-Fe! 
C. J. Meechan 
Atomics International, Canoga Park, California 
(Received January 11, 1960) 

Several years ago, Kuczynski’ introduced a for each point on the InDvs 1/T curve, the sta- | 
method for determining self-diffusion coefficients tistics will be favorable to enable the determina- , 
without the aid of radioactive tracers. An ex- tion of D within a factor of 2 or 3 with about an 
pression was derived which relates the dimen - 80% confidence limit for most metals. 
sions of a sintered interface between two parti- An adaptation of this method has been used to 
cles to D, the coefficient of volume self-diffusion. determine the self-diffusion parameters in both ' 
Although this method will never supplant the face-centered cubic (fcc) and body-centered 
tracer method for determining diffusion para- cubic (bec) thorium,” where no suitable isotope f 
meters, the importance of its use should not be for tracer studies is available. In this latter ‘ 
overlooked in diffusion studies (a) where no suit- work, it was found that near the transition tem- t 
able isotopes are available and (b) where it is _ perature, D in bcc thorium is approximately 4 I 
desired to distinguish various diffusion phenom- factor of five less than D in fcc thorium. In 
ena which have different causes. As long as 50 contrast, radioactive tracer studies*»* show that 
to 100 measurements of separate fillets are made D in bec iron is approximately 400 times larger W 
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than D in fcc iron, near the transition tempera- 
ture. It was therefore suggested’ that sintering 
studies be conducted on iron in an attempt to 
determine the cause of this extremely large dif- 
ference in D for the fcc and bcc structures. 

It is important to note that the volume diffusion 
mechanism responsible for sintering phenomena 
must involve lattice defects, e.g., interstitials 
or vacancies, since a net transport of matter 
must take place to establish the sintered fillet. 
On the other hand, the total diffusion which is 
observed through use of radioactive tracers may 
have occurred either by such a defect mechanism 
or by a direct interchange or ring mechanism, ° 
or possibly by a combination of both kinds of 
mechanisms. The direct interchange mechanism 
cannot contribute to the formation or growth of 
the sintered fillet since it does not result in any 
net mass transport. Consequently, a comparison 
of diffusion data obtained from tracer and sinter- 
ing studies should allow one to determine the 
relative contributions of those diffusion mech- 
anisms which can produce net mass transport 
and those mechanisms which cannot. 

In the present study, 0.020-inch diameter 
wires were fabricated from Johnson- Matthey 
(99.99% pure) iron. These wires were wrapped 
in four spools, each containing about four layers, 
in such a way that approximately one hundred 
possible sintered interfaces were available for 
examination for each cross-sectional cut. These 
spools were pre-sintered in hydrogen at 950°C 
for one hour. Two of the spools were sectioned 
and measurements were made on the sintered 
interfaces. These results established the di- 
mensions of the fillets at “zero” time; this pro- 
cedure also eliminates errors due to distorting 
the interface region during wrapping of the spools. 
One of the sectioned spools and one of the whole 
spools were then heated in vacuum (< 10°°mm Hg) 
for 100 hours at 935°C; the other spools were 
treated similarly at 885°C. Following these 
treatments, final measurements of the sintered 
interfaces were made and the growth of each 
fillet in this time period was determined by sub- 
tracting the initial from the final values. Ap- 
proximately 100 separate fillets were measured 
for each temperature and these values were av- 
eraged to obtain the value from which D was deter- 
mined. The previously derived formula? relating 
D to the fillet measurements is 


D=ka‘T(2655ty)~*y*f(y), 


where k = Boltzmann constant, a=wire radius, 


5=interatomic distance, T=absolute temperature, 
t=time, y=surface energy, y = ratio of fillet 
width and wire diameter, and f(y) =0.0637 
+ 0.1096y+ 0.1426y? + 0.1668y5+... 

For y = 2000 dynes/cm, the values of D deter- 
mined from this expression for the two temper- 
atures are: 


Dogs = 3-5 10°-** cm?/sec (fcc), 
Dogs < 6X 10°** cm? /sec (bec). 


Although some sintering did occur at 885°C, the 
diffusion rate was so slow that after the 100-hour 
treatment, it was only possible to place an upper 
limit on D,,,. These values are plotted in Fig. 1 
along with the previously mentioned data obtained 
using radioactive tracers. More recent tracer 
data® on a-iron are also shown. Approximately 
85% of the present individual measurements of 
D4, fall within the error bar shown in Fig. 1. 
The sintered fillet is a result of contributions 
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FIG. 1. Plot showing comparison of InDvs 1/T 


for a-Fe and y-Fe as determined from radioactive 
tracer data and the present sintering data (,7). 
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from several diffusion mechanisms which employ by a change in the dominant diffusion mechanism, 


a defect mechanism of mass transport. If vol- It seems most likely that this change is from a 
ume diffusion is not predominant in establishing four-atom ring mechanism in bcc iron to a va- 
the fillet, then the value of D,,, corresponding cancy mechanism in fcc iron. This conclusion 
to volume diffusion must be even smaller than is further substantiated by recent diffusion studi 
that given in the paper, thus increasing the dis- on a-iron in a high-temperature gradient, ® and 
crepancy between D,,, determined from sintering is direct experimental support for the theoretica| 
data and the corresponding value determined treatment of Le Claire® on diffusion mechanisms 
from tracer data. This strengthens the main in bcc metals. 


argument, namely, that different mechanisms of 
diffusion are dominant in a-Fe and y-Fe. 





From the preceding discussion and the fact This work was performed under contract with the 
that D,,, from sintering data is in good agree- U. S. Atomic Energy Commission. 
ment with the tracer data, it is concluded that 1G. C. Kuczynski, J. Appl. Phys. 21, 632 (1950). 
(a) the sintered fillets in the fcc structure were *C. J. Meechan, Advances in Nuclear Engineering, 





edited by J. R. Dunning and B. R. Prentice (Pergamm 
Press, New York, 1957), p. 209. 
3c. E. Birchenall and R. F. Mehl, Trans. Am, 


established by volume self-diffusion via a defect 
mechanism, presumably vacancies, and (b) a 


diffusion mechanism which does not involve lat- Inst. Mining Met. Petrol. Engrs. 188, 144 (1950). 
tice defects is primarily responsible for volume ‘F. S$. Buffington, I. D. Bakalar, and M. Cohen, 
diffusion in bcc iron. This latter mechanism is Trans. Am. Inst. Mining Met. Petrol. Engrs. 188, 
presumed to be a ring mechanism such as pro- 1375 (1950). ia 
posed by Zener, where two or more neighboring °C. Zener, Acta Cryst. 3, 346 (1950). 


6 
atoms exchange positions simultaneously. The : cod i. Birchenall and R. J. Borg (private commu- 
nication). 


contribution from a four-atom ring seems most 1 ’ 
favorable from theoretical energy considerations.’ C. Zener, Imperfections in Nearly Perfect Crystals, 
by . edited by W. Shockley (John Wiley & Sons, Inc., New 








Therefore, it is concluded that the large change York, 1952), p. 289. i 
in D, which occurs when iron is heated through ’w. G. Brammer, Acta Met. (to be published). 
the alpha-gamma critical temperature, is caused °A, D. Le Claire, Acta Met. 1, 438 (1953). 





SHELL EFFECTS IN THE OPTICAL POTENTIAL 


Atsushi Sugie 
Nuclear Physics Department, Japan Atomic Energy Research Institute, Tokai, Ibaragi-ken, Japan 
(Received December 7, 1959) 


There are a few qualitative disagreements reason because for A ~40, where both proton 
between the observed low-energy neutron scat- and neutron are magic, the experiments agree 
tering data (the strength function (T'°/D) and the fairly well with the theory (Fig. 1).* Phenomena 
scattering length R’) and the predictions from (2) and (3) were attributed to the quadrupole 
the conventional optical potential whose imaginary deformation of the nucleus by Seth,* as was done 
part (W) is independent of the mass number A. for A~150. This again is doubtful since, for 
(1) The observed (I°/D) is too small around A~60, the quadrupole moments are not too large 
A~100,' that is, W should be smaller than usual to be explained by configuration mixing and the 
in this region. (2) The observed R’ is too small low -lying levels are not of rotational type. 
around A ~60,? (Fig. 1), which means that W We propose to explain the above phenomena 
should be larger than usual in this region. (3) by the model of Lane and Wandel.*® Usually this 


The observed (I°/D) is too small around A ~ 60, ” model is applied to an infinite medium. Now we 
(Fig. 1), and this also indicates that W is larger apply it to a finite nucleus. The theoretical 





than usual here. justification for this was given before.® What is 
Lane et al.’ attributed the phenomenon (1) to the difference between the two? The difference 
the fact that for A ~100 either proton or neutron appears in the conservation law for the scatter- 


is near magic. This will not be the principal ing process leading to the absorption. For an 
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FIG. 1. The observed strength function (I/D) and 


the scattering length R’.? The solid curve is the pre- 
diction from the conventional optical potential [V. F. 
Weisskopf, Physica 18, 952 (1956)]. The dashed curve 
is the one from the black nucleus. 


infinite nucleus, we calculate the effective two- 
body scattering cross section in the nuclear 
matter for which the momentum and energy are 
conserved. For a finite system, on the other 
hand, it is the angular momentum, parity, and 
energy which are to be conserved. The parity 
conservation is very important since, according 
to the shell model, the single-particle levels 
with the same parity are grouped together’ (ex- 
cept for one level with the highest angular mo- 
mentum) and the conservation law is either yes 
orno. The conservation of the angular momen- 
tum is not so restrictive because of its vectorial 
nature. 

Figure 2 represents the schematic single- 
particle levels.” In the figure the filled levels 
are indicated by curly brackets and the positions 
of the zero-energy s-wave neutron are indicated 
by dashed lines for A ~60 and 100. For A ~60: 
(i) The incident s neutron can easily jump down 
to the 3s state just below and raise a nucleon in 
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FIG. 2. The schematic single-particle levels.’ 
The filled levels are indicated by curly brackets, and 
the positions of the zero-energy incident neutron are 
indicated by dashed lines for A ~60 andA~100. For 
the latter, the top of the filled levels is different for 
protons and neutrons. 


the target to a nearby state. (ii) If the incident 
neutron jumps down far below and the parity 
change is yes, a nucleon in the target is raised 
high so that the energy is conserved and the 
parity change is again yes. So we expect large 
W. For A~100: (i) When the incident s neutron 
jumps down to the nearby states, the parity 
change is yes, but for the raised nucleon in the 
target the parity change is essentially no. (ii) 
If the incident neutron jumps down far below, 
the nucleon in the target is raised high and the 
parity law is again not satisfied. Thus we ex- 
pect small W. In this way we can understand (1), 
(2), and (3) as the shell effect, but the meaning 
is quite different from that of Lane et al. The 
sharp rise of (T°/D) and R’ near A = 70 (Fig. 1) 
may be due to the filling of the minus parity 
levels.® 

According to the present theory, W for the p 
state must be small for A ~ 60 and large for 
A~100. There are no reliable data for W in the 
p state. Seth,* and also Lane et al.,’ gave some 
experimental data indicating that W is small for 
A~100 in contradiction to our theory, but these 
are not convincing. An exact analysis for W in 
the p state is highly desirable. 

The surface absorption model® also makes 
(T°/D) small for A~100 but, according to our 
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preliminary machine calculations, the reduction 


is not enough. Further, this model does not 
explain (2). The nonlocal nature® of the imagin- 


ary part may also be responsible for making the 
effective local W dependent on A and the incident 


angular momentum. This point has not been 
studied yet. 





1A, M. Lane, J. E. Lynn, E. Melkonian, and E. R. 
Rae, Phys. Rev. Letters 2, 424 (1959). 

"kK. K. Seth, Revs. Modern Phys. 39, 442 (1958). 
Figure 1 is reproduced from a private communication 
from Seth. 

34. M. Lane et al. stated that W is small since the 
level density is small near magic. W depends, how- 
ever, also on the nuclear matrix elements which will 
be large if the level density is small. Further, we do 
not observe large difference between (I/D) for the 
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even nucleus and the odd nucleus although the level 
density differs by about a factor 2. This was pointed 
out to the author by V. F. Weisskopf. 
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near A ~ 50 but, since R’ is normal here, the theoreti- 
cal R’ and (I/D) will be fitted to the experiment by 
slightly changing the depth (or the radius) of the real 
part of the potential, leaving W near the normal value, 

*F. Bjorklund and S. Fernbach, Phys. Rev 109, 
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FURTHER INTERPRETATION OF A MEDIUM ENERGY (, 2p) EXPERIMENT* 


A. J. Kromminga 
University of Minnesota, Minneapolis, Minnesota 


I. E. McCarthy 
University of California, Los Angeles, California 
(Received February 19, 1960) 


In a previous publication,’ the momentum 
transfer distribution observed in the 40-Mev 
(p,2p) angular correlation experiment of Grif- 
fiths and Eisberg’ was interpreted as giving in- 
formation about the momentum distribution of 
nucleons localized in the nuclear surface. After 
unfolding the momentum transfer distribution 
due to the optical potentials in the entrance and 
exit channels, it was found that momentum com- 
ponents corresponding to proton energies of 
about 1 Mev were most likely to be observed for 
the bound particle. It is the purpose of this 


Letter to point out some additional considerations 


which contribute to an understanding of the re- 
action. 

A collision of a nucleon with a stationary free 
nucleon will result in an angle of 90° between the 
momenta of the particles in the final state. In 
reference 1, it is shown that the effect of the 
binding potential is to move the region of small- 
est momentum transfers to an angle less than 
90°. If small momentum transfers are most 
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probable, the result will be an angular correla- 
tion peaked at an angle less than 90° (depending 
on the binding and excitation energy) as is ob- 
served experimentally. 

In reference 1, it was assumed that the mo- 
mentum distribution was isotropic. The momen- 
tum components were estimated from the experi- 
ment. Essentially the peak in the experimental 
angular correlation is so narrow that only small 
momenta (~1 Mev) contribute. The relationship 
between the width of the angular correlation 
curve and the momentum of the struck particle 
is explained in reference 1. 

It was also shown in reference 1 that the re- 
actions are not isotropically distributed over 
the nuclear surface. They come predominantly 
from an equatorial belt on the surface with the 
incident direction as the north-south axis. Within 
the belt, the main contribution comes from re- 
gions A, Fig. 1, since, for one particle scattered 
to the left and the other to the right, one of the 
particles coming from region B is likely to be 
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FIG. 1. Schematic diagram of a nucleus and ex- 
perimental arrangement showing: P, a plane passing 
through the counters and the incident beam that also 
passes through the center of a nucleus; C, counters; 

I, incident beam; A and B, regions of greater and 
lesser importance in the ($,2)) reaction, respectively. 


absorbed. From simple kinematics with parti- 
cles of equal mass, it may be shown that if the 
incident particle collides with a free particle 
moving at 90° with respect to it in the laboratory 
system, the angular correlation between the 
particles after collision is still given by a 6 func- 
tion at 90° in the laboratory system. In the ap- 
proximation that the collision occurs at the equa- 
tor, only the tangential component of the momen- 
tum of the struck particle in the incident direc- 
tion will contribute to a spreading of the angular 
correlation. From the consideration of proton 
flux in an optical potential® it is observed that 

the incident protons travel roughly tangentially 

in the region of interaction. 

Reactions are most likely to be observed in 
which the struck particle is moving tangentially. 
Radial motion will result in the plane of the col- 
lision either intersecting the nucleus in which 
case at least one particle is likely to be absorbed 
(regions B, Fig. 1) or out of the plane of the 
counters (region A, Fig. 1). 

It has been shown by Baker, McCarthy, and 
Porter* that the shape of the momentum distri- 
bution that may be observed for particles local- 
ized at the nuclear surface is determined almost 
completely by the localization rather than by the 
wave function of the bound particles. The tan- 
gential momentum components for particles in 
higher shell model states are found to be much 
less than the radial components, within the re- 
strictions imposed by the localization. The re- 
sults of this calculation agree with the uncertainty 
principle because the region of interaction is less 
restricted by absorption tangentially than radi- 





ally. Thus a more realistic theoretical estimate 
of the effective momentum distribution than the 
isotropic assumption of reference 1 is possible. 
Using the above momentum transfer components 
in the relationship between momentum transfer . 
and angular correlation (Figs. 2 and 3 of refer- 
ence 1), the experimental angular correlation 
curve may be easily understood both in position 
and width. 

Owing to the imperfection of the tangential ap- 
proximation, there will be a contribution from 
the tangential component perpendicular to the 
incident beam and from the radial component. 
Although they cannot appreciably spread the 
angular correlation, they will cause a prepon- 
derance of events in which the energies of the 
final particles differ if the counters are placed 
fairly symmetrically about the incident direction 
[see Fig. 2]. This effect was observed in the 
experiment. 

The small momentum transfers observed in 
the (p, 2p) experiments are essentially due to the 
anisotropy of the momentum distribution and the 
fact that the interaction does not effectively occur 
uniformly throughout the nuclear surface. Thus 








FIG. 2. Sketch showing conservation of momentum 
under certain conditions: Kp = incident momentum; ky 
=radial component of the momentum of a bound particle 
(region B)= the tangential component of the momentum 
of a bound particle that is perpendicular to the incident 
momentum (region A); k, = final momentum; k;, ky = 
momenta of particles in the final state. Conservation 
of momentum requires k,? >k,”. 
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one would suspect that a theory analogous to the 
Butler theory should not predict the correct, 
rather narrow angular correlation since (a) it 


does not contain information about the momentum 


distribution in different regions of the nuclear 
surface, and (b) it does not distinguish between 
the regions of the nuclear surface as to which 
are more likely to contribute to the reaction. 
Calculations using such a theory confirm this 
prediction. 





‘Supported in part by the National Science Foundation 
and the U. S. Atomic Energy Commission. 

‘I, E. McCarthy, E. V. Jezak, and A. J. Kromming, 
Nuclear Phys. 12, 274 (1959). 

*R. J. Griffiths and R. M. Eisberg, Nuclear Phys. 
12, 225 (1959). 

1, E. McCarthy, Nuclear Phys. 10, 583 (1959). 

‘G. A. Baker, I. E. McCarthy, and C. E. Porter 
(to be published) . 





NUCLEAR MATRIX ELEMENTS IN THE BETA DECAY OF Sb'™T 


R. M. Steffen* 
Department of Physics, Purdue University, Lafayette, Indiana 
(Received February 23, 1960) 


After the recent clarification of the beta inter- 
action it has become of interest to study the re- 
lative contributions of the various matrix ele- 
ments to first forbidden 8 transitions. It is the 
purpose of this paper to demonstrate that an un- 
ambiguous determination of matrix elements in 
a nonunique £ transition (e.g., Sb’**) is possible 
on the basis of precise 8-y directional and B-y 
circular polarization correlation measurements, 
if the 8 transition shows appreciable deviation 
from the ~ approximation.'~* 

The — approximation was first introduced by 
Konopinski and Uhlenbeck’ to explain the statis- 
tical shape of most nonunique first-forbidden 
beta spectra. In this approximation the beta 
transition probability is expanded in powers of 
the nuclear radius R and only the leading terms 
are taken into account, which are associated 
with the Coulomb factor §=aZ/2R. Deviations 
from this approximation may be caused by se- 
lection rule effects which inhibit contributions 
from matrix elements other than { Bij. The con- 
tribution of the f Bi term, which is of rank )\ =2 
and which describes the component of the lepton 
field carrying away two units of angular momen- 
tum, may then become very important. The 
spectra of such £ transitions exhibit deviations 
from the statistical shape and their ft values 
(log ft >10) are considerably larger than the 
characteristic ft values of nonunique first for- 
bidden transitions (log ft = 8). 

It was suggested” * that such a selection rule 
effect rather than a mutual cancellation of 
matrix elements explains the large ft value of 
the 2.31-Mev 8 transition of Sb’ (log ft= 10.6). 
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The results of the present investigation confirm 
this hypothesis. 

The angular and energy dependence of the 8,- 
y, directional correlation of Sb’™ involving the 
8, component of 2.31-Mev maximum energy 
(refer to inset of Fig. 1) was measured with the 
vacuum chamber described previously.* The 
directional correlation W, (6, W,=4.8) of the 
B,-y, cascade measured at a fixed average en- 
ergy of W,=4.8 (in units of mc?) is shown in 
Fig. 1. A least-squares fit of the experimental 
points to the correlation function: 


W (8, W =4.8)=1 +A, (4.8)P,(cosé) 


+A ,(4.8)P,,(cosé), (1) 


yielded the following values for the correlation 
coefficients: 


A,(4.8) = - 0.390+ 0.011, 
A, (4.8) = + 0.004 + 0.013. (2) 


The absence of a P,(cos@) term provides further 
evidence against the decay scheme 4*(8,)2*(y,)0°: 
The dependence of the coefficient A,(W) on the 

8B energy is shown in Fig. 2. A simultaneous 
measurement of the energy dependence of the 
82-2 directional correlation made it possible to 
correct the data for the presence of the £,-7, 
directional correlation at 8 energies below the 
maximum energy of the 8, spectrum (W, = 4.15). 
There is, however, some uncertainty in this 
correction due to the fact that the sign of the £2- 
M1 mixing ratio 6 of the y, transition is not 
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known (6=+ 1.00+ 0.085). The error caused by 
this uncertainty is included in the error flags of 


the experimental points corresponding to W<4.15. pressed by 
All available data indicate strongly that the W?-1 


8,-y, cascade of Sb’* follows the decay scheme A,(W) = 
3"(8,)2*(y,)0*. The four matrix elements which 
can contribute to a first forbidden § transition 
with AJ=+1, yes, are the relativistic matrix 
element y=CyJia, the moment type matrix ele- 
ments x =-CyJr, and u=C, Jioxr (all of rank \ 
=1) and the matrix element z =C, {Bi (A =2). 


FIG. 1. Angular de- 
pendence of the B-y di- 
rectional correlation 
involving the 2.31-Mev 
8 transition of Sb'**, 

The measurements were 
made at an average B 
energy of W=4.8 (in 
units mc?). 


FIG. 2. Energy depend- 
ence of the anisotropy fac- 
tor A,(W) of the B,-y,; di- 
rectional correlation. The 
solid line represents A,(W) 
calculated with the param- 
eters U=-0.01, ¥ =0.08, 
Y=0.38. The dashed line 
corresponds to a pure 
{Bij (unique) transition. 


After Kotani* the energy dependence of the di- 
rectional correlation coefficient A,(W) is ex- 


R,+ew (3) - 





Ww 


1+aW+cW*+(b/W) © 


The coefficients R,,e,a,b, andc, which are 
complicated functions of the matrix element 
parameters x,u,z, and Y =y - (u+x) and of the 
maximum energy W, are given in reference 2. 
By a least-squares method the values of the 
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parameter ratios u/z,x/z, and Y/z were deter- 
mined which resulted in a best fit of the data of 
Fig. 2 to a curve given by Eq. (3). As an addi- 
tional condition it was imposed that the set of 
parameters also satisfy the circular polariza- 
tion correlation measurements of Hartwig and 
Schopper,® i.e., P,.(126°)/P,(160°) = -2, where 
P,(6) is the degree of circular polarization of 
the y, radiation measured at the angle ¢. The fit 
of the data yields for the nuclear parameters: 


u=-(0.01+ 0.04)z, 
x= (0.08+0.08)z, 
Y= (0.38+0.12)z. (4) 


This set of parameters agrees satisfactorily 
with the set obtained by Hartwig and Schopper,® 
which was determined on the basis of a some- 
what different approach. The function A,(W) cal- 
culated with the parameters of Eqs. (4) is re- 
presented as solid line in Fig. 2. For compari- 
son the curve A,(W) corresponding to a unique 
8 transition (pure [{ B;j) is included. 

By taking into account the corrected ft value 
of the 2.31- Mev £ transition, ft=10'°*® sec, or 
ft=3.1x10™ in units h=m=c=1, the absolute 
values of the matrix elements involved in this 
8 decay can be computed’: 


if, | /R = (1.204 0.15) x 1072, 
ifFi/r = (1.24 1.2) x 1073, 

if BxFi/R = (0.14 0.4) x 1078, 
| | ia | = (3.14 2.4) x 1074, 


( fia/[B,,) >. (5) 


The values of the matrix elements are given in 
a form which is independent of the chosen sys- 
tem of units (R = nuclear radius). In addition, 


the lack of overlap of the nuclear wave functions 
which occur in the matrix elements is more evi- 
dent in the form of Eqs. (5). If the wave func- 
tions of the initial and final nuclear states would 
overlap perfectly, the values of J By/R, [r/R, 
and fioxr/R would be of order unity whereas the 
relativistic matrix element fia would be of or- 
der Unucleon/© =9-1. It is interesting to note 
that all matrix elements involved in the Sb 3, 
transition are considerably reduced. Compared 
to unique 8 transitions (AJ=+2), where | / Bi \/R 
is of the order ~0.1, the {Bj; involved in the 
Sb‘ 8, transition is reduced by a factor of about 
10. The reduction of {B;;/R, however, is orders 
of magnitude smaller than the reduction of the 
other matrix elements. The cause of the unusual 
predominance of the { Bij matrix element seems 
to be a selection rule effect as suggested by 
Kotani? and by Morita and Morita.° 

The author is greatly indebted to Professor H. 
Schopper for communicating the B-y circular 
polarization correlation results of reference 6 
prior to publication and for illuminating dis- 
cussions. He also wishes to express his grati- 
tude to K. Alder for many suggestions and sti- 
mulating discussions. 
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'E. J. Konopinski and E. G. Uhlenbeck, Phys. Rev. 
60, 308 (1941). 

*T. Kotani, Phys. Rev. 114, 795 (1959). 

3M. Morita and R. S. Morita, Phys. Rev. 109, 
2048 (1959). 

4R. M. Steffen, Phys. Rev. Letters 3, 277 (1959). 

°T. Lindquist and I. Marklund, Nuclear Phys. 4, 
189 (1957). 

8G. Hartwig and H. Schopper, following Letter 
(Phys. Rev. Letters 4, 293 (1960)]. 

'For this computation the following values of the 
coupling constants (obtained on the basis of the ft 
values of 0-0 transitions) were used: Cy=g=(1.41 
+ 0.01) x 10~* erg cm’, or Cy=g=2.97 x 107" in 
units i =m=c=1, and C4 =-1.25Cy. 
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8 -y CIRCULAR POLARIZATION CORRELATION OF Sb’ 


G. Hartwig and H. Schopper 
Institut fir Kernphysik, Universitat Mainz, Mainz, Germany 
(Received February 23, 1960) 


The first forbidden 3~ -2* transition of Sb'** 
with an end-point energy of 2.3 Mev is of special 
interest as it exhibits a considerable deviation 
from an allowed spectrum and a large £6 -y ani- 
sotropy which cannot be explained by an acci- 
dental cancellation of matrix elements as in the 
decay of RaE. The classical experiments (/t 
value, shape, 8 -y anisotropy) alone do not allow 
one to determine the 4 matrix elements Jr, 
figxt, fid@, and Bij which can contribute to a 
AJ=+1, yes transition. It will be shown that a 
unique determination of these quantities can be 
achieved by measuring the 8 -7y circular polari- 
zation correlation. The set of matrix elements 
derived from the experimental results is in 
agreement with the “modified Bij approximation” 
requiring 5;; matrix elements much larger than 
all the others. The transition of Sb’ investigated 
here is the first uniquely established example of 
the “seiection rule effect” which has been dis- 
cussed by several authors.’»” 

The experimental setup was very similar to 
that used in previous measurements of the cir- 
cular polarization of y rays following allowed 
8 decay. Compton forward scattering from a 
magnetized iron cylinder was used to detect the 
y polarization which is given by P, = w(v/c)cosé, 
where v and c are the velocity of the electrons 
and of light, respectively, @ is the angle between 
electron and y momentum, and w is the sym- 
metry coefficient.* For allowed transitions w 
is only a function of the matrix elements and the 
spins whereas in forbidden transitions it depends 
also on the total electron energy W (measured 
in units of m,c*) and 6. As w changes consider - 
ably if 6 is varied, the instrumental angular 
spread must be kept as low as possible. To 
achieve this, three quarters of the scattering 
magnet was filled with lead, although the coin- 
cidence counting rate is decreased appreciably 
by this procedure. As in addition only a small 
region of the 8 spectrum can be used (a second 
8 spectrum starts at W,=4.1), about two months 
of continuous data collecting were required for 
one measurement in order to get sufficient 
Statistical accuracy. 

The experimental results are presented in 
Table I. Measurements were performed at two 






Table I. Experimental results for (wcos@),,. 











WwW 4.3 3.2 
6 = 126° 0.226 + 0.046 0.190 + 0.060 
6 =160° -0.107 + 0.035 -0. 094 + 0.043 





different angles and for two energy ranges. The 
upper range extending from W=3.8 to W=5.5 
(W =4.3) contains only a very small contribution 
from the second 8 spectrum, whereas in the 
lower energy range (W=3.0 to W=3.8) the co- 
incidence rate due to the low-energy branch 
amounts to about 30%. As the details of this 
second transition are not known, it is difficult 
to infer quantitative results from the measure- 
ments at the lower 8 energy. They show, how- 
ever, qualitatively that w does not depend strongly 
on the 8 energy. On the other hand, w changes 
rapidly with 6. Therefore the most accurate 
determination of the matrix elements will be 
achieved if measurements at different angles 
but at a fixed energy are used rather than by 
taking advantage of the energy dependence. 

The equations relating the observed quantities 
w and € (the 8 -y directional correlation coeffi- 
cent) to the matrix elements are rather compli- 
cated’ but a solution can be found easily in the 
following way. The difference W By(160°)w(160°) 
- Wg,(126°)w(126°) depends only on one combi- 
nation of matrix elements i.e., 2x -u. From the 
experiments at W=4.3, one deduces 2x -u=0+0.19. 
Assuming that 2x =u is rigorously valid, one can 
eliminate x from all equations, leaving the two 
parameters ¢,=Y+(u-x)W,/3 andu. Ina ¢, -u 
plane each experimental value determines a conic 
section (Fig. 1). The intersections of the curves 
for w(@ = 160°) and w(@ = 126°) give two sets of solu- 
tions but only one is consistent with the 6 - y 
anisotropy.° 

Within the experimental errors the three curves 
intersect at one point proving that 2x =u is a good 
approximation.® A better intersection than shown 
in Fig. 1 can be achieved by putting 2x -u=+0.1 
although this change is still within the experi- 
mental uncertainties. It may be pointed out that 
all these calculations were performed by aver - 
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FIG. 1. Graphic determination of the parameters 
u and £, from the measured quantities w4g9°, “426°, 
and €. 


aging over the finite sizes of source, 8 counter, 
and scattering magnet and over the energy and 
by taking into account Coulomb corrections. A 
complete description of these caiculations and of 
the experimental details will be given elsewhere. 

In conclusion the following ratios of matrix 
elements (compared to the standard matrix ele- 
ment B;;) can be inferred: 


Y =0.42+0.1, 
u=-0.01+ 0.02, 
x =0.05+ 0.05. 


The absolute value of B;; can be derived from 
the ft value by taking into account the results 
for Y, u, and x and is found to be | B;;|? 
=(8.0+0.8)x107° sec“. Thus all matrix ele- 
ments are known. It turns out that the B;; matrix 
element is much larger than the momentum-type 
matrix elements /r and fia xf and is of the same 
order of magnitude as the relativistic matrix 
element f/i@. Moreover it is found that fia is 
in phase with B;;. The unusual enhancement of 
Bjj can be explained by the “selection rule effect” 
as discussed by Morita and Morita’ and by 
Kotani.” 

The results obtained here can be illustrated 
by comparing them with some special cases in- 
vestigated by Kotani.” In Fig. 2, wcosé@ is plotted 
as a function of cos@. Curve I corresponds to a 
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FIG. 2. Comparison of experimental and theoretical 
values of wcos@. Curve I: “unique” transition (only 
Bij contributes). Curve II: modified Bij approxima- 
tion. Curve III: intermediate case. Curve IV: can- 
cellation of matrix elements. 


“unique” decay (Bj; contributes only) and curve 
IV to a cancellation of matrix elements as in the 
decay of RaE. Another curve shows the “é ap- 
proximation” which is valid in most first for- 
bidden decays. The experimental results clearly 
exclude all these possibilities and are repro- 
duced only by the “modified Bij approximation” 
which requires large B;; but Y #0. 

We would like to thank Dr. H. Appel and Mr. 
Blatter for their help in setting up the equipment 
and taking the data and Dr. Morita and Dr. Kotani 
for communicating their results to us prior to 
publication and for valuable discussions. We 
are indebted to the Bundesministerium fiir 
Atomkernenergie und Wasserwirtschaft for sup- 
porting this work. 








'M. Morita and R. S. Morita, Phys. Rev. 109, 2048 
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(1958). 
*T, Kotani, Phys. Rev. 114, 795 (1959). 
34, Schopper, Phil. Mag. 2, 710 (1957); and Nuclear 
Instr. 3, 158 (1958). 
‘The notation of Kotani? is used throughout this paper. 
Swe are greatly obliged to Professor R. M. Steffen 


for enabling us to use his results prior to publication 
and for helpful discussions. According to his meas- 
urements €(p?/W)~!=-0.09 + 0.001 at W=4.3. 

6A similar curve could be derived from the shape 
factor C(W). However, we did not add this curve as C 
is not known with reliability. 





CIRCULAR POLARIZATION OF INTERNAL BREMSSTRAHLUNG ACCOMPANYING 8 DECAY 


S. Galster and H. Schopper 
University of Mainz, Mainz, Germany 
(Received February 23, 1960) 


As a consequence of parity nonconservation in 
weak interactions the particles emitted in B de- 
cay possess definite helicities. This is true 
even for second order processes. Ina previous 
paper’ it was shown that the photons emitted 
together with electrons and neutrinos are cir- 
cularly polarized. For the K capture of A*’ com- 
plete polarization (with positive helicity) was 
found? for the internal bremsstrahlung. The 
circular polarization of photons with negative 
helicity accompanying the 8 decay of Y®, how- 
ever, showed some deviations from the theoreti- 
cal expectation." This might be explained by 
systematic experimental errors or by the cir- 
cumstance that the decay of Y® is unique for- 
bidden whereas the theory has been developed 
for allowed transitions only. In order to clarify 
the situation, the circular polarization of internal 
bremsstrahlung for different types of 8 transi- 
tions (allowed, first forbidden, unique forbidden) 
was investigated and all corrections were de- 
termined as accurately as possible. 

The experimental setup was very similar to 
that used previously." Compton forward scatter - 
ing from magnetized iron was used to measure 
the circular polarization. The inner diameter 
of the scattering magnet (18 cm) and, at the same 
time, the length of the central lead absorber (14 
cm) were increased in order to achieve a more 
effective absorption of the direct radiation. The 
electrons were stopped in paraffin so that the 
contribution from external bremsstrahlung was 
kept low. No correction was applied as this in- 
fluence is negligible compared to statistical 
errors. 

The difficulties in measuring the circular 
polarization of internal bremsstrahlung are 
twofold. The yield is very small (~10-* photon/ 
8 decay) and thin sources must be used to avoid 


external bremsstrahlung produced in the source. 
This requires long counting times in order to 
get sufficient statistical accuracy. Secondly the 
spectrum of the bremsstrahlung drops very 
rapidly with increasing quantum energy. This 
circumstance, peculiar to bremsstrahlung, 
makes some corrections much more important 
than in other experiments. For example, the 
scattered photons degraded in energy have to 
compete with the direct radiation which pene- 
trates the central lead absorber and has higher 
energy and hence a larger primary intensity. 

In calculating the y polarization from the meas- 
ured change in counting rate (produced by re- 
versing the magnetization), corrections were 
applied (a) for the elastic scattering of y rays 
from the iron cylinder, (b) for Compton scatter - 
ing from the copper of the magnetization coil, 
and (c) for photons penetrating the central lead 
absorber. These corrections are small for low 
photon energies but amount to about 30% at 1.5 
Mev. Apart from these corrections the calcula- 
tions were performed in the usual way.* 

In order to avoid the production of external 
bremsstrahlung the sources and their backings 
had to be thin. Hostaphan foils (2 mg/cm’) were 
used as backings. The sources were thinner 
than 2 mg/cm? and possessed activities of a few 
mC. A special procedure had to be used to pre- 
pare RaE sources.* RaE was deposited on Ni 
foils from RaD solution. As together with RaE 
some RaF which has a weak y transition is de- 
posited, a further separation was necessary. 
This was achieved by evaporating the RaE ina 
hydrogen stream and condensing it on cooled 
Al foils (~10 mg/cm?) used as target backings 
in this case. Due to the short half-life of RaE, 
several sources had to be produced. 

The experimental results for the three decays 
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FIG. 1. The circular polarization P of internal 
bremsstrahlung as a function of the average photon 
energy k for the decays of P™ (allowed), RaE (first 
forbidden), and Y* (unique forbidden). The full and 
dashed curves are calculated for allowed decays with 
d=-1, whereas for the long dashed curve of RaE the 
value d= -0.62 was used. 


investigated are shown in Fig. 1. Different 
photon energies were observed by selecting the 
pulse heights with a single-channel discrimina- 
tor. In order to compare the measured polari- 
zations with theory, the polarization was cal- 
culated from the formula 


P =dJ,(k)/[J,(k) +J,(k)], 


which has been derived by several authors.°~’ 


J, and J, are functions of the photon energy k and 
are obtained by an integration over the electron 
spectrum. The quantity d is a combination of 
coupling constants and has for allowed transi- 
tions the value -1 for a V-A interaction with 
maximum violation of parity conservation. The 





polarizations plotted in the figure were evaluated 
by averaging over the finite energy resolution 
of the scintillation counter and by taking into 
account the presence of the Compton peak in the 
pulse-height distribution of a monoenergetic 7 
ray. Of course, these corrections must be 
applied to the measured polarizations and not 
only to the intensities. Details of these calcula- 
tions will be published elsewhere. 

The polarization of the allowed decay of P* is 
in good agreement with the theoretical predic- 
tion, confirming nonconservation of parity ina 
second order process. As the theory for forbid- 
den decays has not yet been developed, we com- 
pare the experimental results for Y® and RaE 
with polarizations calculated for allowed decays. 
It turns out that the polarization of the unique for- 
bidden decay of Y®° is in agreement whereas the 
results for RaE are consistently lower. Similar 
results were obtained for the longitudinal elec- 
tron polarization. This similarity, however, is 
not trivial as the electron and y polarizations 
need not be proportional to each other. Since a 
detailed theory is lacking, we tried to fit the 
results for RaE by an allowed curve, adjusting 
the parameter d, and obtained d = -0.62+ 0.08. 
This value is somewhat lower than that inferred 
from the electron polarization.® 

We are indebted to Professor Strassmann for 
enabling us to produce the RaE source in his 
institute. It is a pleasure to thank Professor 
Bauer for putting at our disposal the electronic 
computer Zuse. We are obliged to the Bundes- 
ministerium fiir Atomkernenergie for supporting 
this work financially. 





‘H. Schopper and S. Galster, Nuclear Phys. 6, 125 
(1958). 

2G. Hartwig and H. Schopper, Z. Physik 152, 314 
(1958). 

3H. Schopper, Nuclear Instr. 3, 158 (1958). 

‘4 detailed description of this method will be pub- 
lished by Dr. Herrmann and Mr. Roder. It is a 
pleasure to thank Mr. Réder for his help in producing 
the RaE sources. 

°G. W. Ford, Phys. Rev. 107, 321 (1957). 

SA. Pytte, Phys. Rev. 107, 1681 (1957). 

"L. Spruch and W. Gold, Phys. Rev. 113, 1060 
(1959). 

§P.) =-(0.75 + 0.02)(v/c); see H. Wegener, Z. 
Physik 154, 553 (1959). 
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(a, t) REACTIONS NEAR Z =28T 


J. L. Yntema 


Argonne National Laboratory, Lemont, Illinois 
(Received February 10, 1960) 


Observation of the charged-particle spectrum 
resulting from the bombardment of vanadium by 


43-Mev alpha particles shows that there 


appreciable cross section for the (a, ¢) reaction 

at small angles. It therefore seemed of interest 
to investigate some (a, ¢) spectra and angular 
distributions in order to study the reaction mech- 


anism. The equipment and experimental 
cedure have been described previously.’ 
spectra of tritons from ,,Mn**(q, ¢),,Fe™®, 


C0™(a, t)pgNi®, .,Cu®(a, t),,Zn™, Rh'%(a, t)Pd’™, 
and Ta*®*(a, t)W*® and of deuterons’ inelastically 
scattered from Rh and Ta are shown in Fig. 1. 

It is strikingly apparent that the Mn**(a, t)Fe*® 





is an 


pro- 
The 








ground -state transition does not proceed with 
measurable intensity. This is precisely what 
would be expected if the (a, ¢) reaction were a 
stripping process, wherein a proton is captured 
by the target nucleus, the angular distribution 
of the tritons acting as an indicator of the angu- 
lar momentum of the captured proton. The 
Mn**(a, t)Fe®® ground-state transition would then 
be forced, by conservation of angular momen- 
tum, to proceed by capture of an f,,. proton and 
would, accordingly, be j-forbidden. Similar 
considerations would lead one to expect a strong 
Co**(a, t)Ni® ground-state transition in agree- 
ment with the experiment. The transition to the 
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FIG. 1. Triton spectra from the (a,¢) reactions on Mn®, Co*’, Cu®’, Rh! , and Ta'®!, and deu- 
teron spectra from inelastic scattering on Rh! and Ta'®! (taken from reference 2). The statistical 
error on the experimental points varies from 3% to 5% in the case of (a,¢) spectra. The triton 
spectra were obtained at 15° lab, the deuteron spectra at 30° lab. The energies are the excitation 


energies of the residual nucleus. 
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2* 1-Mev level in Zn® is, however, much 
stronger than the Cu®(a, ¢)Zn® ground-state 
transition. 

A sharper test of the “direct” character of the 
(a, ¢) reaction mechanism can be obtained by a 
study of angular distributions of the tritons. In 
making a detailed study of this matter we cannot 
proceed from the outset by using the simple 
Butler* theory, since one cannot a priori assume 
it applicable, at least without serious modifica- 
tions. One can however proceed empirically by 
measuring the angular distributions in some 
cases in which the 7 value is known with some 
confidence from shell-model considerations. If 
this procedure leads to a consistent picture for 
a number of (a, ¢) angular distributions the ori- 
ginal assumption of a “stripping” type (a, #) 
process will have been confirmed. 


The V™(a, t)Cr®, Fe(a,¢t)Co®, and Co™(a, t)Ni® 


ground -state reactions (if they proceed by strip- 
ping) must involve /=3 transitions in contrast to 


the /=1 transition expected in Cu™(a, t)Zn®. The 


angular distributions for Co*® and Cu® are shown 
in Fig. 2, where differential cross sections are 
given in relative units multiplied by (k? + a)’, 
where k is the momentum-transfer vector and 

a@ is a constant determined from the binding 
energies. Experimentally this is apparently a 
reasonably satisfactory form factor. The curves 
j,(RR) and j,(kR)? (for R=7 fermis) are also 
shown and are seen to resemble, if not fit, the 
angular distributions from Cu® and Co™, re- 
spectively. The angular distributions of V™ and 
Fe*’ are quite similar to that of Co. Since the 
angular shape of the transition to the first ex- 
cited (2*) state of Ni® closely resembles the 


o Fe°®,3.5 Mev GROUP 
x zea®® 6s. 





one to the ground state of Ni® we conclude that 
it has predominantly /=3. Further, since the 
addition of an f,,. proton can only reach the 
ground state of Ni®, it must involve the trans- 
fer of an f,, nucleon. The angular distribution 
of the Mn**(q, ¢)Fe®* reaction to the 840-kev 2* 
level clearly has an /=3 shape. The transition 
to the group observed near 3.5 Mev in Fe*™ is 
shown to have a considerable /=1 admixture. 

It is clear that more (a, ¢) angular distribu- 
tions must be measured in order to refine the 
kind of analysis outlined above. The accumulate 
evidence so far seems to point strongly toa 
process of the “stripping” type as the reaction 
mechanism. Therefore it seemed of interest to 
analyze a few of the gross-structure groups ob- 
served in the Mn®®, Fe®’, Co®®, and Cu® re- 
actions to see whether an analysis of the type 
used by Schiffer et al.* for the gross structure 
observed in (d, p) reactions would give results 
consistent with shell-model expectations. The 
result for the Mn**(a, ¢)Fe*®® reaction to the 3.5- 
Mev group is in excellent agreement with ex- 
pectations. A preliminary analysis for the other 
gross structures was also consistent with quali- 
tative shell-model predictions. It seems there- 
fore plausible that the gross structure in (a, f) 
reactions may be interpreted in terms of single- 
particle states of the proton. 

The gross structure of the triton spectra from 
Pd’™ and W*® are seen to be very similar® to 
that of the inelastic deuteron spectra from Rh’® 
and Ta’. This suggests that an interpretation 
in terms of single-particle excitation for the 
gross structure observed in the spectra of in- 
elastically scattered protons,® deuterons, and 


.60 3 
Ni -$. 


nie? 1.33 MEV 


1oor : : E ” FIG. 2. Angular distribu- 
F — 4, (6R) F o Fe°® s40KEv tions of Co*% (a, ¢)Ni® ground 
*. = op, ws — i, (KR ;* state, Co°%(a, t)Ni® 1.33-Mev 
“s { ° r ‘5 3 level, Mn**(a,f)Fe® 820-kev 
Ps “eo*” r he level, Mn*°(a, ¢t)Fe® 3.5-Mev 
= r ©.%o , s we as group, and Cu®(q, t)Zn® 
3 . . ground state. The angles are 
= 10F x ‘ Fee A xx in the center-of-mass system. 
b F {* F Be . The relative cross sections 
2 ro | " F o* have been multiplied by 
~ ax x x a e > 








2 
7” 
- 
7~ 
” 
a - 
-, 
-— 
T 





(k?+a@?)?, where k is the mo- 
mentum-transfer vector and 
@ aconstant. The curves 
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alpha particles’»* may be possible. 

Iam indebted to Dr. M. H. Macfarlane, Dr. 
B. J. Raz, and Dr. B. Zeidman for discussions 
and suggestions in connection with this work, 
to W. J. O’Neill for technical assistance, and to 
W. Ramler and the cyclotron group for their 
cooperation. 
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POSSIBILITY OF A TEST OF THE CONSERVED VECTOR CURRENT THEORY 
IN THE A =8 POLYAD* 


Hans A. Weidenmiillert 
Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, California 
(Received February 23, 1960) 


The measurement of the 8-a angular correla- 
tion function in the beta decay of Li® and B® has 
been proposed’ as a test of the conserved vector 
current theory (C.V.C.T.) of beta decay.? In 
order to interpret recent experimental informa- 
tion,® a knowledge of the unknown M1 and E2 
transition probabilities for the transition from 
the first J=2, T=1 state in Be® to the first ex- 
cited state (J=2, T=0) is necessary.‘ It is 
especially interesting to know how much of the 
M1 transition probability in question is due to 
the spin part of the magnetic moment operator 
since only that part gives an enhancement of the 
corresponding second forbidden matrix element 
in beta decay, compared to the prediction of the 
Fermi theory (F.T.). 

This Letter reports on an attempt to determine 
the quantities of interest by an intermediate - 
coupling calculation.* The over-all agreement 
of Kurath’s® results with the experimental level 
scheme in the A =8 polyad, his calculation of the 
M1 widths of the 17.6-Mev state’ in Be®, and 
especially the close agreement of his calculation® 
of the magnetic moment of Li® with observation 
seem to encourage such an attempt. 

Using the standard notation as defined in 
reference 6, we carried through the calculation 
for force mixtures of the Serber, Kurath, and 
Rosenfeld types, for values of L/K between 5.3 


and 8.3 and of a/K between 0 and 5. We calculated 


the energy spectrum for J=2 states, the mag- 
netic moment of Li®, the log/t of the beta decay, 
the M1 and the £2 transition amplitudes, and 
corrections to the magnetic moment and the M1 
amplitude. These corrections are due to the 
presence of spin-orbit and exchange forces and 
have to be introduced in order that the Hamil- 
tonian be gauge-invariant.®° It-is hoped that 
the calculation of these corrections provides an 
order of magnitude estimate of unknown terms 
in the magnetic moment operator, such as ex- 
change moments, etc. 

The relative position of the low-lying J=2 
states obtained with the three different types of 
forces and different values of L/K differ very 
little from each other over the whole range of 
a/K if one adopts the scale parameter K properly. 
The attempt to fit simultaneously the magnetic 
moment and log/t fails the Serber forces. This 
is due to the fact that these forces produce a 
very strong admixture of [31]D states to the 
J=2, T=1 state (which in L-S coupling is [31]*P). 
These states contribute considerably to the mag- 
netic moment and very little to the beta transi- 
tion. This result is in agreement with the ob- 
servation that one cannot reproduce the level 
scheme of Li® using Serber forces. We will 
therefore not consider Serber forces in the fol- 
lowing. 

The spread of the other curves obtained for p, 
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log ft, and the M1 amplitude when plotted versus 
a/K is surprisingly small and confirms the re- 
liability of the estimate. Since the operator for 
the beta decay seems to be much better known 
than the magnetic moment operator, we tried 

to fit" log ft = 5.7 by choosing a/K properly.” In 
this way one obtains for each set of parameters 
a prediction of the magnetic moment. The devia- 
tions of these numbers from the observed value 
exceed in no case 20%. The range of values of 
a/K for such a fit is 1.25,...2.50, in agreement 
with Kurath’s calculations. From these values 
of a/K and the calculated curves for A(M1), the 
transition strength in the definition of Lane’* as 
shown in Fig. 1, one obtains 0.40 < A(M1) < 0.50, 
or for the M1 width 3 ev <I'(M1) <4 ev. 

Figure 2 shows the composition of the M1 matrix 
element for a typical case, L/K =6.8 and Kurath 
force mixture. In this case, the a/K value which 
fits log ft is 2.00. The matrix element is defined 
so that A(M1) = 3 ||M1||)* =0.5 in this case. Only 
the isovector part of the magnetic moment opera- 
tor, 

it (é)= (eh /2me)[-47 Ta) +37 OW, -u oO), 

v F 4 Zz np 
contributes to the transition, and one sees that 
orbital and spin part give contributions with the 
same sign and similar in magnitude. This state- 





A(MI) 


5}. ——--— SERBER : 
KURATH 
~ L/K=5.3 


------- ROSENFELD 














FIG. 1. The strength of the M1 transition as a 
function of a/K for various choices of parameters as 
indicated in the figure. E in Mev, I'(M1) in ev=2.76 
x107§ xE3 x A(M1). 
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ment holds approximately for all cases calculated, 
The dashed line indicates the value one obtains 
for the M1 matrix element if one neglects the 
anomalous moments, i.e., replaces Ly by 

(eh /2Mc)[-47, 14) - 47, 5(2)]. 

The smallness of the spin-contribution can be 
easily understood by the observation that the 
final state is in pure L-S coupling a [4]#D state. 
For a/K =2.5 the admixtures of other states are 
very small (largest admixture 15% in amplitude), 
Since the J= 2, T=1 states do not contain the 
symmetry class [4], the operator (07) has van- 
ishing matrix elements between any of the T =1 
states and the [4]'D state. This statement does 
not hold, however, for the orbital part of the 
operator. For later use we note that 


3.2< (ILI) oll) <4.18, 


and also that A varies practically as the square 
of a/K and is therefore extremely sensitive to 
this choice. If one chooses slightly higher values 
for a/K, one obtains bigger A’s and a larger 
predicted asymmetry in the 8-a@ angular correla- 
tion. At the same time, however, one would 
predict a larger Gamow-Teller matrix element 
than observed. 








‘ 7 MI-TRANSITION AMPLITUDE 
2.0- 

1.5 

ORBITAL CONTRIBUTION 
7 YF 
—~<— BETA DECAY (F.T) 
o. i 
A SPIN CONTRIBUTION 
0 1 2 3 4 . 


a/K 
FIG. 2. The composition of the M1 matrix element 


as a function of a/K for a typical case: Kurath force 
mixture and L/K=6.8. The dashed line gives the 
matrix element which one has to use in the Fermi 
theory of beta decay; the uppermost line (the total 
M1 amplitude) gives the matrix element in the con- 
served vector current theory. 
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As an order of magnitude estimate of additional 
contributions to the M1 transition amplitude and 
the magnetic moment of Li®, the contributions 
due to the presence of a spin-orbit -coupling 
term in the Hamiltonian were calculated. The 
magnetic dipole operator corresponding to this 
term is given in reference 9: 


ehy (i) r, F, 
i. nt = =, a eS 
Baa” 2Mc a(1 ‘, Nai SS r, i} . 


Here, y=aM(7r") fh? = -§ for a= -2.5 Mev. This 
gives a correction of +10% at a= -2.5 Mev, +25% 
at a= -5 Mev to the magnetic moment and of 

-15% at a=-2.5 Mev, -25% at a=-5 Mev to the 
Ml amplitude. This is within the uncertainty of 
our fit. It changes, however, the above state- 
ment to 1.6 <(¢IIZI|) /{ loll) < 5.0. No attempt was 





(a) Fermi theory: 


w(6, F.T.)=1 -0.008 t+ ol) 


made to estimate the contribution due to exchange 
forces’® since there are radial integrals involved 
on which the expression seems to depend sensi- 
tively and about which no information is avail- 
able. 

For the parameters of interest, calculation of 
the ratio of the transition strengths for the E2 
and M1 transitions gave values smaller than 
6x10. Here, as in the preceding estimate of 
y, we assumed (77) = $1075 cm’. The actual 
value of the ratio of the reduced matrix elements 
depends, however, very sensitively on the pa- 
rameters used; in particular, it can be positive 
or negative. 

With this information we are now in a position 
to discuss the prediction of the 8-a@ angular cor- 
relation functions in B® for F.T. and C.V.C.T., 
assuming zero anisotropy in the Li® correlation 
function. One obtains the following results. 


( IE2\l) 





P 
Ciloll) P,lcoss) aN +2 GER, 





-1.44 


(il ip,[30x, + 3x 0, - 2(c- r)6.,) il) 
(3)¢ ll l+ol) , 


where the reduced matrix elements ( |i O||) are dimensionless quantities. 


(b) Conserved vector current theory: 


W(0, C.V.C.T.)=1-0. 008-7 


Here, W, is the maximum, W the actual energy 
of the electrons in units h=c=m=1, P their 
momentum. In both cases the quantities in 
brackets have a small energy dependence, which 
for practical purposes, however, can be neg- 
lected. 

We know from our calculation that ¢ |l7+o1l) / 
(loll) =1.6,...5; ¢11+6.7oll) /{ Holl) =6,... 9.5; 
and (| E211) /|¢ il Z+o ll) ($)|< 0.03. 

If we neglect the E2 terms in both expressions 


( ll+4.7oll) P 
Jee P (egnd) ww, 1+0.56 





(NE2II) 
(3)7( 11244. Tall) §° 





is purely statistical, and includes no estimates 
of possible systematic errors.® This result 
favors the F. T. over the C.V.C.T., though not 
strongly. 

In conclusion it must be stressed, however, 
that the A =8 polyad does not seem well suited 
for a critical test of the C.V.C.T. because of the 
large contribution of the orbital angular momen- 








tum operator to the M1 transition amplitude, and 
because of the well-known fact that our estimate 
of the E2 amplitude may be too small by a factor 
of three or even ten. 

The author would like to thank Mr. E. Nordberg 
and Mr. F. Morinigo for their help in the numer - 
ical calculation. He is greatly indebted to Pro- 
fessor C. A. Barnes, Professor F. Boehm, and 
Professor R. F. Christy for discussions and to 
Professor B. Stech of Heidelberg University 


for the moment, we obtain 
W(6, F.T.)=1+aP,(cosé@), 
where -0.04 <a < -0.01, 
W(6, C.V.C.T.) =1+a@P,(cosé), 
where -0.08 <a < -0.05. 


This is to be compared with an experimental 
finding of a ~(-1.5+3)%, where the quoted error 
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program for diagonalizing the matrices is grate- angular correlation functions . 


’Recently, D. Kurath reported on a similar calcula- 
tion [D. Kurath, Phys. Rev. Letters 4, 180 (1960)]. 
The numbers given in this reference agree with the 
results of our calculation. 
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R PARAMETER IN p - p SCATTERING AT 142 MEV 


L. Bird, D. N. Edwards,* B. Rose, A. E. Taylor, and E. Wood 
Atomic Energy Research Establishment, Harwell, England 
(Received February 23, 1960) 


As part of the Harwell program on nucleon- beryllium analyzers was determined by compari- 
nucleon scattering a measurement has been made son of the counting rates in a given counter tele- 
of the Wolfenstein R parameter in proton-proton scope for the two directions of solenoid field (or 
triple scattering at 142 Mev. proton spin). 

The proton beam whose direction of polariza- The over-all analyzing power P,P, was deter- 
tion was vertical, and for which P, = 0.46+ 0.01, mined by moving each telescope to 6 =0 and re- 
was passed along the axis of a solenoid of such ducing the energy of the beam with an absorber 
integrated field strength that the proton spins so that it corresponded to the energy of protons 
were rotated into the horizontal plane. The di- scattered from the hydrogen target at the appro- 
rection of horizontal polarization could be re- priate angle @. 
versed by reversing the direction of the magnetic The R experiment, particularly when it is self 
field in the solenoid. The beam was then focussed monitoring and performed with a solenoid, is very 
by a quadrupole pair on to a 6-inch long liquid little subject to systematic errors. There is, 
hydrogen target. . however, one serious possibility for error which 

The horizontal transverse polarization of pro- arises from the focussing properties of the sole- 
tons scattered at angle @ in the horizontal plane noid for protons. It is necessary to position the 
(i.e., in the plane containing the direction of solenoid extremely carefully so that the mean 
polarization and the incident vector) was then position of the beam at the hydrogen target does 
measured by scattering in the up-down plane not move appreciably when the solenoid current 
from targets of beryllium of analyzing power P,. is reversed. In the present experiment, the move- 
Protons scattered up and down were detected in ment was not more than 0.2 mm, which called 
two two-counter telescopes placed in coincidence for corrections much less than the final standard 
with a counter immediately in front of the beryl- errors on all the measurements. 


























lium analyzer. Altogether ten counter telescopes As a final check on the systematic errors in 
were used, which allowed measurements to be the measurements, the experiment was repeated 
made at five different scattering angles 6 sim- with an unpolarized incident beam, and the asym- 
ultaneously, in a self-monitoring system. The metries were found to be consistent with zero in 
asymmetry RP,P, in the scattering from the this case. 
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Table I. Measured asymmetries and derived values of R. 











Asymmetry RP,P,; 


P,P, R 





45 


24. 
32. 


54. 
67. 
76. 
84. 
90. 


oor NnN rR AAO 


-0. 
-0. 
-0. 
-0. 
-0. 
-0. 
-0. 
+0. 


0421 + 0.0095 
0351 + 0.0088 
0249 + 0.0044 
0252 + 0.0050 
0187 + 0.0035 
0100 + 0.0043 
0084 + 0.0080 
0054 + 0.0064 


.188 -0.224+ 0.051 
.173 -0.203 + 0.051 
.140 -0.178 + 0.031 
.119 -0.212 + 0.042 
-088 -0.213 + 0.040 
. 068 -0.147 + 0.063 
.059 -0.142 + 0.136 
.049 +0.110+ 0.131 


ooo coo coco co 
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20 40 60 
CENTER-OF-MASS SCATTERING ANGLE 


FIG. 1. R parameter in P-/P scattering at 142 Mev, 
with earlier theoretical estimates by Signell and 
Marshak and by Gammel and Thaler. 


The values of asymmetries, analyzing powers, 
and of R are given in Table I. The results are 
shown in Fig. 1, together with predictions from 
the Signell- Marshak’ and the Gammel-Thaler? 
potentials. These are both in qualitative accord 
with the trend of the experimental points. 

Preliminary results have been presented at 
the Few-Nucleon Conference, London, July, 1959 
(unpublished) and at the Ninth Annual International 
Conference on High-Energy Physics, Kiev, 1959 
(unpublished). 

A full account of the experiment will be pub- 
lished elsewhere. 





"Now at Liverpool University, Liverpool, England. 

‘Pp, S, Signell and R. E. Marshak, Phys. Rev. 109, 
1229 (1958). 

23, L. Gammel and R. M. Thaler, Phys. Rev. 107, 
291 (1957). 





POSSIBLE NEW RESONANCE IN THE 1° -p SYSTEM 


* 
Peter Carruthers 
Laboratory of Nuclear Studies, Cornell University, Ithaca, New York 
(Received February 19, 1960) 


Within the past year a great deal has been 
learned about the interactions of high-energy 
negative mesons with protons. In particular 
two very striking peaks have been found to exist 
in the s~ -p cross section at lab pion energies 
of about 600 and 900 Mev.'~* Angular distribu- 
tions of elastically scattered pions in this energy 
range have been reported.?** However, experi- 


mental data for 7* -p scattering in the same 


energy range are not nearly so complete. Total 
cross-section measurements’»® have indicated 
that in this energy region the 7* -p cross sec- 


tion reaches a minimum at about 650 Mev, 
thereafter rising slowly to the broad peak near 
1.3 Bev.*® 

In the present note we wish to anticipate future 
investigations of the 1 -p system by drawing 
attention to some features of data already avail- 
able which suggest the existence of a new reso- 
nance in the isotopic spin 3/2 (x* -p) state at an 
energy of about 850-950 Mev.” This evidence 
further suggests that the favored state has nega- 
tive parity and total angular momentum 3/2 (D,,). 
This is precisely the state of angular momentum 
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and parity assigned to the 600-Mev resonance 

of isotopic spin 1/2 (the “second” resonance) 

through the work of Wilson,® Peierls,® Sakurai,’® 
and Stein." This assignment for the second 
resonance is essential to the arguments pre- 
sented here and will be assumed throughout. 
With the notation L,y7 (L is the orbital momen- 
tum, J the total angular momentum, T the total 

isotopic spin) we suggest that there is a D,,, 4. 

resonance as well as the already established® 
Dy, y2 Tesonance. Theoretical arguments are 
given for the particular quantum numbers of the 
suggested resonance, and for the different ener- 
gies of the two D,, resonances, at the end of 
this Letter. 

Unfortunately the hypothetical D,, ,,. resonance 
cannot be expected to reveal itself in sucha 
dramatic fashion as have the well-established 
resonances, for the following reasons. First 
of all at such high energies (near 900 Mev) the 
maximum permissible cross section for a state 
having J=3/2 is only 30 mb, and this will be 
further reduced by the occurrence of inelastic 
processes. Further, since the proposed reso- 
nance lies at the foot of the broad peak at 1.3 
Bev, the amplitudes for the states composing 
this latter peak will already be appreciable be- 
fore the D,,.,42 resonance has subsided, pre- 
venting the appearance of a distinct peak, although 
one may hope that more refined experiments will 
reveal some structure in the cross section. 
(Clearly these observations do not diminish the 
importance of such a resonance for the under - 
standing of the dynamics of meson-nucleon inter - 
actions.) Therefore a measurement is proposed 
to detect the characteristic PD interference in 
the recoil proton polarization [see Eq. (3)]. The 
Saclay group® reported 7* - measurements up 
to about 850 Mev. Their data, besides displaying 
sharper and higher 1~ -p peaks, show both a 
lower minimum in the 7* -p cross section than 
previously found’’® and an abrupt rise in o(7* -p) 
in the interval 700-850 Mev. If this curve is 

continued up to the points of reference 6, a dis- 
tinct shoulder appears in the cross section.” 
We now discuss the possible structure of this 
shoulder. 

First consider the remarkable behavior of the 
elastic charge exchange scattering (17 +p-—7° +n) 
in the 500-800 Mev energy range. Throughout 
this interval the charge exchange scattering is 
small and rapidly decreasing, without any sign 
of a peak, although the elastic scattering (1 ~ - p) 
displays a sharp peak near 600 Mev.” This is 
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rather surprising, since the cross sections are 
given in terms of the amplitudes f,, and /,,. (the 
subscript denotes isotopic spin) by 


Oop = Maal? +4If ial? +4Reia faa) 


Coe E2NFaal?+21fyal?-4RE ie Fan)» (2 
[The proportionality factor is the same for (1) 
and (2).] Since f,. is known to vary in a reso- 
nant manner, the smallness and smoothness of 
the charge-exchange scattering requires that the 
interference term Re(f ~ be large and posi- 
tive, a point already noticed by Peierls.'* Since 
we are concerned with total cross sections, only 
states having the same L and the same J survive 
in the interference term. Now if these T =3/2 
states were small and nonresonant, one would 
expect the charge exchange cross section to rise 
abruptly on the high-energy side of the resonance, 
which is not the case. Since D,, ,,. presumably 
is the most important siate in this energy inter- 
val, the simplest way to prevent both the occur- 
rence of a peak in o,, and an abrupt rise above 
the second resonance is to have the D,,_ ,. phase 
shift positive and growing rapidly above 600 Mev. 
(The relative phase of the states D,,, ,,. and 
D yz, y2 is expected to be large and near 90° be- 
tween about 600 and 900 Mev, but it must be 
rather less than 90° for the present argument to 
work.) Let us consider the other interference 
terms. Since the S-wave phase shifts are not 
too large and probably have opposite signs,?**»* 
their interference will not help in explaining the 
behavior of ogx. Next consider the P,,. inter- 
ference; a@,, and a@,, are both small and negative 
(we use QT 2 for the P phase shifts, and re- 
serve bop 27 for the D waves) at 300 Mev,'*® 
and at least a,, is negative at 500 Mev.’* Hence 
this term will either be small or else have the 
wrong sign. The P,, interference term will have 

the wrong sign as long as @,, is positive.’> In 

any event a,, is varying slowly in the energy 
region under consideration. Accurate T =3/2 

D-wave phase shifts have been obtained at 310 

Mev by Foote et al.'® They find that 6,, is posi- 

tive, indicating an attractive D,, ,,. state, con- 

sistent with the present scheme, and a negative 

Dye, y2 Phase (5,,). The solutions f and g ob- 

tained by Willis also have this character.“ 

Crittenden et al.? report a tentatively negative 

(and small) 6,, at 460 Mev. Thus it seems un- 

likely that D,, interference can have the re- 

quired character. There is some evidence’ that 
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the third resonance (T = 4) may be associated 
with an Fy. 2 State. However, since the pro- 
posed resonance is supposed to have an energy 
comparable to the third resonance, it seems 
impossible for the Fy, 2 - Fs, 1 interference 
to be appreciable in the vicinity of. the second 


resonance. 

Further, the size of the 7* -p cross section at 
900 Mev (= 22 mb) is compatible with a dominant 
Dy, State with an inelastic scattering of about 
§mb. (Since the other states have been ignored, 
this underestimates Oj .) This is to be com- 
pared with o;, =2.8+0.5 mb at 500 Mev"* and 
g,,=10.5+ 1.5 mb at 990 Mev.’” Only a com- 
pletely absorptive J=5/2 wave can have as small 
atotal cross section as 22 mb at 900 Mev. 

Now consider the available 7* -p angular dis- 
tributions. At both 500 Mev"* and 1.1 Bev"® the 
most noticeable feature is the very small amount 
of scattering in the backwards direction. If the 
scattering is really strong in a D,, state around 
850-900 Mev, then one should expect some 
1+3cos?@ component to be detectable in the 
backwards hemisphere. Now in fact, at 990 Mev, 
where the “resonance” should still be appreciable, 
there is in fact an encouraging backwards maxi- 
mum.’” Photoproduction angular distributions 
should show signs of this state, but in the total 
cross section the proposed resonance probably 
will not be distinguishable from the third reso- 
nance, which occurs at about the same energy. 

As mentioned before, the detection of a large 
negative (with respect to the normal to the piane 
of scattering) polarization for large scattering 
angles (forward protons) at an energy of, say, 

800 Mev would be good evidence for the picture 
presented above. This rough choice of energy 
would be low enough to avoid complications from 
other states and the appreciable inelastic scat- 
tering which accompany the 1.3-Bev peak. It is 
low enough for the tail of the P,,,,. resonance 
to be large and high enough for the proposed 

Dye ¥2 resonance to be important. Keeping only 
these two states in computing the polarization, 
one has for the polarization of the recoil proton 
{the amplitudes f are [p exp(2i6) - 1]/2ik; 0<p<1} 


P Hdo/dQ 


=2 sin6(9 cos76 - 1)Im(f, 
¥2, ¥2 


¥ . (3) 
f ie wl 

(kis the pion momentum in units of h.) Note the 

large cos?@ term. One can see that the other 

states involved give interference terms which 





are likely to be small, and if anything, tend to 
offset the above polarization. [Note that if 
Dyz, 42 is replaced by D,,,,,. then the right-hand 
side of (3) is generally smaller and more iso- 
tropic since 9cos*@ - 1 becomes 3(1 +cos76). | 
The actual magnitude of P, depends on a,,. The 
proton energy (near 400 Mev) should make the 
polarization experiment feasible. 

The dispersion relation calculations of Stern- 
heimer’® show that the real part of the forward 
scattering amplitude for 1* -p scattering vanishes 
at about 900 Mev pion energy. Such a circum- 
stance is certainly compatible with the existence 
of a dominant resonance state at this energy (if 
the real part of the phase shift is 90°, then f is 
purely imaginary). At this energy the S waves 
and some of the P waves are probably mostly 
absorptive, making no contribution to Re(f), and 
the cancellation of the amplitudes of the remain- 
ing states is quite plausible. 

Peierls’* has interpreted the second and third 
resonances in terms of two-meson states, utili- 
zing ideas from the isobar model.”° For instance, 
for the negative-parity second resonance 
(Dy, ye) the S -P configuration is supposed to be 
important. This model had some difficulty in 
explaining why there was a resonance in the 
Dyz, y2 State but not in the D,,, 4. state. This 
difficulty is resolved if the point of view of the 
present note is taken. 

Finally it is necessary to consider why the two 
D-wave resonances should be so far apart (250- 
300 Mev in the lab; 150-175 Mev c.m.) in energy. 
As might be expected, this behavior is to be 
blamed on the 33 resonance. A partial explana- 
tion may be contained in the observation that in 
the partial wave dispersion relations the cross 
section for the 33 resonance enters with opposite 
signs for the T=1/2 and 3/2 amplitudes. 

The author would like to thank Professor H. A. 
Bethe for a number of enlightening discussions 
of this work, and Dr. A. Erwin for supplying 
some unpublished data. 
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For many years, the question to how to investi- 
gate the behavior of the weak interactions at high 
energies has been one of considerable interest. 

It is the purpose of this note to show that ex- 
periments pointed in this direction, though not 
quite feasible with presently existing equipment, 
are within the capabilities of present technology 
and should be possible within the next decade. 

We propose the use of high-energy neutrinos 
as a probe to investigate the weak interactions. 

A natural source of high-energy neutrinos are 
high-energy pions. Such pions will produce 
neutrinos whose laboratory energy will range with 
equal probability from zero to 45% of the pion 
energy, and whose direction will tend very much 
toward the pion direction. For example, 1-Bev/c 
pions will emit neutrinos with an average energy 
of ~ 220 Mev in such a way that ~ § of the neutrinos 
will fall within a cone of half-angle 7°. For orien- 
tation purposes, the mean decay distance for 
such a pion would be 50 meters. 

The best known source of pions is a proton 
accelerator where the beam is allowed to im- 
pinge on a target. Let us assume that we have 
available a 3-Bev proton beam and 10 000 kilo- 
grams of material for sensing a neutrino inter - 
action. We may then estimate the proton flux 
necessary to produce one interaction per hour 

with a cross section of cm’. To do this, let 
us consider the simple setup shown in Fig. 1. 
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FIG. 1. Proposed experimental arrangement. 


Let J be the number of incident protons per unit 
time, and let, say, //10 charged pions with en- 
ergy 22 Bev be produced at the target. These 
pions emerge in a cone of about 45° half-angle, 
or in about 2 steradians of solid angle. We now 
let them travel for a distance of 10 meters before 
hitting a 10-meter shielding wall in front of the 
detector. Approximately 10% of the pions will 
decay with an average neutrino energy of about 
400 Mev. Each square centimeter of detector 
subtends a solid angle of }x10~° steradian. 
Hence, the high-energy neutrino flux at the de- 
tector is (¥/)(}x10~*)(4)(4) =1x107°/. If there 
are 10000 kilograms of detector present, the 
number of events per unit time is given by 


N ~(107)(6 x 10*°)(10~*/)0 = 6 x 107Jo. 
For an intensity of [=5 x10 protons/sec =1.8x1" 
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protons/hr, the high-energy neutrino flux is 


~5§000 neutrinos sec“ cm~*. With a cross sec- 


tion o~107**® cm?, the number of counts is N ~1 
per hour in 10000 kg of detector. The estimate 
here given is for neutrinos from high-energy 
pions. There is, as a matter of fact, a much 
greater flux of lower energy neutrinos from lower 
energy pions. However, because the neutrino 
cross section decreases rapidly with decreasing 
energy, the rate is not likely to be improved by 
more than a factor of two. 

This estimate places the experiment outside 
the capabilities of existing machines by one or 
two orders of magnitude. Optimistic estimates 
for accelerators which are currently under con- 
struction, namely the 3-Bev machine at Prince- 
ton and the 10-Bev machine at Argonne, indicate 
that the experiments may be barely feasible in 
the near future. However, for really quantitative 
experiments it will be necessary to use high- 
intensity machines such as the FFAG machine 
proposed by MURA or the 10-Bev linear proton 
accelerator discussed by Blewett at Brookhaven. 
In these machines, one hopes to attain a beam 
intensity of the order of 10** protons/sec at an 
energy of about 10 Bev. 

The higher energy of the primary beam of pro- 


tons makes the experiment easier because of the 
increased multiplicity of pions, the more con- 
centrated forward distribution of the pions, and 
the increased cross section for neutrino reac- 
tions. Balanced against these is the fact that the | 
percentage of higher energy pions that decay in 
10 meters is smaller. The net result is likely 

to give a counting rate per primary proton that 
probably increases more than linearly with the 
primary proton energy. 

Thus, a high-intensity 10-Bev proton machine 
with a beam intensity ~10** protons/sec may 
give a counting rate of more than 10° per hour, 
using the experimental setup described above. 

If that proves to be the case, it is perhaps de- 
sirable to have magnetic lenses to analyze and 
focus the pions so as to obtain more monoener - 
getic neutrino beams. 

I would like to express my gratitude to Dr. T. D. 
Lee and Dr. C. N. Yang for many stimulating dis- 
cussions which led to the above proposal. 

Note added in proof. The author’s attention has been 
called to a somewhat related paper which has just 
appeared: B. Pontecorvo, J. Exptl. Theoret. Phys. 
(U.S.S.R.) 37, 1751 (1959). 
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The weak interaction so far most extensively 
studied is 8 decay, in which the momentum trans- 
fer is of the order of a few Mev. In yp decay and 
capture, momentum transfers of the order of 
100 Mev are involved. In the theory of these 
processes, because of the limited region of mo- 
mentum transfer studied, the phenomena can be 
described by a few parameters usually called 
coupling constants. For larger momentum trans- 
fers it is obvious that the weak interactions can- 
not continue to be described by these constants, 
because of the clothed structure of the nucleons 
due to the strong interaction, and also because 
of the reasonable expectation that the weak inter- 
actions, even without the interference of the 
strong interactions, may not be of the simple 
four-spinor product form in Fermi’s theory. 


it is easy to see that v, and v, are the same par- 






In the preceding Letter,’ Schwartz points out 
that the neutrinos from the decay of high-energy 
mesons can be used to study weak interactions. 
We have investigated the theoretical implication 
of such possible experiments. Efforts are made 
to separate and dissociate the inferences that 
can be drawn from different assumptions con- 
cerning the weak interactions. In this Letter we 
report briefly on this work. 

1. The identity of the neutrinos. In the pro- 
cesses 





at—=pt+ Pas (7 decay) (1) 
(u capture) (2) 


Z~(Z-1)+e* + vg, (a* decay) (3) 


LL +p—-n+ Vp, 
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ticle. Experimentally it is known that v, and v, 
both have helicity -1. It is simplest to assume 
that v, and vy, are also the same particle. How- 
ever, a test of this assumption is clearly desir- 
able. To obtain such a test it is necessary to do 


some kind of capture experiment on the neutrinos 


or antineutrinos. For example, if v, and vy, are 
different particles, then the reaction 


n+v,~pr+e” (4) 


does not occur. 
2. Conservation of leptons. The conservation 
of leptons can be studied with neutrino capture 





experiments. For example, if both 
v,+p—A°+e* (5) 
and 
D,+p~-A°+e* (6) 


occur, there would be a violation of lepton con- 
servation. While it is possible to study lepton 


conservation by helicity measurements, neutrino 


capture experiments seem to be the most direct 
and clean cut for such purposes. 

In the rest of this Letter we shall assume that 
V, = V,=Vv,=v and that the conservation of leptons 
holds. 

3. Possible existence of a neutral lepton cur- 





rent. By using high-energy neutrinos it becomes 


possible to study whether reactions such as 
v+p—-v+p and v+n—-vin 


exist or not, and if they exist, whether there is 


any similarity between these “neutral lepton cur- 


rents” and the electromagnetic field. [See also 
the discussion in Sec. 8. ] 


4. Point structure of the lepton current. In the 





present theory of 8 decay, yu capture, etc., one 
assumes that all the weak reactions that contain 
both leptons and heavy particles can be repre- 
sented by an effective Lagrangian of the type 


4 
£ oe? 2, Ly, i, Ce) J, "GH, "I, (7) 
where 
j, () - -ilv,"y47, (1 +7. 1, (8) 
j,'e)= ily Nyy, (1 +y,)¥)1, (9) 


I stands for either e™ or yu’, y, and y; are the 
field operators for v and/, and J,(x) and J)’(x) 
are operators that act on heavy particles (in- 

cluding the pions, K mesons) only. Because of 
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eff, we have 


‘= tT 
Jy mdy's (10) 


the Hermiticity of £ 


where 


n, =+1 for A=1,2,3, 
lee for A=4. (11) 


The nature of J, and/J,’ is known so far only 
in the nonrelativistic region. In the low-energy 
limit, the matrix elements of these heavy-par- 
ticle current operators in the case of 8 decay 
are of the form 


(Pld, In)=(G/N2)u Ty gy, (Gy - Garg» (12) 


v7 &47%5n 


(nid, "\p)=(i/V2)u Vy gY(Gy" - Gy "vu p> (13) 
where u, and uy are the spinor solutions of the 
free Dirac equations with the same 4-momenta 
as the physical neutron and proton; Gy and G, 
are the Fermi and Gamow-Teller coupling con- 
stants. Due to the presence of strong interac- 
tions it is expected that in the high-energy re- 
gion (12) and (13) do not hold. But one expects 
(7), which represents a “point interaction” for 
the leptons, to have a wider range of applicabil- 
ity. 

The mere assumption that in the effective 
Lagrangian (7) the lepton current acts only at a 
single space-time point introduces rather strong 
restrictions on the forms of the cross sections 
for all neutrino and antineutrino reactions. For 
example, in either 


v+n—-p+l +pions (14) 
or 
v+p—n +1* + pions, (15) 


suppose one measures in | the laboratory system 
the incoming momentum k,, and the outgoing 
lepton momentum ky and does not measure the 
other kinematic quantities describing the reac- 
tion. The experimental cross section is then in 
a general case a function of the three real vari- 
ables k,, k;, and @ (= the angle between kj and 
k, y)- Independently of the form of J,, assump- 
tion (7) restricts this function to a sum of three 
structure functions each of which has an unknown 
dependence on only two real variables: E =k,-*; 
and P= lk, - kyl, which represent the energy 
transfer and the magnitude of the momentum 
transfer between the leptons and the strongly in- 
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teracting particles. More explicitly, assump- 
tion (7) implies? that the cross section for (14) 
[also for (1 5)] is of the form 


do= dk d(cosé Mank )-*k [(k, +k)? - P?] 
x [xA+x°'A_ +B}, (16) 


where Xx =(k)+k,-P)(kj+k,+P)~* andA,, A_, 
and B are functions of E and P only. 

To test the validity of (16) it is not necessary 
to perform a detailed experiment for specific 
values of E and P. One could perform a capture 
experiment with a neutrino beam with a known 
spectrum J(k ,,)dk, and measure for each event 
the values of x, P, andE. If N(x, P,E)dxdPdE is 
the number of events, then (16) implies that 


Ie we, P,B)=[A, Age Age 1-0), (17) 


where A; (t=1,2,3) are functions of P andE. In- 
tegrating (17) over P and E, one obtains 


it Tae (a, +a,x+agx")(1-x), (18) 
where @,, @,, @, are numerical constants, and 
the sum extends over all events with fixed x. To 
test the validity of (18), less than a thousand 
events could be enough provided that they do not 
cluster around one value of x. 

5. Universality of weak interactions involving 
e* and u*. Neglecting the mass of the » meson, 
the assumption of the universality of the weak 
interactions implies equal differential cross sec- 
tions for u* and e* production and for y~ and e~ 
production. If the mass of the » meson is not 
neglected, comparison should be made between 
u and e production processes in which the en- 
ergy transfer E and the magnitude P of the mo- 
mentum transfer from the leptons to the strongly 
interacting particles are fixed, provided the 
point structure of the lepton current discussed 
in Sec. 4 is valid. 

§. S-symmetry. In (16) the structure func- 
tions determined by using v are related to the 
appropriate matrix elements of J, in Eq. (7) 
while those determined by using v are related to 
that of J,’. Thus, unless J, and J,’ obey some 
further symmetry property the reaction rates of, 
&g., (14) and (15) in general are not related to 
each other in any simple way. A symmetry 
which will be called S-symmetry is of a type so 
a8 to link J and J’. 

To explain the meaning of S-symmetry we ob- 





serve that by using (12) and (13) at the low- 
energy limit it is readily verified that J, and 
J,’ satisfy the following relation: 


Is -1 
J\'=SJ,S", (19) 


where S is the product of a 180° rotation along 
the y axis [exp(iz/,)] in the isotopic spin space 
multiplied by the time-reversal operator T, 
i.e., 


S= [exp(énl,) ]T. (20) 


Condition (20) will be defined as the condition 
for S-symmetry.* If it is satisfied then there 
should exist identities among the matrix ele- 
ments of J andJ’ and consequently relations be- 
tween reaction rates caused by v and v. For ex- 
ample, if the mass of the lepton is neglected, 
the structure functions [defined in (16)] for v and 
v processes are related to each other by 


4.) =@)-, 4) =4)-, (B) =(B)5- (21) 


Experimentally it is not known at present 
whether S-symmetry is satisfied at energies of 
the order of 100 Mev and up. Theoretically it 
has been customary* to assume bare particle 
universal Fermi interactions for all Fermions. 
If such an assumption is made, S-symmetry 
follows. 

7. Conserved vector current and proportion- 
ality with the electromagnetic current. The 
heavy-particle current J, can be written as a 
sum, 








J,=V, +A), (22) 


where Vy and A, are, respectively, its vector 
part and axial-vector part. Recently, Feynman 
and Gell-Mann® proposed that the vector part 
V, satisfies the conservation law 


eV, /e, =0. (23) 


Furthermore it is proposed that V, is equal to 
the corresponding isotopic vector part of the 
electromagnetic current times a constant (here- 
after called the proportional vector current hy- 
pothesis). A sensitive test of the proposal can 
be given by studying reactions such as 


v+n—-e +p, (24) 
p+p—e* +n, (25) 
and comparing them with existing data on elec - 


tron scattering by nucleons® at the same momen- 
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tum transfer to the nucleons. 

Making the proportional vector current hypo- 
thesis, the electron scattering experiments® 
show that V, is proportional to 
FQ@’)=[1+ 540°)" 
where q’ is the invariant four-momentum trans- 
fer squared. We determine the constant of pro- 
portionality by the condition Gy=10-°/M*. For 
orientation purposes we assumed that A, is 
given by 


(PIA, In)=1.2 GyPGN2)u Ty 7 5ty 


~ @ =0.8x1072° cm, 


and calculated the cross sections for (24) and 
(25). The results are exhibited in Fig. 1. 

8. Possible existence of a weakly coupled 
Boson W*. The question whether the weak inter- 
actions are “transmitted” by a Boson field was 
already discussed in Yukawa’s original work on 
the meson. If such a Boson field W~ exists, it 
must have spin 1, and one can also conclude that 
its mass my, is 2m. The nonlocality of the 
weak interactions implied by the finite mass of 
the transmitting field W has been discussed be- 
fore.’ For yu decay, the change in the Michel 
parameter p is given by 





ms =i é 
p 0.75=3(m /my) : 


A value of my,>m x is thus consistent with exist- 
ing experiments. 
The coupling of the W* to the leptons is char- 
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FIG. 1. “Elastic” neutrino cross sections. The 
dashed line represents the limit of o ask, ~-@. 
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acterized by a coupling constant, 


g?/4n =(nV2)™ Gym 7 <6.4x107". 


The decay rate of W~e+ v is found to be 


= 3 ee 
Awees p> Cy My (60v2)", (26) 
and the ratio of the decay rate of W~+ v to that 
of W-e+vis 


2 3 
av (S+u )/(l+v ) - (27) 


where v Mm is the velocity of the muon in the rest 
system of W. Decays of W into pions are also 
possible. The lifetime of W is <10~*" sec. 

While the existence of a W* as a virtual par- 
ticle does not change in any essential way the 
considerations of the previous sections, its pro- 
duction makes possible a much higher cross sec- 
tion for neutrino reactions through the pair crea- 
tion of W* and J” in the Coulomb field of a target 
nucleus: 


v+Z—-W4+1° +Z. (28) 


The cross section for (28) is large compared 
with those without W production. At values of 

k |, >>my’/2q,=2 Bev, where g,=h/nucleon radi- 
us, the cross section is given by 


o= (61V2)"*(137)-?Z7G_ [In(2k qo/m,,?)]*, (29) 


which, for Z=26, is of the order of 107° cm’ 
The W* produced are easily identifiable through 
their decay products e+v or w+v. If experi- 
mentally no W* is found, it would be possible to 
set a lower limit on the value of my: 

The existence of the intermediate Boson W* 
has® been discussed in connection with the ques- 
tion of the absence of u.*~e*+y. On reasonably 
sure grounds one could conclude that the exist- 
ence of W* requires that in the notation of (1), 
(2), (3), v,#v5. In other words, the existence of 
W would imply that reaction (4) does not occur. 

For processes without W production, the dis- 
cussions of Sec. 4 above concerning a point struc 
ture of the lepton current remain unchanged, but 
the current J, now includes the effect of the pro- 
pagator of W. Furthermore a corresponding 
change in the proportional vector current hypo- 
thesis is necessary. 

The question of a neutral W° will not be ex- 
amined here. 

9. Interactions with extremely large momen- 
tum transfers. For momentum transfers of the 












order of (Gy)-’*~300 Bev/c, quantities other 
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than the lowest power of Gy become important. 
For example, for e~-+v-yu +v, the cross sec- 
tion predicted by the lowest perturbation formula 
approaches the limit set by the unitary condition 
as the neutrino momentum in the center-of- 
mass system, p,~300 Bev. For processes in- 
volving pion clouds, high-order terms in Gy 
presumably also become important as p y~300 
Bev. The description of the weak interaction at 
such high energies would seem to need ideas 
radically different from the current picture. 

We wish to thank Professor M. Schwartz for 
many discussions. 





"Work supported in part by the U. S. Atomic Energy 
Commission. 
iM. Schwartz, preceding paper [Phys. Rev. Letters, 


4, 306 (1960). 

2In Eq. (16), we neglect the mass of the lepton. For 
reactions in which the mass of the lepton may not be 
neglected, equations similar to but slightly more com- 
plicated than (16) result. For reactions in which 
strange particles are produced, (16) remains valid. 

3By using the CPT theorem it is easy to see that the 
S -symmetry is closely related to the classification of 
weak interactions given by S. Weinberg, Phys. Rev. 
112, 1375 (1958). 

‘See, e.g., M. L. Goldberger and S. B. Treiman, 
Phys. Rev. 111, 354 (1958). 

5R. P. Feynman and M. Gell-Mann, Phys. Rev. 
109, 193 (1958). 

®See, e.g., R. Hofstadter, F. Bumiller, and M. R. 
Yearian, Revs. Modern Phys. 30, 482 (1958). 

'T, D. Lee and C. N. Yang, Phys. Rev. 108, 1611 
(1957). 

8G. Feinberg, Phys. Rev. 110, 1482 (1958); J. 
Schwinger, Ann. Phys. 2, 407 (1957). 





ENHANCEMENT OF BREMSSTRAHLUNG PRODUCED BY 575-Mev ELECTRONS 
IN A SINGLE CRYSTAL OF SILICON* 
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(Received February 25, 1960) 


Experiments have been performed’ ° to ob- 
serve the expected enhancement of the brems- 
strahlung produced by high-energy electrons 
when they pass through a single crystal.* The 
theory of Uberall has been critically examined, 
and improved calculations have been made by 
Schiff.° The purpose of this report is to present 
the preliminary results of the revised experi- 
ment performed at Stanford, using a more sen- 
sitive technique of detection of the enhancement. 

A single crystal of silicon, 0.013 in. thick and 
cut along the [110] plane, was aligned normal to 
the x-ray beam as defined by a 0.025-in. colli- 
mator, by the Laue back-reflection method. The 
crystal in its holder was mounted in a remote- 
controlled double-goniometer in which it could 
be rotated accurately about two coplanar per- 
pendicular axes normal to the electron beam of 
the accelerator. The normal position of the [110] 
plane of the crystal with respect to the electron 
beam was established by using optical methods 
to transfer the x-ray alignment. Electrons of 
315 Mev produced bremsstrahlung in the single 
crystal and were then deflected. The deflected 
electron beam was monitored by a standard 
secondary emission monitor.* The brems- 






strahlung was monitored by a double -ionization 
chamber consisting of two identical ionization 
chambers constructed in one envelope, each 
preceded by appropriate beryllium converters. 
This unit was mounted accurately so that the 
axis of both the ionization chambers was in line 
with the incident beam direction. Continuous 
flow of purified hydrogen was maintained in the 
chamber at a constant pressure. The sensitivity 
of the first ionization chamber was a maximum 
around a photon energy of 30 Mev. The second 
ionization chamber looked at pairs which on the 
average have penetrated a larger range of beryl- 
lium and are therefore produced by a relatively 
harder component of the bremsstrahlung. Since 
the enhancement is large in the soft component 
of the bremsstrahlung, the relative increase of 
the charge in the first ionization chamber over 
that in the second ionization chamber would be 
quite pronounced. The charges collected in the 
two ionization chambers were fed separately to 
two de amplifiers, the output of which was 
balanced using two potentiometers. The increase 
of the measured ratio as the crystal plane ap- 
proached the correct angle with respect to the 
electron beam, would indicate the enhancement 
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FIG. 1. The ratio of the charges, a measure of the enhancement of the soft component of 


the bremsstrahlung, is plotted against the angle between the direction of the incident elec- 
tron beam and the [110] axis of the Si crystal, in which the bremsstrahlung is produced, 


about the horizontal and vertical axes of rotation. 


6p is the angle of rotation of the crystal 


about the horizontal axis, and 6, is the angle of rotation about the vertical axis. 


of the soft component of the bremsstrahlung. 
Thus, this method provides a more sensitive 
way to detecting the enhancement than that used 
previously,’ using the photoproduction of mesons. 
The ratio of the charges, a measure of the 
enhancement, was observed fcr various angles 
as the crystal was rocked back and forth about 
the vertical and horizontal axes. An angular 
grid of about +16’ of 6, the rotation about the 
vertical axis, and about + 48’ of 6), the rotation 
about the horizontal axis, was covered in steps 
of 4’. The data are shown in Fig. 1 in a three- 
dimensional plot of the ratio vs 6,, and 6p. It is 
clear that the soft component of the brems- 
strahlung is definitely enhanced. The theoreti- 
cally predicted dip in the center of the peak 
could not be observed. Calculations are in prog- 
ress to estimate quantitatively the agreement 
between the theory and the experimental data; 
these calculations depend in detail on the pre- 
cise energy response curve of the instrument 


to photons of various energies. 

I am happy to acknowledge the help of Pro- 
fessor W. K. H. Panofsky, Professor H. C. 
DeStaebler, and Professor L. I. Schiff, and of 
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Theoretical work on nucleon structure’? has 
revealed the important role of the photon -(» -pion) 


In particular the y - 27 vertex is thought to be 

the most important one contributing to the iso- 
topic vector part of the structure, and the » -37 
vertex the most important one for the isotopic 
scalar part. Recent technical developments 
showing the feasibility of colliding beam experi- 
ments make it appealing to think of possible 
direct measurements of the photon-pion vertices 
through processes of the sort 





(1) 


We shall here discuss the reactions (1) in their 
lowest electromagnetic approximation. Examina- 
tin of the higher order terms in the electro- 
magnetic coupling constant may become nec- 
essary when detailed experiments are carried 
out. 

Let us consider the reaction (1) in the center - 
i-mass frame. In the lowest electromagnetic 
approximation only the *S, and *D, states for the 
initial positron-electron state will contribute to 
(i) and the final m-pion state must have parity 
minus, charge-conjugation quantum number 
ninus, total angular momentum one, and total 
isotopic spin one if » is even, zero if ” is odd. 
hparticular the process is forbidden in this 


+ - 
@ +e -n pions. 








a 


o= {a 
328 “(20 oa -5 





5g). a5 aol! 4a 


there a is the fine structure constant (another 
tctor @ is contained in |JI?), E is the electron 
t positron) energy, uw) is the energy of the _ 
thpion, and @ is the angle between the vector J 
ad the line of collision. The dependence on 6, 
ssin*@, is therefore a direct consequence of 
uge invariance independent of any knowledge 
tthe y-pion vertices. In the case of two final 
ons, with only one vector available, J must be 
Woportional to p(1) = -p(2) and the angular dis- 
tibution is thus completely determined of the 


yertices in the theory of the nucleon form factors. 


(2), 






PION FORM FACTORS FROM POSSIBLE HIGH-ENERGY ELECTRON-POSITRON EXPERIMENTS 


N. Cabibbo and R. Gatto 
Istituti di Fisica delle Universita di Roma e di Cagliari 
e Scuola di Perfezionamento in Fisica Nucleare dell’ Universita di Roma, Roma, Italia 
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order if the final pions are all neutrals. The 
S-matrix element for (1) is given by 


(1), (2) 


= 0p 4p) 4... +p 


“+p “-e 


(+) (-) 


-e ) 


Se 


(1) (2) gt) 


<% (ole “yy, ule ))<p 17 (0) 10), 


(2) 


where p(t), p(2).. -p™) are the final pion mo- 
mentum four mens e(+) and e-) the positron 
and electron momenta, respectively, K = e(+) +e(-) 
v(e*) and u(e~) are the Dirac spinors, and the 
matrix element of the electric current operator 
j (x) is taken between the vacuum state and the 
final state of m outgoing pions. 

We define 


wep) 


3n /2 


1 (0, p?), 


(1) (2) 40) 


=(27) ° <p Ij,,(0) 10). 


(3) 


It follows from gauge invariance that, in the 
center -of-mass system for (1), Jy has only the 
space component J. From parity conservation 
J is a polar vector formed from the n vectors 
p(1), 5(2), ... 5) if n is even, an axial vector 
formed from these vectors if n is odd. Inserting 
(3) into (2) one finds for the total cross section 


wo” -onyo@ 45? ....5)iF12sin26, (4) 





form sin*@. In the case of three final pions, J 
must be proportional to the only available axial 
vector p(1) xp(2) = -p(1) x5(3), ete, and therefore 
the distribution of the angle 6 between the normal 
to the plane in which the three pions are pro- 
duced and the line of collision is uniquely given 
by sin’6. 

The final state of two pions (necessarily one 
positive, one negative) produced according to (1) 
is necessarily a p state of relative orbital angu- 


313 

































VoLtuME 4, NUMBER 6 


PHYSICAL REVIEW LETTERS 





Marcu 15, 19% 





lar momentum. The matrix element (3) for two 
pions is expressed in terms of the pion factor 
M(K*) by’ 

J (D™, p™) =e(4w Mw) -M(K?)( pp -p™) . (5) 


The differential cross section, always in the 





center-of-mass frame, is given by 
do= a | M(-4E”) |?p°sin*@d(cosé), (6) 


where p= |p(1); = |p(2)|. Frazer and Fulco‘ pro- 
pose a resonant form of M(K*) as a simplest 
explanation of the isotopic vector part of the 
nucleon structure. According to (6) one finds 
for E =230 Mev, for which M(-4E*) comes close 
to its maximum, a total cross section for 
e++e-~n++n7 of 4.6x10™ cm?, which is =17 
times bigger than the value one would get from 
perturbation theory (M=1). The cross section 
then falls down very rapidly at higher energies. 
The final state of three pions (necessarily one 
positive, one negative, one neutral) produced 
according to (1) is a superposition of states with 
l=L=1,l=L=3, 1=L=5, etc., where / is the 
relative 7*z~ angular momentum and L the angu- 
lar momentum of the 7° relative to 7*n~. The 
matrix element (3) for three pions can be ex- 
pressed in terms of a form factor H depending 
on three independent scalars which can be chosen 
as E, w, (the energy of the positive pion), and 
w_ (the energy of the negative pion),? 
1 6), g, 9) 
-1/2,,* vpot, (+), (-), (0) 
H€ Py p 's -. 


. (7) 


= ~i(8w _%o) 


The vector J is then given by 


5) af ) 3) 


I@ ) 


=2(80 0 wp) /*H "EG xB). (8) 





With our choice of the independent kinematical] 
parameters the differential cross section can lk 
written as 


do i a 6! xp ~(-),2 
dw 44%. (27)? 64E? 





) sta” @d(cosé). (9 


Measurement of this cross section will thus in- 
form us directly about |H |? which plays a rele- 
vant role in the nucleon structure problem. At 
present there is no knowledge at all about this 
form factor. Even a rough knowledge of its 
average magnitude would be very important for 
our understanding of the nucleon structure. 
Bosco and De Alfaro have assumed a constant 
approximation for H, H =(Ae)/u°, and they claim 
that double pion photoproduction puts an upper 
limit =6 for A.° Ata given E the maximum of 
do is obtained for those configurations of the 
three final pions for which the area of the tri- 
angle formed with their momenta is maximum. 
For such configurations and for E =500 Mev, 
we find, using a constant H and the above value 
for A, d°o/dw,dw_=1.3 x10~* cm? Mev™. 

We are indebted to Professor W. K. H. Panof- 
sky for a stimulating seminar on the possibilities 
of colliding-beam experiments. 
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EFFECTS OF PION-PION INTERACTION IN ELECTROMAGNETIC PROCESSES 


L. M. Brown* 
Istituto di Fisica dell’ Universita, Roma, Italia 


and 


F. Calogero 
Istituto di Fisica dell’ Universita, Roma, Italia, 
and Istituto Nazionale di Fisica Nucleare, Sezione di Roma, Italia 
(Received February 5, 1960) 


Resonances in pion-pion scattering have been 
proposed by some authors.’ Recently Frazer 
and Fulco,” following the approach of Chew and 
Mandelstam, and assuming a resonance in the 
J=1, [=1 state of two pions, derived an elec- 
tromagnetic form factor for the pion which has 
a resonance -like behavior as a function of the 
invariant momentum transfer. They show that 
the position and width of the resonance can be 
adjusted to give a good fit to the isotopic vector 
part of the nucleon electromagnetic structure. 

We have calculated the effect of this pion form 
factor on the photon propagator, and have ob- 
tained corrections of order e” to several purely 
electromagnetic processes. In spite of the large 
mass of the intermediate pion pair, these cor- 
rections turn out to be perhaps detectable in 
electron scattering experiments at center -of- 
mass energies smaller than those expected from 
machines now under construction.* Such experi- 
ments would have the obvious advantage of giving 
information on the pion form factor independently 
of the complications of nuclear structure. 

Our method of calculation parallels that used 
by Kallén in his discussion of charge renormali- 
zation. Our modifications consist of using in- 
termediate bosons instead of fermions and of 
inserting the pion form factor. 

We will outline here only the scheme of cal- 
culation, reserving the details and further re- 
sults for publication elsewhere. We write first, 
following KAllén, the renormalized photon 
propagator in the momentum representation: 


Dy (P)= (2 -ie) "40, +(6,., Db, /P*) 
x[R(p?) +4 p?)}, (1) 
with 
R( p”) = 11(0) - Ti( p?), 
I( p?) = - 1 1( p?). 
The imaginary part J( p”), which involves the 


(2a) 
(2b) 






product of two current operators, is expressed 
as a sum over physical intermediate states 
through‘ 


II p*) = -(3p”)* (O1j le) te Ij, 10). (3) 
p\2)=p 


Restricting the sum to intermediate states con- 
sisting of two charged physical pions, we mul- 
tiply the usual first-order matrix element by 
the pion form factor® F,(p*) and carry out the 
sum. We then construct the real part 


p?Tl(-a) 
a( p? +a) 
where P designates the principal value. 

As the form factor of FF has a rather com- 


plicated form for analytical work, we use a 
simple analytic fit: 


IF (py? =1F (0%) 1?, 0>p?>p,? 


R(p*)= Pf da, (4) 


=1, py >? (5) 


where 
IF '(b*) |? =A[( p?+C)? +B?]", 
’ 2- , 2\;2- 
lf (0)! IF (Py )l? =1. (6) 
The numerical results which we give below for 
electron scattering are based on the choice 
A =66.25, B=1.5, C=8, (7) 


which corresponds to the FF form factor having 
a peak at p?=-8, using always units of the pion 
mass. (The other FF form factors give slightly 
smaller corrections.) 

The final result is 


II( p?) = (a /6n)z**0(-p? - 4) IF (p*) FP, (8) 
R(p”) =(a/6n){[1 - IF (bp?) I? VS (2/25) -f(z) 


+p? IF" p”)|?[-y + 6(C +p?)/B]}, 
(9) 


315 

































VotuME 4, NUMBER 6 PHYSICAL REVIEW LETTERS MarcH 15, 19%) 




















where ee Ss ee 
0.5 1 
a@ =1/137, x | 500 Mev (a) | 
z=1+4/p?, 2z,=1+4/p,?, r 250 Mev ] 
O(x>0)=1, O6(x<0)=0, ¥ J 
, a nd 7, «@ «@ 2 2 6 2 2 2 san 
F(z) = 5 +22 +2~”*in za] ° (10) 1.0 

% (b) J 
For real positive z, the absolute value of the J 

argument of the logarithm must be taken in (10). 220 Mev 


For negative z, one must replace z™In[(z¥ -1)/ 
(z* +1)] by -2yarc coty, with y =(-z)”*. For 
complex values of z, f(z) is defined in a plane 
cut in the negative real axis. The constants y 
and 6 are defined by® 

y +46 = 4(1 -2’)z,**f (z’/z,) (11) 
with 

z’=1-4/(C +iB). 


We have applied this result to evaluate the 
correction to Mgller and Bhabha scattering of 











electrons due to pion vacuum polarization. The - r 1 
- a L 1 1 1 aa ee 4 i i Doel calle 

fractional correction k(E,cos@), where E is the 0 82 86 ~# 2 0 2 4 7. 

center-of-mass energy of each electron and @ cos § 

the c.m. angle of scattering, takes the form in FIG. 1. Percentage correction as a function of c.m. 

the extreme relativistic case: angle for fixed c.m. energy of each electron: (a) 


Mgller: k=(1+X)R(p,”) +(1 -X)R(p,”), (12) Mller scattering and (b) Bhabha scattering. 


Bhabha: k=(1+Y)R(p,”)+(1-Y)R(p,”), (13) 





with 
X =4c0s6(3 +cos*6), (14) 
Y =(1+cos@)(3 - cos@)(3 +cos*6)", (15) 
and with 


p,” =p,” = 2E?(1 - cos6), 
p,” = 2E?(1 +cos8@), 
p= -4E?. (16) 





In Figs. 1 and 2 we have plotted some examples 
of these results. 

The largest corrections are obtained for large- 
angle Bhabha scattering, which is dominated by 
the timelike momentum transfer p,’. For this 
reason the large-angle plot in Fig. 2 reproduces 
essentially the correction to the photon propaga - 





a a cr | 


. * v7 
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tor for timelike p?. From an experimental curve 0 100 200 300 oe ., © 
showing this behavior could be deduced unam- 

biguously the position and width of the resonant FIG. 2. Percentage correction to Bhabha scatter- 
part of the pion form factor, free of the uncer- ing as a function of c.m. energy of each electron for 


tain interpretations of nuclear structure. fixed c.m. angle. 
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Of course annihilation processes leading to 
muon pairs or to two or more pions will become 
of interest at these energies. For the muon case 
we get a small correction similar to that for 
large-angle Bhabha scattering. For the very 
important case of annihilation into pion pairs 
one is studying directly the pion form factor. 

For example, at 90° in the c.m. system one gets 


satsabnay (fe) IF, 48" (17) 


do(Bhabha) '*\E? 
where p and E are the momentum and energy of 
each pion. This is certainly, therefore, the 
best means of investigating the pion form factor. 
Ifa resonant pion form factor is found in this 
way, then the small effects we have calculated 
become of relevance in the interpretation of 
scattering experiments in terms of a possible 
failure of quantum electrodynamics at small 
distances.” 

However, while the positron-electron experi- 
ment appear to be the most promising for eluci- 
dating the pion form factor, electron-electron 
experiments will probably be done first and the 
above corrections for Mgller scattering may be 
of interest. We are investigating other electro- 


magnetic processes in which pion-pion inter- 
action plays a role. 

It is a pleasure to thank Professor M. Cini for 
encouragement and advice. 
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INVARIANT COMMUTATORS FOR THE QUANTIZED GRAVITATIONAL FIELD* 


Bryce S. DeWitt 
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(Received February 8, 1960) 


The problem of developing a formalism within 
which a union of the general theory of relativity 
with the quantum theory might successfully be 
brought about has been attacked by many authors.’ 
The difficulties encountered in this work have 
seemed to be twofold: (1) the nonlinearity of the 
general theory, and (2) its coordinate invariance, 
which leads to constraints on the Cauchy data for 
the dynamical equations. The real difficulty may 
lie elsewhere, however. Approaches to the pro- 
blem have almost without exception been made 
via a Hamiltonian or quasi- Hamiltonian canonical 
formalism, with an attendant loss of manifest 
covariance. It is the purpose of this Letter to 
raise the question of the suitability of canonical 
procedures and to suggest, by obtaining, without 
their aid, an explicit covariant expression for 
the commutators of the theory, that it may be 
Possible to avoid them entirely. 





We begin by pointing out that, apart from pro- 
blems of factor ordering, the nonlinearity of the 
dynamical equations offers no complications for 
the commutators, contrary to a widespread im- 
pression. As Peierls? has shown, the commu- 
tator of two dynamical variables is determined 
by the variation in one due to an infinitesimal 
change in the action proportional to the other. 
But infinitesimal variations are propagated by 
means of linear equations, and standard tech- 
niques are available for handling the associated 
propagation functions. 

Secondly, we note that in the general theory of 
relativity, well-defined commutators can be ob- 
tained only between absolute invariants. The in- 
finitesimal coordinate transformation law, 


, Gi infinitesimal) (1) 


wv "psy vey 
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leads to the characterization 


(64/dg, |), =9, (6B/dg |), ,=0 (2) 


of any two absolute invariants A and B, which 
are constructed solely out of the metric field 

Su yp’ Equations (2) may hold either as identities 
or in virtue of the field equations. 


Consider now the following change in the action: 


S-S+€eB, (3) 


where € is an infinitesimal constant. Such an 
alteration will produce a change in the invariant 
A, the precise nature of which depends upon the 
boundary conditions selected. For example, we 
may adopt advanced or retarded boundary con- 
ditions, in which the dynamical states of the 
system before and after the alteration (3) are 
taken to coincide in the remote future or remote 
past, respectively. Denoting the corresponding 
changes in A by 5*A and 6"A, respectively, we 
may express the commutator of A and B in the 
form 


5°A - 5°A=-ie[A, B]. (4) 


(We choose units in which h=c=167G=1.) That 
Eq. (4) follows, in the case of normal canonical 
systems, from the usual canonical definition of 
the commutator (Poisson bracket) has been 
proved by Peierls.” In the present case we shall 
adopt Eq. (4) itself as the definition of the com- 
mutator, the utility of this procedure resting on 
the fact that an independent evaluation of 6*A-67A 
is possible. 

Ignoring, for the present, problems of factor 
ordering, we begin by writing 


ee + 
ofA =[(0A/0g_,)8 gtx. (5) 


Since g,,, is not an invariant its variations 
guy pond determined only up to a coordinate 
transformation (1). In virtue of (2), however, 

the variations 6*A remain well defined. The dif- 
ferential equations satisfied by the b*g y are 
obtained by varying the original field equations, 


0=0S/og, =G=g(RY-agh*R). (6) 
(Here g is the negative of the determinant of 
&uy and R#” is the Ricci tensor.) Under the al- 
teration (3) the new field equations lead to 

+ pv 
6 =- ° 
G €5B/ By (7) 
It is not hard to show that when the field equa- 
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tions (6) are satisfied, the left side of (7) be- | 
comes® | 


<_< es HO VT aghY OT, pr gt 
6 =ig” 
G 28 "(gg g lg” (6 OT*pr 
+ + + 
6 -6 3 
: 8p x-oT Bop .2-° 85-00) (8) 


One may then easily verify that the variation (1), 
if substituted into (8), gives zero when (6) is 
satisfied. The general solutions of Eqs. (7) are 
therefore obtained by adding (1) to particular 
solutions determined by appropriate boundary 
and supplementary conditions. 

As supplementary conditions it is convenient 
to choose 


po J a 1 af. OT 
(g 28 )o~ a. (9) 


Equations (7) then take the form 


oO VT V OT A,# 
ging ght gh gg PX5 + Bnae 


MovT.+ " 
2g"*R ae 2e5B/dg (10) 


of which the solutions are 


6B 
og =€ G (x,z).——d@*z, (11) 
pv Lvap 08, 8 
where G* are the advanced and retarded 


Green’s functions for Eq. (10), satisfying 


OT + CG Tz 


paper’ *,, v a oTap “(8 


ad pa® vB 


o g « - (4) “14 
Fue va ,Fap® V4(~)5(x,z)g“*(z). (12) 


Green’s functions of this type have been studied 
by Brehme and the author,’ who have called 
them bi-tensors on account of their transforma- 
tion properties at two different space-time 
points. The notation of reference 7 is followed 
here. Indices from the first part of the Greek 
alphabet refer to the point z and those from the 
latter part to the point x. The quantity Bua is 
the bi-vector which effects parallel displace- 
ments along the geodesic between x andz. The 
4-dimensional delta function is regarded as a 
scalar density of weight $ at both x andz. The 
general nature of the propagation characterized 
by the above Green’s functions has been de- 
scribed by Hadamard.® Explicit forms for such 
Green’s functions are given in reference 7. 
Typically, they are singular on the light cone 
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(delta function of the square of the geodetic in- 
terval), nonvanishing but finite inside, and 
vanishing outside. 

It is now important to check that the solutions 
(11) satisfy the supplementary conditions (9) 
which were invoked in order to get.them in the 
first place. By using the techniques of reference 
1, it may be shown that 


po vT , wv o7) + 


(g & -28 oTaBev 


GH _g*# (13) 


a-p B-a’ 
vhere G*4, are the advanced and retarded bi- 
vector Green’s functions satisfying 


vo +p 
8 G a-vo 
== Bh gM x)0%(x,2)¢7 2). (14) 
From (11) and (13) it then follows that 
Mo vT 4 pV OT. 
(gng’? - agg 8 _ 
(15) 


+p 6B 
--26f6 (x,z) ——d@tz, 
ap Bg 


which, after an integration by parts and use of 
Eq. (2), leads back to (9), showing the complete 
self-consistency of these supplementary condi- 
tions. 

Finally, therefore, upon substituting (11) into 
(5) and making use of (4), we obtain, for the 
commutator, 
is, B]=if atx d*z ——G 


(x,2) 5 —— 00 


Su vep 8g 


where 


+ - 

ee poe” © vas (17) 
To the extent that it is possible to ignore the 
honcommutativity of factors on the right-hand 
side of (16), this commutator can be shown to 
satisfy all the usual identities. Antisymmetry 
follows immediately from the reciprocity rela- 
tion satisfied by the propagation function’: 


Gi, vag 


G yag'*2)®- G4, 2 >*)- (18) 
The identity, 
[A, BC]=[A, B]C+BlA,C], (19) 


is straightforward. The proof of the Jacobi 
identity, on the other hand, is quite lengthy since 
it involves variational derivatives of the propa- 
gation function. These derivatives can be com- 
puted by varying the defining equations (12); the 
proof then carries through, after a considerable 
amount of labor. 

The first step in attempting to go beyond these 
results by taking into account the noncommuta- 
tivity of factors consists of establishing an algo- 
rithm for the evaluation of the right-hand side of 
Eq. (16). One must require that the propagation 
function G,,, 9g is first to be inserted as a re- 
placement for 58a in all the places in which it 
occurs in the variation 5B, and that the result- 
ing “product” is then inserted as a replacement 
for 5g,,,, in the variation 6A, or, alternatively, 
that the process of insertion is first performed 
in 5A and then in 6B. That the two procedures 
are equivalent is evident from the familiar pro- 
perties of commutator brackets when one takes 
note of the fact that in computing the commutator 
of two absolute invariants one may work directly 
with the factors (the metric tensor and its deri- 
vatives) out of which they are built, proceeding 
as if the metric tensor had the simple commuta- 
tion rule, 


ey Fag! *E vag" 
It is clear from these considerations that the 
question of the consistency of a quantum theory 
of gravitation falls squarely on the propagation 
function, which, contrary to the case of linear 
field theories, is here an operator. The problem 
of a consistent definition of this “operator-pro- 
pagator” is under active consideration at the 
present time. 


(20) 





"Work supported by the U. S. Air Force Office of 
Scientific Research. 

‘Pp. A. M. Dirac, Proc. Roy. Soc. (London) A246, 
333 (1958); Phys. Rev. 114, 924 (1959). S. Deser and 
R. Arnowitt, Phys. Rev. 113, 745 (1959); also subse- 
quent papers to appear in collaboration with C. W. 
Misner. For references to earlier work see P. G. 
Bergmann, Suppl. Helv. Phys. Acta, IV, 79 (1956). 

*R. E. Peierls, Proc. Roy. Soc. (London) A214, 
143 (1952). 

3Covariant differentiation is here indicated by the 
dot followed by an index. The notation for variational 
differentiation appearing in Eq. (2) is the standard 
adopted for most work in quantum field theory. [See, 
for example, N. N. Bogoliubov and D. V. Shirkov, 
Introduction to the Theory of Quantized Fields (Inter- 
science Publishers, New York, 1959), pp. 202, 203.] 
Its meaning is explained by Eq. (5). 
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‘The concepts “remote past” and “remote future” 
are assumed to have an absolute invariant significance. 
For space-times in which this is not the case, altera- 
tions will be required in our results. 

5Here the covariant derivatives are taken with res- 
pect to the original metric. 

6G” has the coordinate transformation law, 


BY _/)AMY,o, oh ov _,v po 
6G =(G'e"), 8 Gene” Gh. 


—, 


7B. S. DeWitt and R. W. Brehme, Ann. Phys. (to 
be published). 

83. Hadamard, Lecture’s on Cauchy’s Problem in 
Linear Partial Differential Equations (Yale University 
Press, New Haven, 1923). 

®*The reciprocity relations for the laws of propaga- 
tion are proved in reference 8. The corresponding 
identities for the explicit propagation functions are 
proved in reference 7. 
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SEARCH FOR THE REACTION p* +e" —y+y. 
C. M. York, C. O. Kim, and W. Kernan [Phys. 
Rev. Letters 3, 288 (1959). 


Equation (4) should read: 
R =(n/2 4 2 P. 
P (1/ jaa, , pe/(na,) 


This changes the numerical result for a, in 
Eq. (5) to a, <1.1x10°~". The authors are in- 
debted to Mr. J. Good and Professor Y. Nambu 
for calling their attention to this error. 


PARTIAL-WAVE DISPERSION RELATIONS FOR 
MESON-NUCLEON SCATTERING. Reinhard 
Oehme [Phys. Rev. Letters 4, 246 (1960)]. 


The second line below Eq. (3) should read: “due 
to s>(m+u)*; (nN) from 0 to + in and...” 


In Eq. (8), N and M should be interchanged. The 
second line of reference 5 should read: “become 
complex for positive or negative values of 8. 
pe 

FERROMAGNETIC RELAXATION MECHANISM 
for Mz IN YTTRIUM IRON GARNET. M. Sparks 
and C. Kittel [Phys. Rev. Letters 4, 232 (1960)]. 


The line below Eq. (3) should read: 
where F(@) =1+(17/2)sin?@ - (35/4)sin‘é. 


In paragraph 6, T>>0.3°K should read T >>1.3°K. 
We would like to thank Dr. R. C. Le Craw for 
calling our attention to this error. In the sev- 
enth paragraph, the corrected values of 1/7, 
are 0.98 x10° Tsp sec™ and 0.8[8]x10° T,, sec”, 
respectively. 
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ABSTRACTS 


In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 











jore the appearance of the correspondin g Article, ref- 
erence should be made to “Physical Review (to be 
published )”” rather than to this Journal. 


QUANTUM MECHANICAL TRANSPORT THEORY. 
L INCOHERENT PROCESSES. Kenneth M. Wat- 
son, Department of Physics, University of Cali- 
fornia, Berkeley, California (Received January 
14, 1958). 


The transport of particles through a scattering 
medium is studied. A generalization of a tech- 
nique due to Placzek and Wick is used to handle 
sums over states of excitation of the medium. 
The collision processes which occur are classi- 
fiedas “inelastic,” “elastic, ” and “quasi-elastic” 
and correspond to different orderings of the 
Placzek-Wick series. The inelastic scatterings 
are described by an essentially classical trans- 
port equation and the elastic scatterings by as- 
signing a refractive index to the medium. The 
“quasi-elastic” scattering involves the excita- 
tion of low-lying states of the scattering system. 
The coherent interference of waves scattered 
from near-by scatterers is important in this 
case and depends upon the structure of the me- 
dium. In this paper the general theory is devel- 
oped in terms of a systematic sequence of ap- 
proximations, of which the first gives just the 
classical form of transport theory. The correc- 
tion terms then appear as quantum- mechanical 
corrections to the classical transport problem. 


THERMAL CONDUCTIVITY OF ISOTOPIC MIX- 
TURES OF SOLID HELIUM. Edward J. Walker 
and Henry A. Fairbank, Laboratory of Marine 
Physics, Yale University, New Haven, Connect- 
icut (Received December 2, 1959). 


The thermal conductivity of solid He* containing 
zero, 0.56, 1.38, and 2.8% He® has been meas- 
ured as a function of temperature from 1.1 to 
2.1°K. For pure He‘, the results are in satis- 










factory agreement with those of Webb, Wilkinson, 
and Wilks. The addition of 1.38% He*, at a sam- 


ple density p= 0.208 g cm~*, caused a decrease 





in thermal conductivity by a factor of 5 at 1.1°K. 
The additional thermal resistance caused by 
adding He* was substantially independent of tem- 
perature, rather than proportional to tempera- 
ture, as would be expected from isotropically 
distributed point scatterers. From this fact and 
from the magnitude of the isotope resistance it 
is inferred, using Klemens theory, that the He*® 
is arranged on lines in the solid. A possible 
mechanism is that the He® has been attracted to, 
and deposited along, dislocations in the He* 
matrix. 


VAPOR PRESSURES OF He® - He* MIXTURES. 
S. G. Sydoriak and T. R. Roberts, Los Alamos 
Scientific Laboratory, University of California, 
Los Alamos, New Mexico (Received December 
7, 1959). 


Vapor pressures, P,, of He* - He* mixtures 
ranging in liquid mole fraction, X, from 0.1 to 
0.9 have been measured between 0.6° and 2.4°K 
versus the vapor pressure, P,°, of liquid He’. 
Except for sharp breaks in the vicinity of the 
lambda and stratification temperatures of some 
of the mixtures, P,/?S is found to vary slowly 
and smoothly with X and with temperature. 
Contrary to much of the work of other authors, 
there are no breaks at the He* lambda tempera- 
ture. 

A comprehensive smoothed table of P,/P,° is 
derived. Using this as a reference it is possible, 
for the first time, to intercompare all of the 
previously existing data on P,. The data of 
some authors are in excellent agreement with 
our table but other data are in serious quantita- 
tive disagreement near 1.2°K and much previous 
data are in qualitative disagreement at the He* 
lambda temperature. 


STABILITY OF UNIFORM PLASMAS WITH 
RESPECT TO LONGITUDINAL OSCILLATIONS. 
Peter D. Noerdlinger, Normal Bridge Labora- 
tory, California Institute of Technology, Pasa- 
dena, California (Received December 2, 1959). 


It is possible to relate the dispersion formula 
for longitudinal oscillations in an infinite, uni- 
form, collision-free plasma with no magnetic 
field to the complex potential of a line charge 
distribution on the real axis of the phase velocity 
(u=w/k) plane. If the initial velocity distribu- 
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tion averaged over directions orthogonal to the 
direction of propagation is f,(v), the plasma is 
stable if and only if 


viuy-p [ L’wdae 
v 


-u 
-& 
is negative at the minima of /,(v) on the real 
axis, with unimportant exceptions. In particular 
it is shown that single -peaked distributions are 
stable, while those with very sharp (e.g., non- 
differentiable) minima or with a zero of f, be- 
tween two peaks are not. The charge analogy 
yields information on the wavelengths for which 
oscillations can grow and on rates of growth. 
Examples are given, including the case of two 
identical interpenetrating hot plasmas. A lim- 
ited generalization to transverse oscillations is 
given. 


NEGATIVE CURRENT-VOLTAGE CHARACTER- 
ISTICS IN HYDROGEN AT HIGH PRESSURE 
USING PLANE PARALLEL ELECTRODES. D. J. 
De Bitetto* and L. H. Fisher, Department of 
Physics, New York University, University 
Heights, New York, New York, and A. L. Ward, 
Diamond Ordnance Fuze Laboratories, Wash- 
ington, D. C. (Received December 17, 1959). 


In conjunction with measurements of current- 
voltage characteristics in hydrogen, a few char- 
acteristics have been obtained which include a 
region with negative slope. The latter charac- 
teristics were obtained with plane parallel 
electrodes at an electrode separation of 2 cm at 
a pressure of 400 mm Hg and with three values 
of externally initiated cathode current. The 
initial currents ranged from about 107" to 10~° 
amp, and the amplified currents reached values 
as high as 10-* amp. The characteristics cor- 
responding to the larger initial currents become 
negative at large currents (~10-° amp). The 
voltage at which a characteristic becomes neg- 
ative, i.e., the maximum attainable voltage 
across the electrodes, decreases slightly with 
increasing initial current. The circuit included 
a series resistor of 20 megohms. 

These characteristics can be explained quan- 
titatively on the basis of the first and second 
Townsend coefficients (previously measured 
with the same apparatus) acting in conjunction 
with space charge, if a not unreasonable dis- 
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charge area is assumed. These calculations 
were carried out on an IBM 704 computer. 


*Now at Philips Laboratories, Irvington-on-Hudsop, 
New York. 


RANDOM-WALK INTERPRETATION AND GEN. 
ERALIZATION OF LINEAR BOLTZMANN EQus. 
TIONS, PARTICULARLY FOR NEUTRON 
TRANSPORT. E. Guth and E. Inonu, Oak Ridge 
National Laboratory, Oak Ridge, Tennessee 
(Received December 21, 1959). 


The connection between linear recurrence 
relations which define generalized random walks 
and the related linear Boltzmann equations is 
clarified. The probability distributions f,,(s) for 
“the state s” reached by a “random walker” 
ae n steps satisfies the recurrence relation 

(4= Sf,(s)P(s, t)ds, where the non-negative 
yg? t) is the probability for a transition from 
s tot. The Boltzmann distribution is given by 

f(s)=U,,-ofn(s). In general, F,{S) contains more 
information than f(s). Moreover, F,{S) is the 
nth term in the iteration series solution of the 
Boltzmann equation and therefore can also be 
obtained from the solution of an associated 
Boltzmann equation which contains an additional 
parameter. As an example, the well-known 
integral Boltzmann equation for neutron trans- 
port in a nonmultiplying infinite medium is 
derived from a P(s,¢) which involves a transi- 
tion in a seven-dimensional phase-time space. 
Brownian motion and Rayleigh’s problem (re- 
lated to neutron thermalization) may be treated 
similarly. 


ZEEMAN EFFECT OF IMPURITY LEVELS IN 
SILICON. Solomon Zwerdling, Kenneth J. Button, 
and Benjamin Lax, Lincoln Laboratory, Massa- 
chusetts Institute of Technology, Lexington, 
Massachusetts (Received December 28, 1959). 


Completely resolved Zeeman spectra for the 
bismuth donor in silicon including optical transi- 
tions from the 1s donor ground state to the ex- 
cited states 2p,, 2p1, 3p), 4P5, 3p4, Spo, 4h 
and 5p, are presented. The transitions were 
observed at liquid helium temperature, using 
linearly polarized radiation alternately parallel 
and perpendicular to the magnetic field, and 
field intensities up to 38.9 kilogauss oriented 
along each of the three principal crystallographic 
axes. Both linear splitting of the p, states and 
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, quadratic dependence on field were observed. 
The use of impurity Zeeman spectra is demon- 
strated for evaluating effective mass parameters, 
determining the nature of energy bands and find- 
ing and identifying impurity excited states. The 
transverse effective mass for the electron in 
silicon was found to be (0.186 + 0.006)m, in 
agreement with recent cyclotron resonance re- 
suits. From Zeeman splitting, electron effec - 
tive masses up to 0.5m, can be measured to 
within +4% at infrared frequencies in a field 

of 40 kilogauss. The behavior in a magnetic 

field of the first two donor excited states could 
be explained by treating the magnetic terms of 
the Hamiltonian as a perturbation to first order. 
Interactions among the higher closely-spaced 
Zeeman levels were observed above 20 kilo- 
gauss and were evaluated with a second-order 
treatment. The Zeeman structure for the 
aluminum acceptor reflected the complexity of 
the valence bands and the acceptor ground state 
and was in qualitative agreement with the theo- 
retical results of Kohn and Schechter. Transi- 
tions were observed to eight excited states 
converging to the series limit. Evidence is 

given for the degeneracy of each state. 


CROSS RELAXATION STUDIES IN DIAMOND. 

P. P. Sorokin, G. J. Lasher, and I. L. Gelles, 
Research Laboratory, International Business 
Machines Corporation, Poughkeepsie, New York 
(Received December 18, 1959). 


A microwave double resonance experiment 
performed on the paramagnetic nitrogen centers 
in diamond shows that in this system cross re- 
laxation occurs via a four-spin-flip mechanism 
which exactly conserves Zeeman energy. In 
this process, which was first postulated by 
Bloembergen and co-workers in their paper on 
cross relaxation, two spins of the center line 
make a downward transition while a spin belong- 
ing to each satellite makes an upward transition. 
Simple rate considerations for this process 
indicate that if a saturating microwave field is 
suddenly applied to one of the three lines of the 
nitrogen spectrum, a weak probing microwave 
signal at either of the two lines should register 
adefinite change in absorption in a time T,,. 
Specifically, if T,, is much less than other re- 
laxation times of the system, then setting the 
pimp upon the center line should force the ab- 
sorption at either satellite to drop to zero. Set- 





ting the pump at the position of one of the sat- 
ellites, on the other hand, should reduce the 
center line absorption to 3/5 its thermal equi- 
librium value but should increase the absorption 
measured at the other satellite by the factor 6/5. 
This behavior was precisely observed at 1.6°K, 
using a bi-modal cavity. 

By resolving the rate at which a satellite de- 
cays to zero when the pump is set on the center 
line, T,, is measured for all five satellites in 
the three principal orientations: H,||[ 100], 
H,|\[110], and H,||[111]. The measured anisot- 
ropy is discussed. 

It is shown that the four-spin-flip transition 
may be used in special cases to establish cw 
maser operation by inverting the population of 
one of the satellite lines. Steady-state inversion 
of one of the nitrogen satellites is incidentally 
observed in a number of diamonds. 


INFLUENCE OF COLD WORK ON THE RESIST- 
IVITY OF DILUTE COPPER ALLOYS. P. G. 
Klemens,* R. J. Tainsh, and G. K. White, Divi- 
sion of Physics, Commonwealth Scientific and 
Industrial Research Organization, Sydney, Aus- 
tralia (Received November 30, 1959). 


Measurements are reported of the residual 
electrical resistivity of a copper + 1 atomic per- 
cent antimony alloy which had been cast and 
rolled. A large increase was observed on anneal- 
ing, as noted by Kropschot, Garber, and Blatt. 
However, metallographic examination and the 
smallness of the change observed on subsequent 
rerolling suggest that their proposed mechanism 
of “migration” during cold work plays a much 
smaller part than lack of homogeneity in the 
cast alloy. 


“Present address: Westinghouse Research Labora- 
tories, Pittsburgh, Pennsylvania. 


TRANSVERSE COLLECTIVE EXCITATIONS IN 
SUPERCONDUCTORS AND ELECTROMAGNETIC 
ABSORPTION. T. Tsuneto, Department of 
Physics, University of Illinois, Urbana, Illinois 
(Received December 24, 1959). 


With use of the generalized random phase 
approximation an attempt is made to estimate 
the absorption of photons with energy less than 
the energy gap due to transverse collective ex- 
citations. The ratio of the surface resistance 
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due to transverse collective excitations to that 
of normal metals in the extreme anomalous 


limit, calculated within the weak-coupling theory, 


turns out to be too small to explain the observed 
data for superconducting lead and mercury. The 
interpretation of the collective excitations as 
bound pair states is briefly discussed. 


MAGNETIC SCATTERING OF NEUTRONS BY 
EXCHANGE-COUPLED LATTICES. A. W. 
Sdenz, U. S. Naval Research Laboratory, Wash- 
ington, D. C. (Received September 21, 1959). 


The magnetic scattering of slow neutrons of 
arbitrary initial polarization by an extensive 
class of magnetically-coupled lattices is treated 
by a time-dependent operator approach. This 
magnetic scattering is carefully divided into 
purely magnetic and magnetovibrational scatter- 
ing, the types thereof involving, respectively, 
only zero-phonon processes and solely nonzero- 
phonon processes, and general formulas for 
these two types of scattering are obtained. 
These formulas are applied in temperature re- 
gions which are sufficiently large (I) or suffi- 
ciently small (II) compared with the temperature 
above which paramagnetism obtains. In region 
I, for the purely magnetic scattering and under 
certain invariance requirements on the above 
magnetic coupling, we analyze the energy spec- 
trum of outgoing initially unpolarized neutrons 
of sufficiently high incident energy by a moment 
method. We thus obtain general formulas for 
the energy-integrated effective differential cross 
section defined in this paper and for the mo- 
ments of energy transfer defined therein. These 
formulas involve certain spin averages, ex- 
plicit equations being given for a wide range 
of these averages for exchange- coupled lattices. 
These results are illustrated numerically and 
compared with experiment for the case of poly- 
crystalline MnF,. In region II, we discuss cer- 
tain broad features of the purely magnetic one- 
magnon scattering of arbitrarily polarized neu- 
trons by exchange-coupled lattices of the class 
alluded to previously and by more complex ones, 
studying this scattering in detail for ferromag- 
nets and certain antiferromagnets. A new spin- 
wave effect is pointed out for polarized neutrons 
incident on ferromagnets. Brief treatments of 
the magnetovibrational scattering in regions I 
and II are also given. 
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THEORY OF NUCLEAR QUADRUPOLE INTER. 

ACTION IN BERYLLIUM METAL. M. Pomerany 
and T. P. Das,* Department of Physics, Univer. 
sity of California, Berkeley, California (Receiyy 
December 29, 1959). 


The theory of the origin of the field gradient a 
nuclei in metals has been analyzed. The contri- 
butions of the ion cores and conduction electrons 
have been separately considered. In the case of 
beryllium metal, using orthogonalized plane 
wave functions, the conduction electrons are 
shown to enhance, by about eight percent, the 
field gradient due to the ion cores. Combining 
the results of our calculations with Knight’s ex- 
perimental value of 48 kc/sec for the Be® coup- 
ling constant e?¢Q/h, a value of Q =0.029x10™™ 
cm? is obtained. The dependence of the potential 
for the conduction electrons on the model chosen 
is analyzed in some detail. The various uncer- 
tainties in our field-gradient calculation and the 
theoretical value of the Knight shift in beryllium 
metal are discussed. 


“On leave of absence from Saha Institute of Nuclear 
Physics, Calcutta, India. Present address: Depart- 
ment of Chemistry, Columbia University, New York, 
New York. 


MEASUREMENT OF LIFETIME IN Ge FROM 
NOISE. Susumu Okazaki and Hiroshi Oki, 
Department of Physics, Faculty of Science, 
Okayama University, Okayama, Japan (Received 
December 14, 1959). 


The lifetime of the minority carrier can be ob- 
tained from a noise measurement whose method 
consists of liberating hole-electron pairs by 
light. A convenient experimental arrangement 
is suggested and a typical plot of the experi - 
mental results is shown. The curves of the re- 
sults are compared with the calculated curves. 


THE ELECTRONIC STRUCTURE OF TIN IN- 
VESTIGATED BY ULTRASONIC ATTENUATION. 
Tore Olsen, * Department of Physics, Brown 
University, Providence, Rhode Island, and 
Fysisk Institutt, Universitetet i Oslo, Blindern, 
Norway (Received December 1, 1959). 


The magnetic field dependence of the ultrasonic 
attenuation has been measured in very pure tin 
single crystals. Oscillations were found that 
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can be explained as a result of resonant condi- 
tions between the electron orbit diameter and 
the periodic field set up by the sound wave. 
These oscillations yield information about the 
fermi momentum, and the general features of 
apossible electron distribution in tin are sug- 


dient at gested. 
a "present address: Fysisk Institutt, Universitetet 
ase of i{Oslo, Blindern, Norway. 
ine 
are 
the | SATURATION AND RECOVERY OF A PORTION 
ining | OF THE ELECTRON SPIN RESONANCE OF F 
’s ex. | CENTERS IN KCl AT 4°K. Gordon A. Noble, 
coup- | Research Department, Zenith Radio Corporation, 
10°" | Chicago, Illinois (Received December 14, 1959). 
tential The inhomogeneously broadened magnetic rese- 
chosen nance of F centers in additively colored potassium 
ncer~ | chloride was photographed using a field modula- 
- the tion which swept rapidly across the resonance 
yilium as the absorption was displayed on the oscillo- 
scope. The time between sweeps was long and 
iclear | adjustable. A concentration dependance of the 
spart- | spin lattice relaxation time was observed at 
York, § 4.2°K. In crystals with a concentration near 
1x10" cm™$, a portion of the resonance could 
be saturated without affecting the rest of the 
broad spectrum. The recovery time of this 
»M “hole burned in the resonance” is about a quarter 
ofa minute. The results are compared with 
present theories on spin diffusion. 
eived 
e ob- § EFFECTS OF HYDROSTATIC PRESSURE ON 
sthod § THE PIEZORESISTANCE OF SEMICONDUCTORS: 
y i-InSb, p-Ge, p-InSb, AND n-GaSb. R. W. 
ent | Keyes and M. Pollak,* Westinghouse Research 
: laboratories, Pittsburgh, Pennsylvania (Re- 
a ceived December 10, 1959). 
es. A method for measuring the piezoresistance 
ofa sample under high hydrostatic pressure by 
comparison with the piezoresistance of intrinsic 
iiSb is described. The method is tested by a 
» measurement of the piezoresistance of p-type 
'ION. } germanium and p-type InSb up to 13000 kg/cm?. 
; The piezoresistance of these materials is found 
fo be independent of pressure, in good agree- 
rn, _} tent with predictions based on other experi- 
ments. Measurement of the piezoresistance of 
sonic § "-GaSb as a function of pressure up to 12 000 kg/ 












tm’ confirms the correctness of the model of 
the conduction band which has been proposed by 


Sagar. Values for certain parameters of the 
conduction band of GaSb are deduced. 


*Now at Bell Telephone Laboratories, Murray Hill, 
New Jersey. Permanent address: Israel Institute of 
Technology, Haifa, Israel. 


SURFACE IMPEDANCE OF SUPERCONDUCTORS. 
Piotr B. Miller, Department of Physics, Univer- 
sity of Illinois, Urbana, Illinois (Received 
December 21, 1959). 


A detailed calculation of the surface impedance 
of superconductors is given based on the general 
theory of the anomalous skin effect in normal 
and superconducting metals given by Mattis and 
Bardeen. It is found that there are large cor- 
rections to the extreme anomalous limit value 
of the superconducting to normal surface re- 
sistance ratio; corrections to the surface react- 
ance ratio are much smaller. The theory is 
compared with recent experiments on the sur- 
face impedance of aluminum and of tin. It is 
found that the theory gives satisfactory agree- 
ment with experimental data on the surface 
impedance, both in absolute value and in its 
temperature and frequency dependence over a 
wide range of temperatures and frequencies. 


HEAT CAPACITY OF SODIUM AND POTASSIUM 
AT TEMPERATURES BELOW 1°K. William H. 
Lien and Norman E. Phillips, Department of 
Chemistry and Lawrence Radiation Laboratory, 
University of California, Berkeley, California 
(Received December 17, 1959). 


The heat capacities of sodium and potassium 
have been measured from about 0.15°K to just 
above 1°K. No anomaly was observed in either 
case: the heat capacity could be represented by 
the sum of a term linear in temperature and a 
cubic term. Values of the Debye temperatures 
and effective masses of the conduction electrons 
are given. 


TRANSPORT NUMBER IN SOLID CESIUM BRO- 
MIDE. Neal Laurance, Department of Physics, 
University of Illinois, Urbana, Illinois (Re- 
ceived December 9, 1959). 


The transport number of cesium bromide sin- 
gle crystals has been measured over the tem- 
perature range 350° to 450°C by the method of 
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Tubandt. No temperature dependence is ob- 
served, and the average value of the cation 
transport number is 0.49+0.05. A preferential 
growth phenomenon is described, and its pos- 
sible interference with the measurement is 
discussed. It is estimated that the true cation 
transport number may be as low as 0.3 because 
of this phenomenon. 


RECTIFICATION WITHOUT INJECTION AT 
METAL-TO-SEMICONDUCTOR CONTACTS. 

N. J. Harrick, Philips Laboratories, Irvington- 
on-Hudson, New York (Received December 7, 
1959). 


A recent publication shows that extraction in 
the semiconductor bulk may occur for either 
direction of current flow through the same metal- 
to-semiconductor contact when an insulating 
layer separates the metal and the semiconductor 
and the field-effect determines the surface bar - 
rier under the metal. It is shown here that 
strong rectification, whose direction depends 
only on the bulk type, may occur for such con- 
tacts to extrinsic, but not intrinsic, semicon- 
ductors. Thus, rectification, without injection, 
may occur at the metal-to-semiconductor con- 
tact for the two-carrier system. 


DIFFUSION OF CADMIUM AND ZINC IN GAL- 
LIUM ARSENIDE. Bernard Goldstein, Radio 
Corporation of America Laboratories, Princeton, 
New Jersey (Received December 4, 1959). 


The diffusion of Cd and Zn in GaAs has been 
studied by using radioactive isotopes of these 
elements as tracers. The diffusion of Cd fol- 
lows the correct solution to the diffusion equation 
and its temperature dependence is of the cus- 
tomary form, D=D,exp(-E/kT), where E, the 
activation energy, is 2.43 ev and D, is 0.05 cm?/ 
sec. The diffusion of Zn from the vapor cannot, 
however, be described in terms of a single dif- 
fusion constant. The penetration curves de- 
crease much more sharply than they theoreti- 
cally should. Hall measurements indicate that 
all the Zn is substitutional and that it forms an 
impurity conduction band merging with the 
valence band. When Zn diffuses from a thin 
electroplated layer of radio-zinc, then the pene- 
tration profiles do correspond to the proper solu- 
tion to the diffusion equation. The diffusion con- 
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stants so determined have the usual temperaty, 
dependence given by D=D, exp(-E/kT), when 
D, is 15 cm?/sec and E is the same as that foun 
for Cd. From the work reported here and that 
of others, it is suggested that the diffusion of 
Cd and Zn in GaAs proceeds via vacancy migra. 
tion within the gallium sublattice. 


SURFACE TRANSPORT IN SEMICONDUCTORS, 
R. F. Greene, United States Naval Ordnance 
Laboratory, White Oak, Silver Spring, Mary- 
land, D. R. Frankl, Central Research Labora- 
tories, Sylvania Electric Products, Incorporated 
Bayside, New York, and Jay Zemel, United 
States Naval Ordnance Laboratory, White Oak, 
Silver Spring, Maryland (Received November 
13, 1959). 


A transport theory is given for electrons and 
holes in space-charge layers at semiconductor 
surfaces. For diffuse surface scattering, the 
effective surface mobilities may differ signifi- 
cantly from the bulk mobility for any strength of 
space-charge layer. Agreement with Schrieffer’s 
formulas is found only for strong space-charge 
layers, and the discrepancy is explained. The 
results are extended to cover an arbitrary de- 
gree of diffuseness of surface scattering and to 
cover samples of small thickness. 


FAR INFRARED TRANSMISSION THROUGH 
SUPERCONDUCTING FILMS. D. M. Ginsberg 
and M. Tinkham, Department of Physics, Uni- 
versity of California, Berkeley, California 
(Received December 30, 1959). 


The far infrared transmission through films 
of superconducting and normal lead, tin, indium, 
and mercury has been measured in the wave- 
length region between 0.1 and 1.1 mm. The 
transmission data have been analyzed to find 
the ratio of the complex conductivity in the 
superconducting state to that in the normal state, 
as a function of frequency. The width of the 
energy gap at 0°K may be estimated from the 
frequency of the extrapolated cutoff of the real 
part, o,(w), of the superconducting conductivity. 
The values so obtained are 4.00.5, 3.30.2, 
and 3.9+0.3 kT, for lead, tin, and indium, re- 
spectively. These values are in good agreement 
with those obtained in other experiments on bulk 
samples. The frequency dependence of 9,(w)/y 
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ig in qualitative agreement with the results of 

a calculation by Mattis and Bardeen based on 

the theory of Bardeen, Cooper, and Schrieffer, 
except for an unexpected hump in 0,(w) for lead 
and (tentatively) mercury at low frequencies. 
This hump may be due to the production of col- 
lective excitations or an anisotropy in the energy 
gap. It has also been found that a magnetic field 
as high as 8000 gauss applied in the plane of a 
lead film about 12 A thick has only a very small 
effect on the electromagnetic properties of the 
film. This is not surprising, in view of the re- 
sults of the microwave experiments of Pippard 
and of Spiewak. 


ATOMIC HEAT OF SODIUM METAL FROM 0.4 

TO 2°K. R. E. Gaumer and O. V. Heer, Depart- 
ment of Physics, Ohio State University, Columbus, 
Ohio (Received November 9, 1959). 


The atomic heat of sodium metal was measured 
from 0.4 to 2°K using the magnetic refrigerator 
calorimeter. The experimental data are given by 
C=1.32T +0.4857° millijoule/mole-°K. No anom- 
aly is observed in this temperature range. The 
experimental value of the electronic specific 
heat and the Pauli spin paramagnetism measured 
by Schumacher and Slichter are used for com- 
parison with the electronic band theory of sodium 
and those modifications introduced by electron- 
electron and lattice-electron interactions. 


TWO-PHONON INDIRECT TRANSITIONS AND 
LATTICE SCATTERING IN Si. W. P. Dumke, 
International Business Machines Research 
Laboratory, Poughkeepsie, New York (Received 
December 30, 1959). 


The probability of indirect transitions with the 
emission of two phonons has been calculated. 
Several of the intensity maxima in the intrinsic 
low-temperature emission spectrum of Si are 
explainable in terms of these transitions. Scat- 
tering matrix elements obtained from an analy- 
sis of the observed emission spectrum indicate 
that intervalley scattering is the dominant scat- 
tering mechanism for electrons in Si, with 
0.023-ev and 0.046-ev longitudinal acoustic mode 
phonons umklapp scattering electrons between 
valleys on the same and on different crystal 
aes, respectively. The valleys are approxi- 
mately 82% of the way from the center to the 





edge of the Brillouin zone. Optical mode scat- 
tering in the valence band is largely responsible 
for the anomalous temperature dependence of 
the intrinsic hole mobility. 


VARLEY MECHANISM FOR DEFECT FORMA- 
TION IN ALKALI HALIDES. D. L. Dexter, Uni- 
versity of Rochester, Rochester, New York (Re- 
ceived December 2, 1959). 


The Varley mechanism is examined, accord- 
ing to which Frenkel defects are produced in the 
halogen sublattice of alkali halides subsequent 
to multiple ionization of the halide ions. Argu- 
ments are presented to show that the lifetime of 
a positive halogen ion against recapture of elec- 
trons is orders of magnitude smaller than the 
ejection time of the halogen, and thus that the 
Varley mechanism is inoperative. The argu- 
ments may not be applicable for inner shells 
alone, but experimental evidence is adduced to 
eliminate this case. 


SOME OPTICAL PROPERTIES OF CADMIUM 
TELLURIDE. P. W. Davis and T. S. Shilliday, 
Battelle Memorial Institute, Columbus, Ohio 
(Received December 28, 1959). 


Some optical properties of CdTe samples hav- 
ing a carrier concentration of the order of 10'° 
cm™~* were measured. The index of refraction 
in the infrared was determined to be 2.61. Val- 
ues of the absorption coefficient a as a function 
of frequency were obtained at various tempera- 
tures. By plotting a” and a? against frequency, 
it was found that at room temperature the mini- 
mum energy for direct transitions is 1.50 ev; 
that for indirect transitions is 1.44 ev. The 
temperature dependencies between 77°K and 
373°K are given by (1.66 -5.6x10~“T) ev for 
direct transitions and (1.56 -4.1x10~T) ev for 
indirect transitions. 


LOW-TEMPERATURE PROPERTIES OF H, V, 
AND F CENTERS IN KCl AND KBr. John Cape 
and Gilbert Jacobs, * Department of Physics, 
University of Illinois, Urbana, Illinois (Re- 
ceived December 9, 1959). 


Potassium chloride and potassium bromide 
crystals were exposed to x rays at 10°K. The 
optical absorption produced by this irradiation 
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and the changes in optical absorption produced 
by subsequent annealing at higher temperatures 
were measured. The temperatures at which 
changes in optical absorption occurred were 
correlated with the temperatures at which free 
electrical charge appeared, and thermolumines- 
cence was observed. The absorption band at 
345 mu, in KCl, which has been named the H 
band, is shown to result from the superposition 
of two or more bands, one of which is the ab- 
sorption band due to self-trapped holes. In 
KCl, the self-trapped hole band bleaches ther- 
mally at 43°K with a release of free electrical 
charge. H centers disappear at 56°K with a 
release of free charge. The optical absorption 
band of the H center is shown to have its maxi- 
mum at 335 mu. In KBr, the thermal release 
of free charge at 30°K is attributed to the disap- 
pearance of H centers. No charge burst was 
observed in KBr, which may be attributed to 
the destruction of self-trapped holes. 


*Now at University of Ghent, Ghent, Belgium. 


FREQUENCY FACTOR AND ENERGY DISTRIBU - 
TION OF SHALLOW TRAPS IN CADMIUM SUL- 
FIDE. James J. Brophy and Robert J. Robinson, 
Physics Division, Armour Research Foundation, 
Chicago, Illinois (Received December 10, 1959). 


Current noise and photoconductivity measure- 
ments taken under uniform 5200A illumination 
on CdS single crystals are used to derive the 
energy distribution and frequency factor of 
shallow traps in the range 0.3 to 0.6 electron 
volt below the conduction band for samples of 
different CuCl impurity content. Trap densities 
varying from 10 to 10'7 cm~ ev™ and total 
trap concentrations of 10°* cm™ with discrete 
levels at 0.36, 0.43, and 0.60 ev below the con- 
duction band are observed. In a moderately 
doped, good photosensitive crystal, the traps 
also have a continuous distribution in energy 
and all have the same frequency factor, 10" 
sec™', which suggests that the traps are struc- 
turally similar. The results imply that a photo- 
electron may experience several thousand re- 
trapping transitions on the average before re- 
combining. It is possible to account semiquanti- 
tatively for the 1/f noise spectrum observed in 
some crystals at high frequencies in terms of 
the near -exponential trap distributions and con- 
stant frequency factor derived from low-fre- 
quency noise measurements. 
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CUBIC FIELD SPLITTING OF D LEVELS IN 
METALS. Joseph Callaway, Department of 
Physics, University of Miami, Coral Gables, 
Florida, and D. M. Edwards, Department of 
Mathematics, Queen Mary College, University 
of London, London, England (Received Decem- 
ber 8, 1959). 


The splitting of the fivefold degeneracy of free. 
atom d electron states by nonspherical com- 
ponents of a crystalline field is calculated. The 
crystal potential employed is that of a lattice of 
positive point charges screened by a uniform 
distribution of electrons. The calculation is 
done to first order in the cubic field, using hy- 
drogenic electron wave functions. The triply 
degenerate d state is lowered with respect to 
the doubly degenerate one in both body-centered 
and face-centered cubic lattices. Numerical 
results are given for both lattices. Finally, 
analytic atomic wave functions are used to esti- 
mate the splitting in iron and copper at the ob- 
served lattice spacing. The crystal field split- 
ting of these levels is found to be much smaller 
than the overlap splitting as obtained in pre- 
vious calculations for both materials. 


NEUTRON DIFFRACTION STUDIES OF ANTI- 

FERROMAGNETISM IN CrF, AND CrCl,. J. W. 
Cable, M. K. Wilkinson, and E. O. Wollan, Oak 
Ridge National Laboratory, Oak Ridge, Tennes- 
see (Received December 11, 1959). 


Neutron diffraction observations have been 
made on powder samples of CrF, and CrCl, 
from 298° to 4.2°K. These materials, which 
have crystal structures similar to rutile but dis- 
torted from tetragonality, become antiferromag- 
netic at low temperatures with different mag- 
netic structures. For CrF, the magnetic unit 
cell is identical with the chemical cell and the 
moments at the corner sites are directed op- 
positely to those at the center of the cell. The 
magnetic unit cell of CrCl, requires doubling of 
the b and c axes of the orthorhombic chemical 
cell and the structure consists of ferromagnetic 
(011) planes with adjacent planes antiparallel. 
In neither case were the intensities compatible 
with a magnetic axis directed along a simple 
crystallographic direction. Unique magnetic 
axes could not be definitely established but the 
data suggest that they lie parallel to the longest 
Cr-F and Cr-Cl bonds. Néel temperatures of 
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*rsity GROUND-STATE ENERGY AND GREEN’S 

cem- | -yNCTION FOR REDUCED HAMILTONIAN FOR 
§UPERCONDUCTIVITY. J. Bardeen and G. Rick- 
of free. ayzen, Department of Physics, University of 
m- lllinois, Urbana, Illinois (Received December 

|. The | 14, 1959). 


ice of In their theory of superconductivity, Bardeen, 
TM | cooper, and Schrieffer made use of a reduced 
= Hamiltonian which included only scattering of 

g hy- girs of particles of opposite momentum and 
Ply} jin, It is shown that the solution they obtained 


t to bya variational method is correct to O(1/n) for 
ntered |. large system. The single-particle Green’s 
cal function is derived and used to calculate the 

y; interaction energy. 

) esti- 


2 ob- "Present address: Department of Physics, Univer- 


split- sity of Liverpool, Liverpool, England. 


aller 
e- 
IRON SERIES HARTREE-FOCK CALCULATIONS. 
R. E. Watson, Solid State and Molecular Theory 
Group, Massachusetts Institute of Technology, 
Cambridge, Massachusetts (Received November 
\TI- 10, 1959). 


J. Seventy-six Hartree-Fock calculations have 
Oak FP been completed for the iron series (Sc to Cu) 
nnes- } atoms and ions. All calculations are for the 
(3d)" type of configuration (i.e., no 4s electrons 
en are present). The results are discussed, but 
, due to lack of space, are not presented in full 
h here. Comparisons are made with the experi- 


xt dis- — Mental ionization and multiplet spectra. Agree- 
omag- — "ent is poor due to limitations in the Hartree- 


g- Fock formalism. The results are used in an 

nit effort to gain information concerning correlation 
the energies. 

p- 

The 


ng of DSPERSION RELATIONS IN ATOMIC SCAT- 
cal TERING PROBLEMS. E. Gerjuoy and N. A. 
netic § Krall, John Jay Hopkins Laboratory for Pure 
el. ad Applied Science, General Atomic, San 


ible Diego, California (Received December 17, 
e 1959), 








Dispersion relations appropriate to the scat- 
ring of electrons by hydrogen atoms are de- 
tueed, and applied to actual measurements in 


the 
gest 
of 






the 0- 10 ev energy range. Two such experi- 
ments exist, yielding quite different results. 
Dispersion relations indicate that only certain 
angular distributions at low energy, and certain 
scattering amplitudes at high energy are con- 
sistent with these low-energy total cross-sec- 
tion measurements; this suggests experiments 
which could be used as checks on the accuracy 
of the existing measurements. 


BETA-GAMMA ANGULAR CORRELATION 
MEASUREMENTS ON Au”. I. DIRECTIONAL 
AND CIRCULAR POLARIZATION CORRELA- 
TION. Rolf M. Steffen, Department of Physics, 
Purdue University, Lafayette, Indiana (Received 
September 2, 1959). 


The beta-gamma directional correlation of 
Au’ has been measured as a function of the en- 
ergy of the beta particles. The anisotropy fac- 
tor A,(W) in the beta-gamma directional corre- 
lation Wa (6,W) =1+A,(W)P,(cosé) is propor- 
tional to Be / W and its value near the maximum 
beta energy is A,(2.8)=0.029+0.001. The beta 
directional gamma circular polarization corre- 
lation of Au’®® was measured using the back- 
scattering of the circularly polarized gamma ra- 
diation from a radially magnetized iron disk. 
The correlation is of the form: W, 1, W) 
= 1+ (0.454 0.07)(p/W)P,(cosé). The shape cor- 
rection factor of the first-forbidden beta spec- 
trum is independent of the beta energy above 
W=1.6 within 5%. 

The analysis of the data shows that the é-ap- 
proximation for first-forbidden nonunique beta 
transitions represents the Au’ results in a 
very satisfactory manner. 


BETA-GAMMA ANGULAR CORRELATION 
MEASUREMENTS ON Au”. Il. TRANSVERSE 
POLARIZATION OF THE BETA PARTICLES. 
P. C. Simms” and R. M. Steffen, Department of 
Physics, Purdue University, Lafayette, Indiana 
(Received October 15, 1959). 


As a consequence of the nonconservation of 
parity in beta decay, beta particles emitted in 
first-forbidden beta transitions exhibit a small 
degree of polarization transverse to their mo- 
mentum. The direction of the transverse polar- 
ization is defined with respect to a plane which 
is introduced by observing the direction of emis- 
sion of the beta particles and the direction of 
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emission of a gamma ray following the beta 
transition. The degree of the transverse polar- 
ization parallel to the beta-gamma plane Pr \ 
and the degree of polarization perpendicular to 
the beta-gamma plane P7, have been measured 
for Au’®*. The beta polarization has been de- 
tected by means of the left-right asymmetry in 
a Mott scattering process for an average electron 
energy W=2.0 mc? and for an angle@= 135° be- 
tween the beta momentum and the gamma direc- 
tion. The results of the measurements, Pry 
=+0.011+0.005 and Py, = 0.03+0.008, agree 
satisfactorily with the values calculated on the 
basis of the — approximation from the anisotropy 
of the Au’ beta-gamma directional correlation. 


"Present address: Pupin Physics Laboratory, 
Columbia University, New York, New York. 


ELASTIC SCATTERING OF 5.6-Mev NEUTRONS 
FROM CARBON. J. E. Braley and C. W. Cook, 
Convair, San Diego, California (Received 
October 26, 1959). 


Energetic neutrons produced from the D(d, n)He® 
reaction have been scattered from a cylindrical 
carbon sample to study the angular distribution 
of elastically scattered 5.6-Mev neutrons. The 
differential elastic scattering cross sections for 
carbon have been obtained for angles in the range 
of 30 to 150 degrees. A thin plastic neutron-proton 
recoil detector was used in the measurements to 
permit discrimination against 4.4-Mev gamma 
rays from the carbon sample, and other gamma- 
ray backgrounds. 


RADIOACTIVE DECAY OF Lu’”. R. G. Wilson 
and M. L. Pool, Department of Physics and 
Astronomy, Ohio State University, Columbus, 
Ohio (Received October 5, 1959; revised manu- 
script received January 11, 1960). 


Ytterbium oxide enriched to 95.9% in the 172 
mass number was irradiated with 6-Mev protons. 
An activity decaying by electron capture with a 
half-life of 6.70+ 0.04 days was produced and its 
assignment to Lu’™ confirmed by the identifica - 
tion of the ytterbium K x ray and by comparison 
with the activities produced by similar proton 
irradiations of the other enriched isotopes of 
ytterbium. The 4.0-hour positron activity pre- 
viously assigned to Lu'” was not observed and 
is best attributed to an impurity. The observed 
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activity of Lu’” consists of the L and K x rays 
of ytterbium and gamma rays with energies of 
79, 91, 113, 182, 203, 270, 324, 373, 490, 527 
697, 809, 900, 912, 1093, 1402, and 1583 key, 
Because no positron radiation exists in the 
activity of Lu’”, the mode of decay is solely by 
electron capture to Yb'”. Gamma-gamma 
coincidence measurements, energy considera- 
tions, and consideration of the relative numbers 
of the radiations observed in the activity of Ly” 
have led to the assignment of energy levels at 
530 (6+), 1172 (3+), 1263 (4+), 1375 (5+), 1662 
(3-), (1699), and 2072 (4+) kev in Yb in addition 


to the previously known 78.7 (2+)- and 260.2 (44). 


kev levels. The positions of all of the observed 
radiations and some observed only in conversin 
electron measurements are shown in a propose 
energy level scheme for the decay of Lu”. 
Approximate branching ratios for the electron 
capture disintegrations of Lu'” are also shown 
in the level scheme. Few if any electron capture 
transitions of Lu'” occur to the ground and first 
excited states of Yb’. Of the two predicted 
spins for the ground state of Lu’”, 4- is more 
consistent with the proposed energy level scheme 


MULTIPLE SCATTERING OF POLARIZED 
ELECTRONS. M. K. Sundaresan, Department 
of Physics, Panjab University, Chandigarh, 
India (Received December 29, 1960). 


The correct theory of multiple scattering is 
used to evaluate the numerical magnitudes of the 
depolarization produced due to multiple scatter- 
ing in gold foil of thickness 1 mg/cm’ for po- 
larized electrons of velocities: v/c=0.6, 0.7, 
0.8,0.9. The depolarization effect is found to be 
extremely small. The correction due to mul- 
tiple scattering to the electrostatic rotation of 
spins is also computed. 


DECAY OF THE i,,,. STATE IN Pb”, R. Stock- 
endal, Nobel Institute of Physics, Stockholm, 
Sweden (Received December 21, 1959). 


The Ly, Ly, Ly, Mp My, My» Myy, Np Mm 
and Oy conversion lines of a 26.22+ 0.01 kev 


transition in Pb*® have been found from studies 
with electromagnetically separated samples of 
Bi?®. The transition, having M2 character, most 
probably takes place from a 4-msec i,,, state 
of 1014.0 kev to the 987.8-kev 9/2- state, thus 
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causing the latter state to appear metastable, 

as has earlier been reported. This suggestion 

is also supported by strong evidence for the 
existence of a 310.5-kev E3 transition between 

the isomeric state and the 703.3-kev state. The 
energy 1014.0 kev of the i,,, state coincides with 
atransition energy in Pb?® earlier reported, 

and is found to fit well into the energy systematics 
af the i, ~/s2 M4 transitions in the odd lead 


isotopes. 


ISOMERIC TRANSITION IN Pb”. David E. 
Alburger, Brookhaven National Laboratory, 
Upton, New York (Received December 30, 1959). 


A (26+1)-kev transition in Pb*® occurring in 
the electron capture decay of Bi?® has been 
identified from its L, M, and N internal con- 
version electrons measured in an intermediate- 
image beta-ray spectrometer. By using Nal 
scintillation detectors behind the source in the 
spectrometer, the M and N lines of the 26-kev 
transition are found to be not in coincidence 
with electron-capture K x rays but they are in 
coincidence with a principal gamma ray of 1 Mev 
and weak components of ~0.7 and ~0.3 Mev. 
Since the Bi*® gamma rays of 987.8 and 284.2 
kev and a fraction of the 703.3-kev gamma rays 
have been shown by Vegors and Heath to be de- 
layed with respect to electron capture K x rays 
with a half-life of 4.8 msec, the present coinci- 
dence results indicate that the delayed radiation 
is associated with the 26-kev transition origin- 
ating from an isomeric state in Pb*® at 1013.8 
kev. The 26-kev transition is probably a quad- 
rupole and a possible assignment is M2 if the 
spin-parity assignment of the 987.8-kev level 
is 9/2- and if the 1013.8-kev level is the 13/2+ 
state predicted at 1.1 Mev by Pryce. The pos- 
sibility that all or part of the known 1014.2-kev 
gamma radiation constitutes the “missing” M4 
transition in Pb*® is discussed. 


DECAY OF Sm***. R. E. Sund, R. G. Arns, and 
M. L. Wiedenbeck, Harrison M. Randall Labo- 
ratory of Physics, The University of Michigan, 


Ann Arbor, Michigan (Received November 25, 
1959). 


Sm'* was found to decay to Eu’® with a half- 
life of 21.9+0.2 minutes. Internal conversion 
electrons corresponding to the gamma rays in 


Eu’ were observed in a magnetic spectrometer. 
The gamma-ray spectrum was studied with a 
well crystal. The spectra of gamma rays coin- 
cident with the x ray, 104-kev, 142-kev, and 
246-kev gamma rays were observed. A number 
of new, weak transitions are proposed. Direc- 
tional correlation measurements were made on 
the 142-kev—104-kev cascade. Possible spin 
assignments are discussed. 


RANGE -ENERGY RELATIONS FOR PROTONS 
IN VARIOUS SUBSTANCES. R. M. Sternheimer, 
Brookhaven National Laboratory, Upton, New 
York (Received December 30, 1959). 


An expression is obtained for the range-energy 
relation R(T») for protons (Tp =proton kinetic 
energy) as a function of the mean excitation 
potential J which enters into the Bethe-Bloch 
formula for the ionization loss dE /dx. The ex- 
pression for R(Tp) is obtained by an interpola- 
tion of the previously calculated range-energy 
relations for Be, Al, Cu, and Pb. The resulting 
expression for R(T p) can be used for any sub- 
stance, provided an appropriate value of 7 is 
assumed. Values are also obtained for the 
quantity g=(I/R)(dR/dI) which gives the frac- 
tional change of A for a small variation of the 
excitation potential /. 


DIRECTIONAL CORRELATION OF THE GAM- 
MA RAYS IN W"®. G. D. Hickman* and M. L. 

Wiedenbeck, Harrison M. Randall Laboratory, 
University of Michigan, Ann Arbor, Michigan 

(Received December 14, 1959). 


Directional correlation measurements were 
performed on seven cascades involving the gam- 
ma rays in W*®. The angular momenta of the 
excited states and the character of the gamma 
transitions which were determined by the cor- 
relation measurements are in good agreement 
with the values as determined by previous in- 
vestigators. There is evidence for the exist- 
ence of three rotational bands with the possibil- 
ity of a fourth band in W*™. The nature of these 
rotational bands is explained in terms of the 
Bohr-Mottelson model for spheroidal nuclei. 
The levels are characterized by the quantum 
numbers (K, J, 7). After correction for x-ray 
interference the 1.222 Mev-0.068 Mev correla- 
tion results in a spin of 3 to the 1.290-Mev 
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level. The 0.068-Mev gamma ray is essentially 
pure dipole radiation. The results of the 1.231 
Mev-0.100 Mev correlation are in good agree- 
ment with spin 3 to the 1.331-Mev level with a 
quadrupole content of (2+ 0.5) % for the 1.231- 
Mev radiation. The 0.152 Mev-1.222 Mev cor- 
relation is only in agreement with an assignment 
of spin 3 to the 1.374-Mev level with a quad- 
rupole content of less than 0.5% for the 0.152- 
Mev transition. The 0.152 Mev-1.122 Mev cor- 
relation yields a quadrupole content of (3- 11% 
or 94 - 99%) for the 1.122-Mev transition. The 
0.222 Mev-1.231 Mev, 0.264 Mev-1.222 Mev 
correlations are in good agreement with the 
assignment of spin 4 to the 1.554-Mev level. 
However, spins of 3 or 5 for the 1.554-Mev 
level are also in agreement with the correlation 
data. 


7 
Present address: Knolls Atomic Power Laboratory, 
Schenectady, New York. 


PROTON -PROTON SCATTERING AT 25 Mev. 

T. H. Jeong, L. H. Johnston, and D. E. Young, * 
University of Minnesota, Minneapolis, Minnesota, 
and C. N. Waddell, University of Southern Cali- 
fornia, Los Angeles, California (Received 
December 28, 1959). 


The differential cross section for proton- 
proton scattering has been measured for 23 
center -of-mass angles from 10° to 90°, with 
+ 0.8% absolute probable error at angles greater 
than 14°. The incident proton energy was 25.63 
Mev lab. The 90° cross section is 18.59 milli- 
barns, and the interference minimum of 17.09 
mb occurs at 24° c.m. A set of phase shifts 
which fit the data are: *S,, 49.5°; °P,, 8.2°; 
SP, -4.2°; *P,, 2.0°; *D,, 0.62°. 

*Now at Midwestern Universities Research Asso- 
ciation, Madison, Wisconsin. 


PROTON POLARIZATION IN p-d SCATTERING. 
S. M. Shafroth,* R. A. Chalmers, and E. N. 
Strait, Northwestern University, Evanston, 
Illinois, and R. E. Segel, Aeronautical Research 
Laboratory, Dayton, Ohio (Received November 
13, 1959). 


The proton polarization resulting from proton- 
deuteron elastic scattering has been measured 
in a double-scattering experiment. The first 
scattering took place in helium, which served 
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as a polarizer, and the left-right asymmetry 
observed in a second scattering in deuterium, 
Spurious asymmetries were checked by sub- 
stituting xenon and, separately, helium for 
deuterium as the second scatterer. Measure- 
ments taken at a proton-deuteron scattering 
energy of E, =3.4 Mev, @=45° and 90°, and 
E,=3.74 Mev, 6=45° all yielded results con- 
sfstent with no polarization. From these data 
it is concluded that the proton polarization in 
p-d elastic scattering is <~10% in this energy 
region. 


“Present address: Naval Research Laboratory, 
Washington, D. C. 


NONELASTIC SCATTERING OF FAST NEU- 
TRONS. J. T. Prud’homme, P. L. Okhuysen,* 
and I. L. Morgan, Texas Nuclear Corporation, 
Austin, Texas (Received December 2, 1959). 


The relative angular distributions of neutrons 
inelastically scattered from iron, yttrium, 
zirconium, radiogenic lead (88% Pb), lead, 
and bismuth were measured for neutrons in the 
region from 3.7 to 4.7 Mev. The relative angu- 
lar distributions of the low-energy (0.5 to 4 Mer 
neutrons resulting from nonelastic scattering of 
15.2- Mev neutrons were also measured. In eaci 
case the distributions were found to be isotropic 
within experimental error (+ 15%), therefore 
supporting earlier evidence of compound nucleuw 
formation as the predominant interaction mech- 
anism. 


*Now at The University of Texas, Austin, Texas. 


DECAY OF 45-DAY Fe. R. L. Heath, C. W. 
Reich, and D. G. Proctor, Atomic Energy Divi- 
sion, Phillips Petroleum Company, Idaho Falls, 
Idaho (Received December 14, 1959). 


The gamma rays following the decay of 45-day 
Fe have been studied using the techniques of 
gamma-ray scintillation spectrometry, including 
y-y coincidence and y-y directional correlation 
measurements. In addition to the previously re 
ported 0.192-, 1.10-, and 1.29-Mev gamma rays, 
two additional gamma rays having energies of 
0.145 and 0.337 Mev were observed. These 
latter two arise from a state at 1.43 Mev in Co™. 
All gamma rays were observed to decay with 
half-life of 4545 days. Directional correlation 
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measurements on the 0.14-1.29 Mev and 0.19- 

1.10 Mev cascades were performed. The results 
of these measurements, together with the gamma- 
ray relative intensities, suggest an assignment 

of }- to the 1.43- Mev state in Co. 


SURVEY OF (p,d) REACTIONS AT 22 MEV. 
¢, D. Goodman and J. B. Ball, Oak Ridge Na- 
tional Laboratory, Oak Ridge, Tennessee (Re- 
ceived December 11, 1959). 


Energy spectra of deuterons from (p,d) reac- 
tions on medium and heavy weight elements 
were surveyed. The experimental method of 
particle identification is described. The spec- 
tra show gross structure indicative of strong 
selection rules. The gross structure can be 
correlated with nuclear shell structure, and the 
levels which are most strongly excited are 
those which have the same shell configurations 
as the target with one neutron missing. Angular 
distributions confirm the shell assignments. 
This leads to a picture of the reaction mechanism 
for (p,d) reactions in which the incoming proton 
interacts principally with a single neutron rather 
than with the nucleus as a whole. 


MOMENTS OF INERTIA OF EVEN-EVEN 

EARTH NUCLEI. J. J. Griffin and M. Rich, Los 
Alamos Scientific Laboratory, University of 
California, Los Alamos, New Mexico (Received 
November 16, 1959). 


Moments of inertia of even-even nuclei are 
computed using the Nilsson model for deformed 
nuclei and the moment formula derived from 
the superconductor theory of nuclei. Values for 
the energy gap and the deformation are obtained 
from appropriate experimental data. Good 
agreement is found between the computed and 
observed energies of the first excited states of 
wenty-six rare earth nuclei. This success 
lends strong support to the superconductor 
theory of the nucleus. 


CHARGE INDEPENDENCE IN THE REACTIONS 
+d~2°+He® AND p+d—7*+H?® at 450 Mev. 

4.V. Crewe,* B. Ledley,t E. Lillethun, S. M. 
Marcowitz, and C. Rey, The Enrico Fermi 
stitute for Nuclear Studies, The University 













of Chicago, Chicago, Illinois (Received Decem- 
ber 30, 1959). 

An experiment has been performed to measure 
the branching ratio in the production of He* and 
H® in p-d collisions at 450 Mev at 140° in the 
c.m. system. The result is 2.134 0.15 which 
is in agreement with the prediction of a ratio 
2 on the basis of charge independence alone. 
The production cross sections were found to be 
(do/dQ) 4,3 = 5.414 0.29 ub/sr and (do/d2)ys 
= 11.55+ 0.49 ub/sr. 


"Also at Argonne National Laboratory, Lemont, 
Illinois, and the Department of Physics, The University 
of Chicago, Chicago, Illinois. 

TNow at the University of Sydney, Sydney, Australia. 


SEARCH FOR AN ELECTRIC DIPOLE MOMENT 
STRUCTURE OF THE MUON. David Berley and 
George Gidal, Columbia University, New York, 
New York (Received December 8, 1959). 


A search was made for an electric dipole mo- 
ment in the muon with a sensitivity of 0.1% of a 
muon Compton wavelength times the electronic 
charge. The motivation for this investigation is 
provided by the interest in finding some property 
of the muon which would indicate a structure 
different from that of the electron, even though 
such a structure would violate both parity con- 
servation and time-reversal invariance. The 
muons pass through the fringe field of the cyclo- 
tron and an additional system of magnets pro- 
ducing an electric field in their rest frame. Any 
electric dipole moment would precess about this 
field, producing a vertical plane component of 
spin transverse to the momentum. The absence 
of such a component within the stated sensitivity 
gives an upper limit to the electric dipole mo- 
ment of the muon as 2107** cm x the charge of 
the electron. 


NULL ELECTROMAGNETIC FIELDS IN GEN- 
ERAL RELATIVITY THEORY. Asher Peres, 
Department of Physics, Israel Institute of Tech- 
nology, Haifa, Israel (Received December 17, 
1959). 

The existence of null electromagnetic fields, 
within the frame of the general relativity theory, 
has been subject to some doubt. We here con- 
struct and investigate two solutions of the 
Maxwell- Einstein equations, representing null 
electromagnetic fields, and corresponding to 
unidirectional radiation flow. One of these fields 
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